UNCLASSIFIED 

>2 


OKLAHOMA  STATE  UNIV  STILLWATER  DEPT  OF  PHYSICS  F/9  7/4 

CHARACTERIZATION  Of  THE  OPTICAL  PROPERTIES  OF  LASER  AND  PMOSPHO— ETCOJl 
JAN  S2  R  C  POWELL  0AA929-B0-C-00S2 

AR0-161S2.1I-P  NL 


DUC  Fl£  COW  ADA110886 


! 

L 


REPORT  DOCUMENTATION  PAGE 

’  READ  WSTKUCnW'S 

BEFORE  COMPLETIHG~FORM 

t.  REPORT  NUMBER 

FINAL  REPORT  ARO  16182.1 1-1 

2.  GOVT  ACCESSION  NO. 

d  (  ■  '1  //  \  ^ 

I  RECIPIENT'S  CATALOG  NUMBER 

*.  TITLE  (an*  SunttUa) 

Characterization  of  the  Optical  Properties  of 
Laser  and  Phosphor  Crystals 


r.  author*  •> 

Richard  C.  Powell 

I.  PERFORMING  ORGANIZATION  NAME  AMO  AOORCSS 

Oklahoma  State  University 

Physics  Department 

Stillwater,  Oklahoma  74078 _ 

I.  CONTROLLING  OFFICE  MAKE  ANO  AOORESS 

U.S.  Army  Research  Office 

Post  Office  Box  12211 

Research  Triangle  Park,  NC  27709 

4.  MONITORING  AGENCY  NAME  *  AOORE SSfll  SSnilii  I 


Final  Report  | 

Jan.  1  ,  1979  -  Jan.  14,  1982 

1.  PERFORMING  ORG.  REPORT  NUMBER  j 

I 

a.  coh tract  on  grant  NUMaear*; 

DAAG29  79  G  0048: 

DAAG  29-80-C-0052 


a  Hoc  a  am  c 
AREA  «  VO  I 


.CtfCMT.  PftOJCCT,  T A 
X  UNIT  NUMSCRS 


ilUn|  Otilmm) 


12.  report  oate 

January  14,  1982 

)L  NUMBER  OP  PAGES 

112 _ _ 

lk  *4<iuRI+Y  CLASS.  *•*  '/.  I 


U  OfCLAttlPieA'hON/OOWNGftAOlN< 
SCHEDULE 


l«.  0ISTRI8UT1ON  STATEMEWT  (at  UUa  Rapart) 

.  .A  ppi  o\  o<.  lot'  j > 1 1 •  < * i o  ;  ili.'Strilmtion' 

unlimitcil. 


I  17.  OISTRIBUTIOH  STATEMENT  (at  i ft.  akaaramt  MlwW  In  Blank  39,  It  i 


IB.  SUPPLEMENTARY  NOTES 

THE  VIEW,  OPINIONS,  AND/OR  FINDINGS  CONTAINED  IN  THIS  REPORT 
ARE  THOSE  OF  THE  AUTHOR(S)  AMO  SHOULD  NOT  BE  CONSTRUED  AS 
AN  OFFICIAL  DEPARTMENT  OF  THE  AO'.  Y  "  JE’TIO.Y.  POLICY,  OR  DE¬ 
CISION,  UNLESS  SO  DESIGNATED  LY  OTHER  DOCUMENTATION. 

19.  KEY  WORDS  (Cantina*  on  rawaraa  at4*  It  naa***arr  and  Idantlty  *7  kl*ak  mmkar) 

LASERS 
RARE  EARTHS 

EXCITON  v 


r~X\0 

/llECTEI 


This  report  covers  the  work  done  in  the  Department  of  Physics  of  Oklahoma 
State  University  under  contract  number  DAAG  29-80-C-0052  supported  by  the  U.S. 
Army  Research  Office,  Research  Triangle  Park,  North  Carolina  from  January  15, 

1979  through  January  14,  1982.  SFbefresearch  involves  the  utilization  of  laser 
optics  techniques  such  as  time-resolved  site-selection  spectroscopy,  photo¬ 
acoustic  spectrosopy,  and  four-wave  mixing  spectroscopy  to  characterize 
dynamical  optical  properties  such  as  energy  transfer  and  radiationless  relaxatior 


do 


-j  €" 

Lc  - 


EDITION  OP  •  NOV  «S  Ik  1 

il  p  o 


SECURITY  CLASSIFICATION  OF  THIS  PA SI/RMa  0<M 


/Processes  in  phosphor  and  laser  matjhals.  The  materials  investigated  in¬ 
cluded  stoichiometric  laser  materials  such  as  neodymium  pentaphosphate 
and  neodymium  aluminum  borate,  and  lightly  doped  laser  materials  such  as 
neodymium  doped  yttrium  vanadate,  garnets  and  glasses.  In  addition,  several 
theoretical  approaches  and  computer  modeling  techniques  were  developed 
for  analyzing  the  properties  of  energy  transfer  among  Ions  in  solids. 

Some  of  the  important  results  obtained  in  this  work  include  the  direct 
measurement  of  the  exciton  diffusion  coefficient  in  NdxLa1_xP50^  crystals, 
the  measurement  of  the  quantum  efficiency  of  Nd^+in  various  hosts, 
and  the  development  of  Monte  Carlo  modeling  techniques  for  describing 
energy  migration  on  a  random  lattice. 


,  -r ‘K  f? 

U;  2anr.:-ord  r- 

|  Juntifio.-.tion. _  ~ 

{  Hv _ 

joiitrilmtlon/  ~ 

f  ability'  Codes 

I  ^  -4vaii.  anr' /or 
P'*  7?  ipeci.-=i 


SECURITY  CLASSIFICATION  OP  this  FAOtrWfcm  Oat*  I ataraO 


OKLAHOMA  STATE  UNIVERSITY 


Department  of  Physics 


CHARACTERIZATION  OF  THE  OPTICAL  PROPERTIES  OF 
LASER  AND  PHOSPHOR  CRYSTALS 

Richard  C.  Powell,  Ph.D. 

Principal  Investigator 


FINAL  REPORT 

Contract  Number  DAAG  29-80-C-0052 
(ARO  Proposal  Number  DRXRO-PH-l6l82-P) 
U.S.  Army  Research  Office 
Research  Triangle  Park,  North  Carolina 


January  15,  1979  -  January  14,  1982 


ABSTRACT 


This  report  covers  the  work  done  in  the  Department  of  Physics  of 

Oklahoma  State  University  under  contract  number  DAAG  2Q-30-C-0052  supported 

by  the  U.S.  Army  Research  Office,  Research  Triangle  Park.  North  Carolina 

from  January  15,  1079  through  January  14,  1932.  The  research  involves  the 

utilization  of  laser  optics  techniques  such  as  time-resolved  site-selection 

spectroscopy,  photoacoustic  spectroscopy,  and  four-wave  mixing  spectroscopy 

to  characterize  dynamical  optical  properties  such  as  energy  transfer  and 

radiationless  relaxation  processes  in  phosphor  and  laser  materials.  The 

materials  investigated  included  stoichiometric  laser  materials  such  as 

neodymium  pentaphosphate  and  neodymium  aluminum  borate,  and  lightly  doped 

laser  materials  such  *3  neodymium  doped  yttrium  vanadate,  garnets  and 

glasses.  In  addition,  several  theoretical  approaches  and  computer  modeling 

techniques  were  developed  for  analyzing  the  properties  of  energy  transfer 

among  ions  in  solids.  Some  of  the  important  results  obtained  in  this  work 

include  the  direct  measurement  of  the  exciton  diffusion  coefficient  in 

\d  La,  P-0,,  crvstals,  the  measurement  of  the  ouantum  efficiency  of  Nd" ~ 

X  t-x  3  14 

in  various  hosts,  and  the  development  of  Monte  Carlo  modeling  techniques 
for  describing  energy  migration  on  a  random  lattice. 
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I.  INTRODUCTION 


The  purpose  of  c.iis  research  was  to  enhance  our  understanding  of  the 
dynamical  properties  of  optically  pumped  solids  having  potential  applications 
as  laser  and  phosphc-r  materials.  The  major  research  goals  stated  in  the 
original  proposal  have  all  essentially  been  accomplished.  The  important 
results  obtained  during  the  three  years  of  this  contract  are  briefly  outlined 
in  this  section  and  presented  in  detail  in  the  remainder  of  the  report. 

1.1  Summary  of  Research  Accomplishments 
The  research  work  performed  under  this  contract  can  be  divided  into 
four  thrust  areas  as  listed  in  Table  I.  The  first  of  these  is  Materials 
Characterization.  It  is  important  to  compile  as  much  information  as  possible 
about  the  properties  of  optically  active  materials  to  be  available  when 
specific  design  requirements  are  generated.  This  is  especially  important 
in  a  field  such  as  electro-optic  technology  which  is  advancing  rapidly  and 
finding  many  new  applications.  The  second  thrust  area  is  to  obtain  a  better 
understanding  of  the  fundamental  physical  processes  underlying  the  optical 
properties  of  active  materials.  Of  special  interest  are  processes  which 
lead  to  the  enhancement  or  quenching  of  the  optical  emission.  The  third 
thrust  area  is  the  development  of  new  theoretical  approaches  and  computer 
modeling  techniques  to  provide  greater  accuracy  in  interpreting  experimental 
data  obtained  on  energy  transfer  among  ions  in  solids.  This  has  become 
especially  important  with  the  availability  of  microcomputers  which  can  be 
easily  interfaced  to  laboratory  equipment  for  data  acquisition  and  analysis. 
The  ability  to  use  a  three-pronged  attack  on  a  problem  involving  experi¬ 
mental.  theoretical,  and  modeling  work  provides  a  powerful  method  for  in¬ 
vestigating  the  physical  properties  of  solids.  The  final  thrust  area  is 
the  development  of  new  experimental  techniques.  These  involve  using  the 
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Table  1 


CHARACTERIZATION  OF  MATERIALS 


Nd  Y,  P.0,, 
x  1-x  5  14 

NdxLal-xP5°14 
mai3(bo3)4 

YVO^:Nd3+ 
Vl5°12:Nd3* 

I3(A10.;0a0.5)5°l2:Kd 


3+ 


Nd3Ga5012 

Nd:Silicate  glass 

Nd  terminate  glass 

Nd: Phosphate  glass 

Nd:Pentaphosphate  glass 

Eu  Y,  P.O, , 
x  1-x  5  H 

Al901:Cr3' 


UNDERSTANDING  OF  PHYSICAL 
PHENOMENA 


Exciton  dynamics 
Ion-ion  interaction 

Radiationless  and  vibronic 
decay  processes 

Concentration  quenching 


Research  Thrust  Areas 


DEVELOPMENT  OF  THEORETICAL  ANfi 
MODELING  TECHNIQUES 


Monte  Carlo  methods  for  energy  migration 
on  a  random  lattice 

Fast  diffusion  energy  transfer  with  direct 
interaction 

Energy  transfer  on  a  discrete  lattice 

Theory  of  degenerate  four-wave  mixing  for 
diffusion  measurements 


DEVELOPMENT  OF  EXPERIMENTAL  TECHNIQUES 

Time-Resolved  Site-Selection  Spectroscopy 
Photoacoustic  Spectroscopy- 

High  Pressure  Spectroscopy 

Four-wave  mixing  spectroscopy 


special  properties  of  lasers  in  new  spectroscopy  methods. 

The  materials  characterization  effort  focused  on  NdJ+-doped  samples, 
especially  mixed  pentaphosphate  crystals,  mixed  garnet  crystals,  and  glasses. 
Other  samples  wh^ch  we  initially  anticipated  studying  were  not  investigated 
because  of  the  lack  of  time  and  the  greater  interest  in  those  materials 
which  were  chosen  to  study.  Three  results  of  major  interest  came  from  this 
work.  The  first  is  a  better  understanding  of  the  low  quantum  efficiency 

'y  , 

of  YAG:NdJ~  laser  crystals.  Our  phtoacoustic  work  described  in  Sec.  V.l 
verified  the  quantum  efficiency  of  controversy  in  the  literature.  This  result 
coupled  with  the  results  of  the  work  performed  on  our  previous  Army  contract 
1DAAG  2Q-76-G-OOQQ)  showing  that  efficient  energy  transfer  takes  place  among 
Nd0  ions  in  nonequivalent  crystal  field  sites  in  this  material  has  led 
Pr.  L.G.  PeShaser  of  Hughes  Research  Laboratories  to  perform  an  investigation 
of  the  effects  of  hvdroxal  ions  in  YAG  crystals  and  he  has  developed  a  model 
attributing  the  low  quantum  efficiency  to  energy  migration  and  transfer 
to  "killer  sites"  which  are  NdJ"  ions  near  to  hydroxyl  ions.  This  points 
the  way  to  new  crystal  growth  and  heat  treatment  methods  which  will  greatly 
improve  the  efficiency  of  YAG  lasers. 

The  second  result  of  major  importance  is  a  better  understanding  of 
the  concentration  quenching  mechanisms  in  stoichiometric  laser  materials 
such  as  NdP -0. , .  The  various  measurements  performed  for  the  characterization 
of  this  type  of  material  are  reported  in  Sec.  III.  The  final  results  indicate 
that  below  about  lOfS  Nd°  concentration,  ion-ion  cross-relaxation  causes 
the  fluorescence  quenching  whereas  at  higher  concentrations  the  quenching 
occurs  by  efficient  migration  to  killer  sites.  These  sites  are  not  associated 
with  hvdroxal  ions  and  approximately  10'-  of  them  may  be  surface  sites.  The 
exact  nature  of  these  defect  centers  has  not  yet  been  determined. 

The  third  important  result  involved  with  material  characterization 
is  described  in  Sec.  IV. 3.  where  the  properties  of  energy  transfer  among 


T  — -  '  '  . . 

ions  in  different  types  of  glass  hosts  are  discussed.  The  efficiency 
of  energy  transfer  is  found  to  vary  with  host  composition  as  the  spectral 
overlap  and  Judd-Ofelt  parameters  change.  Also  the  phonon-assisted  nature 
,  of  the  energy  transfer  is  found  to  be  very  sensitive  to  the  type  of  glass 

host  involved. 

One  of  the  major  goals  in  the  thrust  area  of  understanding  physical 
phenomena  was  to  elucidate  properties  of  energy  migration  in  heavily  doped 
materials.  Studies  of  mixed  neodymium  pentaphosphate  crystals  described 
in  Sec.  III. 3  show  that  at  room  temperature  migration  occurs  by  a  diffusive 
mechanism.  Spatial  migration  takes  place  without  spectral  transfer  and  the 
excitons  travel  over  distances  of  the  order  of  O.J  urn.  Time  did  not  permit 
extending  these  measurements  to  low  temperatures  to  attempt  to  observe 
coherent  migration  but  this  work  will  be  careied  out  in  the  near  future. 

For  other  rare-earth  pentaphosphates  such  as  EuJ  it  was  found  that  spectral 
energy  transfer  takes  place  by  a  single  step  process  at  both  high  and  low- 
temperatures  (see  Sec.  II.2). 

The  two  results  of  greatest  importance  in  the  are  of  theory  and  modelin 
techniques  involved  the  development  of  methods  for  using  microcomputers 
such  as  our  LSI-11  system  for  doing  complex  simulations  of  physical  situa¬ 
tions.  The  first  of  these  uses  a  Monte  Carlo  method  for  describing  energy 
migration  among  ions  on  a  random  lattice.  It  is  shown  in  Sec.  I I . 1  that 
explicitly  accounting  for  the  random  nature  of  the  distribution  of  ions 
in  this  way  can  result  in  values  for  energy  transfer  parameters  which  differ 
by  a  factor  of  two  from  those  obtained  from  analytical  theories  whi'ch  assume 
a  uniform  distribution  of  ions.  The  second  modeling  technique  involves 
explicitly  accounting  for  the  discrete  structure  of  the  lattice  in  describin 
single  step  energy  transfer.  The  work  described  in  Sec.  II. 2  show  that  for 
low  concentrations  of  active  ions  the  results  are  equivalent  to  those  found 
from  analytical  expressions  which  are  obtained  after  making  a  continuum 
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lattice  approximation.  However,  for  highly  concentrated  materials  there 
can  be  more  than  an  order  of  magnitude  difference  in  the  parameters  predicted 
by  the  two  approaches. 

Three  major  areas  were  developed  involving  new  experimental  techniques. 
The  first  of  these  was  the  use  of  photoacoustic  spectroscopy  to  measure 
quantum  efficiencies  of  laser  materials  as  described  in  Sec.  V.l.  Although 
it  was  shown  that  this  can  be  done,  the  technique  is  not  more  accurate  or 
easy  to  perform  than  other  methods  currently  in  use.  In  Secs.  III. 2  and 
V.2  spectroscopy  methods  using  high  pressure  diamond  cells  in  conjunction 
with  laser  excitation  sources  are  discussed.  Pressure  measurements  add 
a  new  parameter  to  vary  in  spectroscopic  studies  along  with  temperature, 
concentration  and  time.  The  ability  to  highly  focus  laser  light  makes  a 
laser  an  ideal  source  to  use  with  small  diamond  anvil  cells.  This  has 
greatly  enhanced  the  potential  importance  of  high  pressure  spectroscopy. 
Four-wave  mixing  spectroscopy  was  developed  to  establish  and  prove  transient 
gratings  in  doped  solids  as  described  in  Sec.  III. 3.  This  has  distinct 
advantages  over  other  coherent  transient  techniques  such  as  free  induction 
decay  since  the  mestastable  state  under  investigation  does  not  have  to  be 
pumped  directly.  This  is  a  powerful  tool  for  measuring  spatial  migration 
of  energy  without  spectral  diffusion  and  our  results  represent  the  first 
direct  measurement  in  the  exciton  diffusion  length  of  a  stoichiometric  laser 
material.  In  addition  to  the  development  of  specific  experimental  methods 
as  discussed  above,  much  time  was  spent  in  interfacing  our  LSI-11  micro¬ 
computer  system  to  our  laboratory  equipment.  This  has  greatly  extended 
the  capabilities  efficiency  with  which  data  can  be  acquired  and  analyzed. 

1.2  Publications  and  Personnel 

The  work  performed  during  the  three  years  of  this  contract  resulted 
in  ten  publications  and  twenty-five  unpublished  presentations  and  colloquia. 
These  are  listed  in  Table  II. 
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The  personnel  making  major  contributions  to  this  work  include  the 
principal  investigator,  Richard  C.  Powell,  professor  of  physics,  and  two 
graduate  research  assistants  C.M.  Lawson  and  D.K.  Sardar  who  received  their 
doctoral  degrees  during  this  time  period.  Their  theses  are  listed  in  Table 
II.  One  other  graduate  student,  J.K.  Tyminski,  as  well  as  three  under¬ 
graduate  students,  D.P.  Neikirk,  E.E.  Freed,  and  G.P.  Quarles,  have  been 
supported  by  this  contract.  Two  visiting  professors  in  the  Physics  Department 
also  worked  on  this  project,  R.C.  Chow  and  L.D.  Merkle.  Also  it  was  a 
pleasure  to  collaborate  with  J.G.  Gualtieri  of  the  Army  Night  Vision  and 
Electro-Optics  Laboratory,  W.K.  Zwicker  of  Philips  Research  Laboratories, 
and  M.J.  Weber  of  Lawrence  Livermore  National  Laboratory  on  certain  aspects 
of  this  project. 

Supplemental  support  for  some  parts  of  this  research  was  provided  by 
a  grant  from  the  Materials  Research  Division  of  the  National  Science 
Foundation . 
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II.  THEORY  AND  COMPUTER  MODELING 


This  section  is  divided  into  three  parts.  In  the  first  section  the  manu¬ 
script  describes  a  Monte  Carlo  technique  for  treating  energy  migration  on 
a  random  lattice,  and  a  new  theoretical  model  for  treating  energy  transfer 
in  the  case  in  which  weak  direct  sensitizer-activator  interaction  acts  as 
a  small  perturbation  on  energy  transfer  by  diffusion  among  sensitizers.  In 
the  second  section  the  model  outlines  the  treatment  of  energy  transfer  among 
ions  on  a  discrete  lattice  and  shows  the  importance  of  using  this  theory 
in  describing  energy  transfer  among  ions  in  heavily  doped  solids.  The  third 
section  describes  the  development  of  the  theory  of  four-wave  mixing  in  the 
configuration  which  is  useful  for  studying  exciton  dynamics  in  solids. 

Appendix  A  reproduces  some  theoretical  aspects  of  this  problem  that  are  impor¬ 
tant  in  developing  the  final  analytical  expressions. 
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PHYSICAL  REVIEW  B  VOLUME  21 ,  NUMBER  9  I  MAY  1 980 

Models  for  energy  (ransfer  in  solids 

H.  C.  Chow  and  Richard  C.  Powell 

Dtpurrim’ni  of  PAiwj.  Oklahoma  Sul  O’  Lou  W'/Wi  Shll*ulfi  Oklahoma  'Ll  I  'V 
(Received  16  October  I97y) 

Two  new  approaches  lor  describing  ihe  dynamics  of  energy  iranslcr  in  solids  have  been 
developed  and  are  described  here  The  first  is  a  method  lor  treating  the  case  in  which  weak 
direct  sensitizer-activator  interaction  acts  as  a  small  perturbation  on  energv  transfer  by  dilfustun 
among  sensitizers  A  technique  involving  a  propagator  expansion  and  the  Born  approximation 
in  the  interaction  strength  is  used  to  solve  this  problem  The  second  approach  is  j  Monte  Carlo 
technique  to  simulate  me  migration  of  energy  on  a  random  distribution  ol  sensitizers  The  pre¬ 
dictions  of  hoih  these  models  are  compared  to  experimental  results  and  to  the  predictions  of 
other  theoretical  models.  In  their  regions  of  validity  they  predict  significantly  different  results 
than  those  of  the  models  commonly  used  The  models  developed  here  are  applicable  to  many 
irnpoitant  materials  such  as  those  used  lor  rare-earth  lasetc. 


I.  INTRODUCTION 

The  transfer  of  electronic  excitation  energy 
between  tons  or  molecules  in  a  solid  has  been  the 
subject  of  many  .  esligalions  for  over  forty  years. 
Recently  there  has  been  increased  interest  in  charac¬ 
terizing  this  process  in  maierials  used  for  laser  and 
phosphor  applications.  Energy  iransfer  can  be  used 
io  increase  the  pumping  efficiency  of  the  active  ions 
or  molecules  in  these  maierials  but  it  can  also  cause 
decreased  fluorescence  emission  through  concentra¬ 
tion  quenching  interactions.  Despite  the  significant 
amount  of  interest  and  research  activity  in  this  area 
there  are  still  some  important  aspects  to  the  problem 
of  energy  iransfer  which  jre  not  well  characterized 
and  undersiood  The  work  described  here  treais  iwo 
aspects  ol  (his  problem  which  have  no!  been  satisfac¬ 
torily  accounted  for  in  previous  developments 

Energy  transfer  is  generally  treated  in  one  of  two 
limiting  cases  The  first  is  a  direct  iransfer  from  an 
excited  "sensitizer"  to  an  unexcited  "activator."  The 
theory  describing  this  "single-siep  process"  has  been 
developed  in  the  classic  papers  ol  Forster'  and 
Dexter'  with  problems  such  as  the  effect  of  random 
disorder  being  attacked  recently  1-5  The  second  case 
is  that  of  energy  transfer  to  activators  after  multistep 
diffusion  among  sensitizers.  I  he  description  of  this 
excum"  diffusion  energy  iransfer  was  first  develop¬ 
ed  by  Frenkel.'1  Trltfaj, 7  and  Forster.1  A  major  prob¬ 
lem  in  studying  diffusional  energy  transler  is  how  to 
separate  the  properties  of  diffusion  among  sensitizers 
and  ihe  properties  of  trapping  at  an  activator  site 
from  the  total  energy-transfer  properties  which  are 
measured  Several  theories  have  been  proposed  to 
account  for  simultaneous  sensitizer  energy  dtlfuston 
and  xcnsitizer-acnvaior  ifiieraciion®''u  but  ihe  linal 
solutions  of  each  of  these  theories  involve  assump¬ 
tions  which  limit  their  validity  to  specific  cases  One 


important  case  to  which  the  currently  available 
theories  do  not  apply  is  one  in  which  the  direct 
sensitizer-activator  interaction  is  smaller  than  the  dif¬ 
fusional  term  but  not  small  enough  io  be  negligible 
A  theoretical  descriplion  of  this  case  is  developed  in 
Sec.  II  of  this  paper  and  xeveral  cases  where  n  may 
apply  are  discussed 

A  second  problem  involving  multistep  migraiion  of 
energy  concerns  the  spatial  distribution  of  sensitizers 
All  ihe  standard  theories  assume  a  uniformly  distri¬ 
buted  lattice  of  sensitizers  so  that  the  random  walk  of 
the  exciton  can  be  described  by  a  single  average  hop¬ 
ping  nme.  This  picture  vhoulu  be  valid  for  "host- 
senxitized"  energy  transfer  or  other  situations  involv¬ 
ing  high  concentrations  of  sensitizers  but  n  is  certain¬ 
ly  a  poor  jpproximjnon  if  ihe  sensitizers  are  random¬ 
ly  distributed  impurities  wah  low  concentrations  In 
Sec.  Ill  of  this  paper  we  describe  a  Monte  Carlo  treat¬ 
ment  of  this  problem  and  compare  ihe  results  to  ihe 
predictions  of  existing  theories  and  experimental 
results 


II  EFFECTS  OF  DIRECT  SENS1TIZER-ACT1V ATOR 
INTERACTION  ON  DIFFUSIONAL 
ENERGY  TRANSFER 

In  this  section  we  consider  the  situation  in  which 
energy  transfer  includes  bolh  diffusion  among  sensi¬ 
tizers  and  single-step  resonant  iransfer  between 
sensitizer-activator  pairs.  Yokola  and  Tammoto® 
treated  the  limiting  case  in  which  the  one-step 
iransfer  is  dominant  and  the  diffusive  contribution  is 
a  small  periurbation.  They  developed  an  interpola¬ 
tion  scheme  for  the  lime  dependence  of  ihe  excited 
sensitizer  concentration  in  this  limn  On  ihe  other 
hand,  if  diffusion  is  ihe  more  important  of  ihe  two 
transfer  processes,  the  usual  treatment10,  !l  is  to  as- 
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sume  some  ilniie  (rapping  radius  lor  (he  activators 
ouiside  of  which  diffusion  lakes  place,  while  ncgleci- 
ing  any  direct,  single-step  sensitizer-activator  interac¬ 
tion  Here  we  will  consider  the  last-diffusion  regime 
and  specifically  examine  what  effects  the  single-step 
transfer  might  have  on  the  time  development  of  the 
excited  sensitizer  concentration  otherwise  depleted  by 
diffusion. 

The  equation  governing  the  excited  sensitizer  con¬ 
centration  «(  T.i)  including  both  diffusion  and 
single-step  transfer  is 

9n  (  f  ,i)/dl  -  -  !3n(T  .1)  +  DV2n(T.l) 

—  2)v(  F  —  R, )  n  (  7,t)  .  (i) 

Here  fi  is  the  intrinsic  decay  rate  of  the  sensitizers,  O 
is  the  diffusion  coefficient,  R,  is  the  position  vector 
for  a  given  activator,  and  u(F  —  R,)  is  the  interaction 
strength  for  a  given  sensitizer-activator  pair.  In  this 
development  the  interaction  is  taken  to  be  of  the 
common  electric  dipole-dipole  type,  and  thus 
1MF-R.7  varies  inversely  as  the  sixth  power  of 
|F  —  R, | .  The  solution  to  Eq.  (I )  must  eventually  be 
averaged  over  the  configuration  of  activator  distribu¬ 
tions,  which  is  a  difficult  task.  With  the  assumption 
of  a  uniform  activator  distribution,  however,  Yokota 
and  Tanimoio1  have  shown  that  the  problem  is  re¬ 
duced  to  averaging  the  solution  of  a  single-center 
problem,  in  which  an  activator  acting  as  an  absorber 
is  located  at  F*-  0.  For  this  case 

d/i  ( r./)/bl  “  -  d"  (r.r)  +  D  V;«  (r.l)  -  v(r)n(  r.i )  . 


where 


(2) 


Fermi  pseudopotentia!  method  to  find  the  lowest 
c.genvalue  and  eigenfunction,  k2  and  respec¬ 

tively,  which  according  to  Eq.  <3)  Jetermme  the 
asymptotic  behavior  of  the  excited  sensitizer  concen¬ 
tration.  In  this  approach  <j>4(r)  thus  determined  is 
identically  zero  for  a  finite  distance,  (iSrS  a.  the 
so-called  scattering  length. 

The  above  remarks  serve  to  justify  the  use  of  an 
approximate  potential  in  lieu  of  the  exact  dipole- 
dipole  interaction.  Namely,  we  let 

v(  r)  =  vo  <  /" )  +  u(r)  (Si 


in  whicn 


vat  r> 


|0.  r  >  a 
00 .  r  XS  u 


u  (  r) 


0.  r  s=  <7  . 

u/rt.  r  >  cl 


'6 ) 


Wnh  the  use  of  this  approximate  potential,  the  kinet¬ 
ic  Eq  (2)  is  solved  exactly  except  for  u(r)  which  is 
treated  in  the  Born  approximation  Such  a  solution 
clearly  corresponds  to  the  case  in  which  energy 
transfer  via  diffusive  migration  is  stronger  than 
transfer  by  single-step  electric  dipole-dipole  resonant 
transfer. 

The  method  of  solution  is  based  on  the  well-known 
propagator  expansion.1’  Upon  writing  n(r.t) 

l)(r,i).  it  is  seen  because  of  the  linearity  of  the 
operations  in  Eq.  (2)  that  the  temporal  and  spatial 
development  of  <i>(r.i)  is  governed  by 

il /(  r.l)  -  J  G(  './|ri),/a)i|/(  ro.r<|)  1  Iro  *?) 

in  which  the  propagator  or  Green’s  function 
G  (  r.t \r0.ia)  obeys  the  equation 

!  -d/bt  +DV}-v(r))G(r.l'ro.l,j) 


is  the  isotropic  Laplacian  operator.  The  general  solu¬ 
tion  to  Eq.  (2)  can  be  written  as 

11 1  r.i)  —  £exp(  -  (il  -  Dk2l)<bi,(r)  (3) 

t 

with  1  r)  obeying  the  eigenvalue  equation 

C7;  -r  k!- D~'v(rl\<t>k{r) -0  <4> 

Because  of  the  assumed  r~*  dependence  of  v(r).  no 
solution  to  Eq  (4)  can  be  found  that  is  regular  at 
r  -0  Nevertheless,  it  is  possible”  to  employ  the 


—  Sir  -r„)S(r -r„)  .  (8) 

where  5(x)  is  the  Dirac  8  function.  If  instead  of  the 
full  electric  dipole-dipole  interaction,  one  solves  for 
the  Green's  function  for  the  approximate  potential 
vo(r) 

I  -3/9/  +OV’-uo(r)]CQ(r,/|r0,/J 

-6(r  -r0)6(f  -/„)  (9) 

then  the  objective  enunciated  in  the  preceding  para¬ 
graph  is  accomplished  by  obtaining 


(  r.  / 1  31  I*u(  r.l)  7-  9>|(  r.l) 

-  J  O  of  r.l  \  r{).  I<)  I  ili(  /  n  Jo)  rfro  —  J  Oyt  r.l  \r  1 ,  / 1 )  u  (  r  \ )  Oi|(  r  |,/|!ri,./(jJi|il  rij.lo)  dr  1  dt  \  f/r^ 
We  determine  G„(  r.i  ir,)./„J  b>  the  method  of  image.  Upon  writing 
(/„<  r.:  \r.,.t., )  -  r~'  G  t  r  1 '  r /„ ) 


Ud) 
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Eg  (8)  reads 

!  -  hi  hi  *■  0  I  h' !  hr1 )  -  i-0(  r)  |G  (  r.l  | /■„./„)  -  r„ftt  r  -  )  fl(  /  -  i„) 

The  effect  of  !•„(/>  is  accounted  for  if  il  is  required  that  </<  •  i  r „  /„)  vanishes  lor  /  sf  </  This  rr.av  be  Jicim- 
phshed  b>  introducing  a  source  ol  strength  ru  at  '  —  r,,  and  an  image  source  at  <  =  hi  -  leading  :.i 

r  ,  it1  1  -  l„t  ,  .  .  . 

G'0t  r.r  r, - 7  icxpl  -  ( r  ■  •  ,r>il)<  ’  -  /..» ! 

•'  1 4  tt /J  t ;  -  i„ )  | 

-  expl  -t  r  -  ,u  -  2u)  ■ AD  l  '  -  / ,  l  ! : 

where  tfl  />  is  a  unit  Heaviside  step  function  Tor  uniform  initial  excitation  of  sensin/ers.  'a1  /  at  ' ,  - 

the  result  is 


i.i  r  ,  1  -  —  -  —  erf 


[  r  r  |  <  d/)/l  ■  *  j  j 

which  .s  identical  to  the  result'  describing  the  diffusive  migration  o!  excited  settsiti/ers  m  :ne  presence  m  a 
of  radius  a 

The  elfect  of  concurrent  single-step  electric  dipole-dipole  resonant  traitsler  from  sensm/er  to  acticamr  ,s  . 
tamed  in  the  next  term  in  the  Horn  approximation 


ill ■  I  r.  I )  —  J*  Jr  |  u  t  r  i  )K ’  •  r.  1 1  r .  1 
w  here 

.sir.'  r.l- - —•  f  - — ‘/r— r  ,  'exp!  -  l  /  -  / 1  4  0  1  ;  -  r  I  i  -  exp! 

t  4  tr  /> ) "'  J "  1  /  -  r)1'- 


•i ,  -l.i  I ■  4  O  '  ~  7  I  i 


I  s-  u  r,  u 

| - -  —  ori  - - 


l  t  4  O  r  > 1  ; 


I  urther  progtess  cjn  be  made  with  the  aid  of  the  following  integrals  'whi.h  are  derivable  from  the  results  m  the 
appendix  .d  Ref  13): 


(  - — — expi  -  W  4 /> •  i  -  r  )  ]  -  2\  /  expl  -  /}|:  4 /).»)  —  —  I  /i  eric  — - 

J  •  -  -  i  '  /)  j  i4/3/i 


J - - —  exp!  -  d  • 

J  ■  '  /  —  i  • 


d  ;  40"  -  r )  I  erf  - 
UO-i 


ten  /.  '  « ’  -  f  ,/, 


fI  j  '  A  Pit'  • 

t  expl  -  -  —1  >i  ■: 

! n  1 

-i--  .l! 

1  <  4/M  ■  1 

\  1  i 

'  r  -  .1  1  .  '  < 

1 :  -  ii~_ 

i  tOI'  ;  t ' 

4/)' 

4  /ti ;  r- 

1  77.01  ^  ■'■--a- 

A.- 

■il  f.'  ♦  *'•*  •  ••-«!• 

ntegraf  Then  1  u 

1  1 5  »  G.in  he  * rmcr1  is 

14-7  0 1  • 


4-  V  (  /  t  r  -  .  )  -  I  I  r  -  r  -  2 I  | 
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where 

F,  (  ' >  - 

FAx)  - 


—  -  1  -t-erl 


(4  Dr)'n 


rr//3)l/2!x|erl 


( 4D/1 1 


(  itD) 


1/2 


\-  f  t  r. 


(  r.  -  a)2 


4  Di 


FA  .x)  - 


—  4 1  x  I 


(  ir D) 


n  X.  -./(«/), A11  ^  e  Xlxl/Ur.|i,  dz  ' 


Rather  than  finding  (he  explicit  space  and  (ime 
dependence  of  the  excited  sensitizer  concentration 
which  entails  a  further  integration  in  Eq  (14),  it  is 
expedient  to  obtain  the  time-dependent  energy- 
transfer  rate  A  ( /)  which  is  also  accessible  experimen¬ 
tally  The  latter  is  related  to  the  total  amount  of  ex¬ 
citation  \(r)  via14 

V(/l  «  V</)e4'-;V0exp(- J\o')<//']  (17) 

It  is  also  customarily  treated  as  (he  quantity  that  ap¬ 
pears  in  the  rate  equation 


time  dependence  with  a  somewhat  different  interpre¬ 
tation  of  the  physical  parameters  involved.  The  last 
factor  in  Eq.  (20)  is  the  integrated  amount  of  excita¬ 
tions  weighted  with  the  probability  per  unit  time  of 
making  a  single-step  transfer  to  an  activator  This 
represents  the  rate  of  energy  transfer  by  a  single-step 
resonant  process.  The  presence  of  such  single-step 
transfer  has  the  consequence  of  diminishing  the 
amount  of  excitations  available  for  transfer  by  dif¬ 
fusive  migration.  This  is  accounted  for  by  the  second 
factor  in  Eq.  (20)  which  is  intrinsically  negative  by 
virtue  of  the  fact  that 


dS{i)/dl  ~-k{i)S0 


(18) 


It  is  clear  that  the  two  definitions  are  equivalent  only 
in  the  limn  that  (  A  (/')<//'  is  much  less  than  unity 
To  proceed  further,  Eq.  (2)  is  integrated  over  the 
volume  of  the  sample  li  and  the  first  identity  of 
Green  invoked  to  obtain 
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Comparison  between  Eqs  (18)  and  (19)  leads  (o  the 
expression  for  the  time-dependent  energy-transfer 
rate  for  the  case  of  S,  activators  per  unit  volume, 
correct  to  the  first  Born  approximation 
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The  three  factors  that  appear  in  Eq  <20)  have  rather 
obvious  physical  interpretations.  The  first  (actor  has 
the  explicit  form 

*.,(/)  -4rrDa,V.|l  +  a  (  I  (21) 

and  is  simply  the  rate  at  which  excitations  arriv,.  dif¬ 
fusively  at  the  surfaces  of  the  activator  traps.  It  has 
been  used  js  the  standard  result1'  14  for  the  time- 
dependent  energy  transfer  m  the  fait  diffusion  re¬ 
gime  The  generalized  random  walk  model'"  which 
includes  an  extended  trapping  -egnm  yields  'he  same 
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Thus  the  time-dependent  energy-transfer  rate  for 
the  case  considered  here  may  be  written  as 
k  ( /)  —  Ac(  t)  +  k  |  ( /)  wuh  k o  being  the  usual  diffusion 
expression  given  in  Eq.  (21)  and  k  i  accounting  for 
the  iwo  effects  of  direct  sensitizer-activaior  interac¬ 
tion  discussed  above.  For  electric  dipole-dipole  in¬ 
teraction  the  latter  term  becomes 
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At  this  point  it  is  necessary  lo  use  numerical  integra¬ 
tion  lo  extract  the  explicit  time  dependence  of  this 
additional  contribution  to  the  energy-transfer  rate 
However,  since  n  is  readily  seen  that  the  magnitude 
of  k  |(  /)  varies  from  4»r  V, n/t  3u! )  at  /  -  0  to 
-  4  7T.V3a/(  1  5iJJ)  as  z  —  °o,  a  measure  ol  the  optimal 
significance  of  ihe  effects  of  single-step  transfer  in 
the  last-diffusion  regime  may  conveniently  be  chosen 
as  the  ratio  A  ;  t  z  •*  0 )  /  Au(  r  —  oo  I  —  n  ( 3  Da1 ) .  Ta  ble 
l1'"'  lists  values  of  this  ratio  for  a  number  of  sys¬ 
tems  for  which  fast  diffusion  has  been  idemTiei 
It  should  he  pointed  out  that  in  the  absence  if  any 
other  theoretical  considerations,  the  values  of  n.  ./ 
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TABLE  I.  Estimated  importance  of  direct  scnsiti/er-atiivjior  interaction  on  energy  transfer  m 
the  last-diffusion  regime 


Sysiem 

<«  ( cm  “/see) 

U  <  cm2/ sec) 

u  ( cm ) 

-0) 

k  qW  '  ) 

rbjAljOu* 

12  x  10'” 

1.25  x  nr" 

2.1  x  I0~7 

1.65  x  102 

YuaAbo  |Ho0 osFj6 

1.8  x  lO'41 

2.2  x  if)-" 

2  03  x  I0-,C 

161 

Eu,  Cr. 

5  0  x  lO'1* 

60*  lO-10 

1  2  x  I0-1 

0.134 

Anihracene  in  fluorene4 

1  74  x  |(j-n 

3  4  x  K)-1 

3  69  x  1 0** 

9  2  *  I0J 

Feiracene  in  jnihracenee 

1  7*  x  |ir’: 

3  1  x  )()-' 

1  84  /  I0-6 

1  67  x  l(J's 

Tetracene  in  naphthalene* 

3  52  x  UT15 

3  1  x  )()-» 

1.37  *  It)*6 

1  0"  *  lo-4 

“Reference  16.  JReference  15 

"Reference  17  'Reference  14 

'Estimated  using  a  -0  6761  from  Ret  8 


and  D  listed  in  Table  I  have  been  determined  by  fit¬ 
ting  the  experimental  data  with  the  standard  diffusion 
theory  and  are  therefore  not  totally  reliable  for  estab¬ 
lishing  or  dismissing  the  importance  of  single-step 
transfer  in  the  fast-diffusion  regime  [Note  the  in¬ 
consistency  indicated  by  the  fact  that  the  parameter 
.a/I  2 Da*)  is  larger  than  unity  in  several  cases. 1  The 
present  calculation  hopefully  can  provide  some  cri¬ 
teria  for  assigning  the  values  of  these  physical  quanti¬ 
ties  and  for  assessing  the  significance  of  resonant 
transfer  either  by  using  an  independently  determined 
set  of  parameters  or  by  comparing  the  full  time 
dependence  of  the  energy-transfer  rate  to  the  experi¬ 
mental  results 


111  HOPPING  ON  A  RANDOM  LATTICE 

Another  phenomenological  approach  to  the  treat¬ 
ment  of  energy  transfer  which  also  enjoys  rather  fre¬ 
quent  use  is  based  on  the  physical  picture  of  excita¬ 
tions  hopping  among  sensitizers  and  from  sensitizers 
to  activators  In  the  limit  of  many  steps  the  results 
of  this  random-walk  approach  are  equivalent  to  those 
of  the  diffusion  approach.  The  basic  premises  of 
random-walk  treatments  will  now  be  briefly  re¬ 
viewed’  11  Let  p/r)  denote  the  probability  that  an 
excited  sensitizer  is  located  at  R,  at  time  i.  the  equa¬ 
tion  for  excitation  migration  is 

- -5><r,-r ,W'>  . 

ill  , 

The  meanings  of  the  other  symbols  are  the  same  as 
those  encountered  in  the  last  section  except  that  R, 
now  denotes  the  position  vector  for  either  in  activa¬ 
tor  or  >enxitizer  The  solution  to  Hq  '24)  must  be 
averaged  over  the  configuration  of  activator-sensitizer 
.tis'nhultous  and  the  re.ult  designated  by  .pi  .•>  may 


be  correlated  to,  for  example,  the  luminescence  emit¬ 
ted  by  the  sensitizer.  It  should  be  noted  that  Eq 
(24)  neglects  the  effects  of  back  transfer  from  activa¬ 
tor  to  sensitizer  ions.  This  is  justified  for  systems  in 
which  relaxation  processes  on  the  activator  take  the 
excitation  out  of  resonance  with  the  sensitizer  which 
is  true  for  many  cases  that  have  been  investigated 
For  systems  where  back  transfer  is  not  a  negligible 
process  an  additional  term  is  present  in  Eq.  (24) 

If  the  excitations  are  incapable  of  migration  among 
the  sensitizers  themselves,  a  is  possible  to  solve  Eq 
(24)  and  carry  out  an  ensemble  average  over  the  uni¬ 
form  distribution  of  activators  and  the  result,  ji(  :l 
for  electric  dipole-dipole  interaction  is1 

p(  i) -exp  [  -  fit  -  (  C,/Cu)<  ir/Jr)  l/:|  ,  (25) 

where  C,  is  the  activator  concentration  and  the  "criti¬ 
cal  concentration"  is  Co  -  (  -  w R  ,J  ) Here  R.,  is  de¬ 
fined  as  the  distance  between  sensitizer  and  activator 
at  which  the  rate  of  energy  transfer  is  equal  lo  ihe  in¬ 
trinsic  decay  rate;  i.e.,  v(  R  )  -  j3(  Rn/R)<>  If,  howev¬ 
er,  the  excitations  may  hop  among  sensitizers,  the  in¬ 
teraction  strength  v(R)  that  appears  in  Eq  (24) 
abruptly  changes  each  time  a  hopping  takes  place  and 
is  therefore  a  random  variable.  By  assuming  that  the 
duration  over  which  the  interaction  strength  v(R) 
does  not  change  is  distributed  according  to 
ro1  exp  (  -  i/t0).  and  thus  identifying  r0  as  the  mean 
hopping  time,  Burshtein’  was  able  to  arrive  at  an 
equation  governing  4>(i) 

■f>(  f)  —  J/( /It’  °  +■  h  lit  r  )p(  i  -  /  )  i-  0  <// 

1 26  i 

The  solution  for  Eq.  ( 26 J  can  be  obtained  by  numeri¬ 
cal  methods1’  jnd  gives  results  equivalent  to  ihose  of 
the  Y  okota-Tanimota*  theory  :n  the  appropriate  limit 
On  t he  other  hand,  if  Eq  (24)  i-,  m  y>e  augmented 
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with  an  additional  term  to  account  for  the  back 
transfer  from  activator  to  sensitizer  ions,  the  problem 
is  significantly  more  complicated4  and  indeed  no  ana¬ 
lytic  solution  has  yet  been  found. 

The  merit  of  Burshtein's  approach  lies  in  its  simpli¬ 
city  anJ  us  relative  ease  for  numerical  solution.  It 
remains,  however,  an  approximation.  For  example, 
it  is  easily  seen  from  the  probability  distribution  for 
dipole-dipole  interaction  in  a  random  system20  that 
the  corresponding  hopping  lime  distribution  is  pro¬ 
portional  to  r_l/2exp(  —x/i.  where  A  is  a  constant  for 
a  given  concentration  of  sensitizers.  For  a  small  or 
moderately  large  concentration  of  sensitizers  this 
function  falls  off  much  slower  than  the  approxima¬ 
tion  used  by  Burshtem.  The  underlying  difficulty  in 
studying  hopping  motion  on  a  random  lattice  as 
described  above  is  this  possible  wide  distribution  of 
hopping  times  In  the  absence  of  better  theoretical 
techniques  to  handle  such  situations,  we  have  resort¬ 
ed  to  Monte  Carlo  methods  to  study  this  problem. 

As  will  be  seen,  our  study  leads  to  rather  different 
results  from  those  predicted  by  Eq.  (26) 

The  essence  of  the  Monte  Carlo  calculation  will 
now  be  described.  We  generate  a  finite  number  of 
excited  sensitizers,  allow  the  excitations  to  hop 
around  on  a  matrix  of  prespecified  concentrations  of 
sensitizer  and  activator  sues,  and  count  the  fractional 
excitations  that  survive  at  various  times  Intrinsic 
decay  of  the  excitations  is  easily  accounted  for  and 
thus  is  not  explicitly  considered  in  this  treatment. 

The  disappearance  of  excitations  therefore  occurs 
soiely  as  the  result  of  a  jump  onto  an  activator  sue 
which  has  a  jump  probability  dependent  upon  the 
fractional  occupation  of  sues  by  activators.  The  phys¬ 
ical  nature  of  the  random  distributions  of  the  sensi¬ 
tizers  and  activators  and  the  ion-ion  interaction 
mechanism  can  be  simulated  by  the  generation  of  a 
weighted  set  of  random  numbers  to  be  used  for  hop¬ 
ping  times.  The  standard  set  of  random  numbers 
generated  by  the  computer  are  uniformly  distributed 
and  must  be  transformed  to  have  the  desired  charac¬ 
teristics  For  example,  the  sensitizers  must  'bey  the 
law  of  distribution  of  the  nearest  neighbor  in  u  ran¬ 
dom  distribution  of  available  sues  (the  Hertzian  dis¬ 
tribution!  11  This  law,  m  the  representation  of  the 
number  of  available  sues,  v.  in  a  sphere  centered  on 
a  sensitizer  sue  and  interior  to  the  nearest- 
neighboring  sensitizer,  is  cexpi  —  it),  where  c  is  the 
fractional  sensitizer  occupancy  of  sues.  The  sequence 
of  numbers  having  such  a  feature  is  obtained  from  a 
sequence  of  uniformly  distributed  random  numbers  r 
\  ia  .he  expression  v  —  t  “ ‘  In  (  I  —  r)  (The  last  rela¬ 
tion  lollows  from  equating  the  cumulative  distribu¬ 
tions  of  the  two  sequences.)  In  a  similar  manner,  we 
require  that  the  sequence  of  numbers  reflects  the  na- 
'ure  of  an  electric  dipole-dipole  interaction  (which 
tails  'll  as  r~*  1  m  addition  to  that  of  a  random  distri¬ 
bution  ol  sensitizers  of  a  given  concentration.  The 


set  of  weighted  random  numbers  constructed  in  this 
way  are  then  used  as  hopping  times  for  the  simulated 
random  walk. 

A  computer  simulation  is  thus  constructed  which 
"observes"  the  hopping  of  ejeh  generated  excitation 
having  first  insured  that  the  excitation  was  not  creat¬ 
ed  on  an  activator  site.  The  time  lor  each  hop  is 
selected  from  the  weighted  set  ol  landom  numbers 
generated  by  the  method  described  above  This  al¬ 
lows  the  excited  sensitizer  on  each  step  of  the  ran¬ 
dom  walk  to  interact  with  any  of  the  other  randomly 
distributed  sensitizer  and  activator  ions  which  reflects 
the  electric  dipole-dipole  nature  ol  the  interaction  In 
this  manner  any  spatial  correlation  is  accounted  for 
Hupping  is  allowed  to  continue  until  cither  ol  the  fol¬ 
lowing  two  events  takes  piace:  an  activator  site  ,s  en¬ 
countered  or  a  particular  hop  takes  a  longer  time  than 
the  lime  of  interest  (in  the  present  study,  about  ten 
times  the  intrinsic  lifetime!.  The  survivjl  time  ol 
each  excitation  is  determined  by  calculating  the  sum 
of  all  its  previous  hopping  times  before  being  ter¬ 
minated  by  one  of  the  two  criteria  mentioned  previ¬ 
ously.  In  the  end,  a  bin  sort  is  performed  to  deter¬ 
mine  the  fractional  number  of  excited  sensitizers  that 
have  succeeded  in  avoiding  any  activator  sue  at  vari¬ 
ous  times.  To  ensure  that  correct  random-number 
sequences  are  used  in  the  calculation,  the  sequences 
are  tested  and  seen  to  obey  the  Hertzian  distribution 
Furthermore,  the  results  are  compared  to  test  wheth¬ 
er  they  are  insensitive  to  different  arbitrary  cutoffs  in 
the  random  numbers  used  in  ihe  calculation,  as  they 
should  be  if  the  proper  distribution  characteristics  are 
built  into  the  hopping  time  and  sensitizer  position 
distributions. 

The  results  of  some  typical  Monte  Carlo  calcula¬ 
tions  of  this  type  are  shown  in  Fig  1.  For  each  run 
an  average  hopping  time  is  determined  and  is  used  as 
an  input  in  Watts's  program14  for  integrating 
Burshtein's  equation.  The  predictions  of  the  Bursh- 
tein  approach,  alter  being  multiplied  by  exp  (/in.  are 
also  shown  in  Fig.  1  for  comparison.  In  general,  for 
the  same  interaction  strength,  the  Burshtem  approach 
predicts  a  faster  decrease  in  the  excited  sensitizer  po¬ 
pulation  than  does  the  Monte  Curio  simulation.  Note 
that  the  concentrations  used  in  curve  <a)  are  those  of 
a  typical  case  of  energy  transfer  among  Nd1*  ions  in 
yttrium  gjllium  garnet  crystals  which  has  been  inves¬ 
tigated  recently.21  The  data  shown  in  Fig.  7  of  Ref 
21  consist  of  the  time  evolution  of  the  intensity  ratios 
of  emission  centers  at  two  different  crystal-field  sites. 
Ions  in  one  type  of  sue  are  selectively  excited  by  the 
specific  laser  frequency  used,  and  energy  transfer  oc¬ 
curs  to  ions  in  the  other  type  of  sue  The  intensity 
ratios  lor  ihe  three  different  sets  of  spectral  transition 
shown  in  the  figure  are  fitted  with  the  same  theoreti¬ 
cal  parameters.  A  good  fit  to  experimental  data 
under  discussion  occurs  for  either  the  Monte  Carlo 
simulation  with  an  interaction  distance  of  20  A  or  ihe 
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FIG  I  Comparison  of  Monte  Carlo  predictions  ol' excita¬ 
tion  survival  probability.  with  ihe  predictions  of 

Hurshiein  s  iheory.  (a)  For  C  ,  -8  33  x  ID"3.  C  „  -  I  o’ 

<  ID--  d'1  -23  *  I  O'4  sec.  and  /< 0  —  20  A.  O  gives  the 

Monte  Carlo  predictions  and - gives  the  Uurshtem-lheory 

prediction.  gives  the  Burshiem  theoretical  prediction  for 
■he  same  concentrations  but  with  Ho”"  A.  lb)  For 
C,  -  3  33  x  I0'!.  C.-167XI0'3  H~' -2  5  x  I  O'4  sec.  and 

H,j  -!U  A.i  gives  Ihe  Monie  Carlo  predictions  and - is 

the  prediction  of  the  Burshtem  theory.  For  ihe  same  con- 
venirjtions  but  Hti  -  12  A  ihe  Burshtem-iheory  prediction  is 

given  by -  (c)  For  C,  -  8  33  x  I0'3  (.  „  - 1  7  x  10-3. 

j'1  -  2  5  •«  I0_4sec,  and  R0-20  A.  ::  represents  the  Monte 

Carlo  results  and  —  gives  ihe  Burshiem  results. - 

gives  ihe  prediction  of  ihe  Burshiem  iheory  lor  ihe  same 
concentrations  but  with  Rq  -  12  A 


Burshiem  theory  with  a  critical  interaction  distance  of 
il  A.  Theoretical  estimates  for  R o  lor  this  system 
are  closer  to  the  larger  value  hut  (his  must  be  con¬ 
sidered  us  only  a  rough  approximation  because  of  the 
complicated  njture  of  the  phonon-assisted  diffusion 
enhanced  energy  transfer  in  this  system.31 


IV  DISCUSSION  AND  CONCLUSIONS 

A  rigorous,  complete  solution  to  the  problem  of 
energy  transfer  among  randomly  distributed  sensitiz¬ 
ers  and  activators  is  still  lacking.  This  .s  especially 
•rue  :t  the  additional  complications  of  back  transfer 
■omi  activjtors  to  sensitizers  and  random  distribu¬ 
tions  of  transition  energies  are  included.  Recently 
there  have  appeared  a  significant  number  of  theoren- 
at  works  on  these  and  related  subjects.  Jmong  which 


may  be  mentioned  the  papers  of  Haan  and  Zwan/ig.' 
Huber,4  and  Holstein,  Lyo,  and  Orbach  3  These  gen 
erally  use  more  powerful  mathematical  techniques 
and  starting  with  the  elemental  interactions  among 
sensitizers  and  activators,  they  seek  to  elucidate  the 
connection  between  the  behavior  on  the  microscopic 
scale  and  the  somewhat  more  phenomenological  dif¬ 
fusion  and  random-walk  models.  The  final  exact 
solution  to  this  complicated  physical  problem  has  not 
yet  been  obtained.  In  the  preceding  sections  we  have 
concerned  ourselves  with  several  aspects  of  both  the 
diffusion  and  random-walk  models  of  the  energy 
probiem  with  the  objective  to  supplemem  and  clarify 
some  aspects  ol  these  existing  models  The  virtue  of 
these  models  has  been  their  capability  of  being 
brought  into  direct  comparison  with  experimental 
results  and  the  more  fundamental  theories4'  by  and 
large  confirm  in  the  appropriate  limits  the  validity  of 
these  phenomenological  models  This  provides  the 
justification  and  possible  usefulness  of  the  present 
work. 

We  have  developed  here  two  important  aspects  of 
the  energy-transfer  problem.  The  first  is  a  theoretical 
expression  for  treating  the  case:  The  direct  sensi- 
tizer-activator  interaction  is  a  small  perturbation  on 
the  transfer  of  energy  by  diffusion  among  the  sensi¬ 
tizers.  The  second  is  a  Monte  Carlo  approach  to 
simulate  the  migration  of  excitations  on  a  random 
distribution  of  sensitizers.  It  should  be  pointed  out 
that  the  approach  to  the  Monte  Carlo  simulation  used 
here  is  similar  to  the  continuous-time  random-walk 
model  developed  to  explain  anomalous  transit-time 
dispersion  for  charge  carriers  in  amorphous 
solids  33  33  It  attempts  to  account  for  the  possibility 
of  transfer  Irom  an  excited  sensitizer  to  any  other 
sensitizer  in  the  system  at  each  step  in  the  random 
walk  by  using  the  configuration-averaged  distribution 
of  hopping  times  at  each  site.  This  is  different  from 
the  Monte  Carlo  procedure  used  previously3  which 
generates  a  specific  lattice  topology,  assumes  only 
first  nearest-neighbor  steps  and  then  forms  a  config¬ 
uration  average  of  the  results.  In  that  study3  which 
correlates  the  decrease  in  fluorescence  emission  Irom 
ions  in  selective  excited  sites  due  to  energy  migration 
to  ions  in  different  types  of  neighboring  sues,  care 
must  be  exercised  in  Jealing  with  excitations  hopping 
back  to  the  original  site  since  this  offsets  the  time- 
dependent  line-narrowing  effect.  This  clearly 
presents  no  particular  problem  in  the  present  study,  a 
return  to  the  original  sue  has  no  distinguishable  ef¬ 
fects  from  hopping  to  a  fresh  sensitizer  sue 

Table  II  compares  the  commonly  used  phenomeno¬ 
logical  models  of  energy  transfer  and  indicates  the 
physical  situation  in  which  each  model  is  valid.  This 
shows  the  important  areas  of  usefulness  for  the  two 
theoretical  treatments  developed  here.  It  should  be 
stressed  that  each  of  these  models  should  be  applied 
only  lo  the  physical  situations  in  which  they  arc  valid. 
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TABLE  II  Phenomenological  models  lor  ihe  dynamics  of 
energy  transfer 


M  ode! 

Region  of  vipplica biluy 

Ref 

Forsier-Dexler 

Direct  i  <j  interaction; 
no  sensitizer  energy  diffusion 

1.  2 

>  okota-Tammota 

Sirong  s  -a  interaction, 
weak  sensmzer  energy  diffusion 

8,  1 

C  ho* -Powell 

Weak  \‘<j  interaction; 
strong  sensitizer  energy  diffusion 

This  paper 

Uurshtein 

Hquivalent  to  Yokota-Tammota 
model  and  diffusion/random-walk 
model  in  the  appropriate  limits 

9 

Soos- Powell 

Sirong  sensitizer  energy  diffusion 
w  uh  ihe  effects  of  extended 

trapping  regions 

10.  14 

0  illusion- 

No  direct  s-u  interaction; 

rjndom  walk 

energy  transfer  by  hopping  onto 
an  activator  site 

7,  1 1 

Munie  Carlo 

No  direct  s -u  interaction; 

3, 

transfer  by  energy  diffusion 
on  a  random  lattice 

this  paper 

ihis  has  noi  generally  been  true  in  the  past  For  ex¬ 
ample.  the  results  of  the  preceding  section  show  that, 
if  data  involving  energy  migration  on  a  random  distri¬ 
bution  of  sensitizers  are  analyzed  with  the  Burshtein 
model,  the  interaction  strength  one  obtains  is  approx¬ 
imately  a  factor  of  2  smaller  than  the  value  found 
from  the  Monte  Carlo  treatment  which  accounts  for 
the  random  distribution  of  hopping  times.  Note  that 
the  Burshtein  model  is  equivalent  to  the  Yokotj- 


Tammota  model  in  the  weak-diffuston  limit  and 
equivalent  to  the  standard  diffusion  result  in  the 
strong-diffusion  limit  However,  it  does  not  correctly 
account  for  the  effects  of  direct  sensitizer-activator 
interaction  when  it  is  not  negligible  in  the  last- 
diffusion  regime  and  it  does  not  correctly  account  for 
the  distribution  of  hopping  times  when  hopping  takes 
place  on  a  random  lattice.  The  models  developed 
here  are  most  applicable  m  these  regions  where  the 
Burshtein  approach  fails.  Also,  of  course,  there  are 
special  situations  where  none  of  ihe  models  in  Table 
II  apply.  An  example  of  this  is  when  high-resolution 
laser  experimental  techniques,  such  as  fluorescence 
line  narrowing  are  used  and  theories  accounting  lor 
differences  in  transition  energies  must  be  employed. 

In  summary,  we  have  developed  two  new  ap¬ 
proaches  to  the  problem  of  energy  transfer  of  solids 
which  allow  for  the  treatment  of  physical  situations 
involving  strong  diffusion  with  weak  direct  sensi- 
tizer-activalor  interaction  and  migration  of  energy  on 
a  random  distribution  of  sensitizers.  Both  ol  these 
situations  have  important  applications  to  impurity- 
doped  laser  materials  and  many  of  the  previous 
results  on  these  materials  should  be  reinterpreted  us¬ 
ing  the  new  models  described  here 
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Recent  results  are  presented  on  the 
the  characteristics  of  energy  transfer 
discrete  lattice  model  must  tie  used  to 


growth  of  rare  earth  pentaphosDha 
in  these  materials.  For  EuPcOiu  i 
account  for  the  observed  results. 
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I .  Introduction 


Turing  tr.e  past  two  years  we  have  been  studying  the  soectroscoDi  o  properties  of 
stoichiometric  laser  materials  and  espeoiallv  rare-earth  pen ta phosona tes .  Or.e  aspect  of 
these  materials  of  significant  importance  in  determining  the  concentration  Quenching 
properties  is  understanding  the  basic  physical  mechanism  of  energy  transfer  among  the  high 
concentrations  of  rare  earth  ions.  At  last  year's  Lasers  '  EO  Meeting  we  reported  on  the 
spatial  diffusion  of  energy  without  spectral  transfer  in  NdxLa i _xPij' • s  crystals.  Since  that 
time  the  major  portion  of  our  work  on  this  svsten  has  centered  around  the  development  of  the 
crystal  growth  of  higher  duality,  larger  sice  samples.  An  example  of  t n e  results  of  this 
-ffert  are  shown  in  Fig.  1  which  pictures  MdPgOj#  and  PrP^O*^  crystals  grown  in  the  :hilics 
L  at  d re t c r ie s  .  These  were  grown  tv  the  previously  described  flux  teo.rnicue  from  hot 
p.n  cstr.jrc  c  acid.  The  crystals  are  greater  tear.  A  cm  in  cross  section,  "he  sclibilitv  of 
-2  i  r  *■  ~  z-i*  ticrioscr.  ^*1 -?3  ap,  sr^s  \s  tin  *5  ins  izp.  ic  r  -  *a  ~  & 

-  i  r  i ; "  .  ;  v  increasing  ere  tencera  t  ure  wit  h  Cecrean  Lr.z  ionic  size,  crvscals  c  f 

ctner  hare  earth  pent aphos p rat a3  can  be  grown  pr.  similar  duality  and  sipe.  "pr.ee  past 
:  -  ■;  p  1  e  s  ..  i :  n  c  r  y  t  p  a  1  growth  nave  peer  an  important  limitation  to  the  '•ommeropal  cbolpoatidn 
f  stc lot pptie t r i c  laser  materials,  this  ability  tc  grow  suen  large  sice,  nign  cualitv 
crvstals  represents  a  significant  advance  in  tr.e  field. 


p  t  p  c  a  1 
p  e  c : ri 


-  ;  *  i  a  c 


pne  pasp  /ear  our  major  experimental  effort  has  centered  pn  the  study  of  the 
n. per- ties  of  mixed  EuPsGig  and  ?r?sO*n  mixed  crystals.  "nlike  the  Nd?gG-p  system, 
energy  diffusion  is  Title  prominent  in  each  of  these  types  of  samples.  Sesultp  o'- 
f  energy  transfer  m  mixed  EuxY  •  _  X?  =  Q  crystals  dene  in  collaboration  wits  J. 

are  presented  tele, 


:are  Earth  -an 


-x 


_  i-  .  r  =  c  •  u 

.(  -  A  ‘  -  ■<  ?  P  -  '  U 

v  x°p3-„ 


"  x  - 


no  ir^r.s:  er 

weak  sneo^ral  transfer 
3 1 r o r. s  noectral  t r a r. s ! 


ire 


»  r  a  ;  r 


-■jsi'.vej  3ite-2el.e2ti.crt  spectrcscccy  tech  nieces  were  us  el 
Ecrrl -u  and  mixed  EuxY • _X?=G < u  samples  containing  ’O'  and 
cradle  dye  laser  provides  excitation  pulses  of  about  5  ns 
half  widtn.  These  pulses  excite  the  sample  mounted  in  a 
a  1  .s  detected  with  a  1-mete-  nor.  och.remeter ,  a  cooled  pr.c 
e-  tc  deserve  tre  spectrum  at  a  specific  time  after  tme  1 
nde:  end  analyse:  wi'n  ar.  LSI-1’  computer  system. 


to  stjdy  energy  trar.s 
■t  Eu ’ *  ions.  A  nitn 
in  duration  and  less 
cryogenic  refrigerator 
t o t u b e ,  and  a  t  o  x  ;  a  -  c 
a3er  pulse.  The  ppeot 


-  -  •  Energy  Transfer  in.  Crystals. 

F i g  u-e  Z  snows  the  fluorescence  spectrum  at  \2  K  for  the  ’GO’  sample  at  short  times  and 
ir.ng  times  after  the  laser  pulse^  Here  we  show  an  expanded  region  of  one  of  the  sets  of 
pransitpens  between  the  ~Tq  and  ■ rm  levels.  The  two  peaks  in  the  spectra  represent  the  same 
transition  fer  pens  in  slightly  different  crystal  field  sites  .  The  narrow  laser  excitation 
line  selectively  excites  a  higher  percentage  of  the  ions  i.n  the  site  with  the  lower 
transition  energy  and  as  time  evolves  -r.is  energy  is  trar.ferred  to  ions  ir.  the  ether  type  of 
ppte.  Gi.milar  res. Its  are  cbt3ir.ed  at  nigh  temperatures  and  the  fluorescent  lifeti-es  'or 
t'.e  ’1GI  and  "0*  samples  change  only  slightly  as  a  function  of  temperature.  The  weak 
temperature  dependence  of  these  results  indicates  tnat  the  energy  transfer  characteristics 
are  associated  with  resonant  interactions  and  not  strongly  dependent  or.  rh.onpr  assisted 
processes. 
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appropriate  expression  for  the  energy  transfer  rate,  t  a  k  i  n  ?  into  account  me  d  i  s  c  r  ■? 
nature  of  the  lattice 

u  .  a  ,  _  .  '  <  ^  .  o  f 

y  -j  ■  a  '  3  ■*  •  J  - 1  ,  To  *< 


k  ;■  ;  a; 


i1'.  tJ-'lj^C. 


o  a 

!<  -  (  t.  /  r  )  v  n  /  ? 

:  dt  a  }e  3  :  4 

4  =  1  4 

The  summations  run  over  the  number  of  lattice  spheres  considered  ’<  and  uip£/  represents  the 
probability  of  occupancy  of  an  activator  ion  on  a  given  spnere.  The  expression  for  tms  is 


mere  r;  represents  tne  total  nutter  of  sites  on  tor.ere 
of  sites  o coup oed  tv  an  activator  ion  on  the  1  spnere. 
orptaoility  :f  site  occupancy  by  an  activator  .on. 

le  nave  developed  a  program  which  can  be  run  on  a  no 
.  union  'he  a  o  t  u  a  1  lattice  s  t  r  v  ;  t  j 
u  . ” "i 2 r  i ^ in  -s  •  5  .'  ■■'<? r  f  c  rn «?  J  .  T’n 

i : !  i  “  :  :  n  a  1  :  c  n  1  r  i  c  u  mm  *.  *>  -  n  e  r.  *?  r «  •/  1  r  a  r.  s  f  r  r  r’  2  *: 

:  m*  i  :m-? :  is  l*?ss  ‘-nr.  *  f»  :f  me  "ver~.ll  i  r  i  r.  3  f  r 

r.  1  f  7  -  t  r.e  e  n  ■ 

:  1  r.  v  1  n  e  •  j  1  t  .1 
1/  :cced  sir.p 
at  1  ice  Tio-je  1 . 
el  are  3  1  j n  i  f 
o :  e  1  are  are 


4  ,  1  r.  J  1  13  me  r  j n b e r 
x _  13  tne  frantic  nil 


? o n r  1: t e r  s  / s  t e r  cun',  as  an  _ H 1  -  ' 
e  n  *  e  r  e  1  ir.i  t  n  a  a  ~  ~  r  n  0  *’*  ’  •  * 
il  stances  r.  ter°  2 r. v 


was  r.  e  n  e  s  s  i  r  •/  t  n  « 


3  n  r  e  t  e  e  n  e  r  2  v  t  r  a  r.  s  f  e  r  "  i  1 0 
r.  1  7  3 1  i  *  n  C  1  ■/  mailer  t  r.  a  c 
r  tne  'Of'  s  a  r.  c  1  e  t  re  /  a  1  j 


.  1 .  .  1 1  e  v  e  1 


t  r.  e  theoretic.. 


?  :  t e c  ,  electric  licole  - dipole  ir. teraction.  J c w e v e r *  a  ■/ e r ’ 
-  n e a v  1 1  v  1  c c e c  3 2 r p  1  e s  it  is  r. e c e z s a r v  to  a c c c u r. t  fir  t - 0 
oe  and  usm?  the  normal  oor.ti.nuu.o  a : o rcxima t ion  lead.-  t  :  or: 
i  in'erootion  distance.  Two  other  attemcbs  have  been  made 
account L n?  for  the  discrete  .nature  of  the  lattice.  Stavel. 
;d  extressipns  for  luminescence  efficiency  but  net  the  er.er, 
tnlv  a  few  lattice  spheres  in  tneir  summation.  This  is  fa: 
v  o  on  car.  t  ra  t  e  d  materials  suer.  33  those  ir.ves  t  i  ea '  ed  -.ere  . 
mi  : 0 workers*  arid  they  derived  an  expression  fm  the  er. erv 
for  tneir  probability  of  .occupancy  they  sinpiy  used  t: 
Of  a  Site  Seine  an  activator,  wr.ioh  is  a  useful  aoorex  mat  • 

0  3  .on  as  tr.ose  investigated  here. 
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Theory  of  Four  Wave  Mixing 

The  four  wave  mixing  process  can  be  thought  cf  as  the  production  and 

reading  out  of  a  holographic  index  of  refraction  grating  in  a  non-linear 

medium.  In  the  four  wave  mixing  configuration  generally  used  to  study 

energy  migration,  a  laser  beam  is  split  into  two  strong  "pump''  beams  of 

wave  vectors  k  and  k,  ,  and  a  weaker  "Drobe"  beam  of  wavevector  k  which 
-a  "b  — 

counterpropagates  against  one  of  the  pump  beams.  The  two  pump  beams 
establish  or  "write”  the  grating  and  the  probe  beam  "reads"  the  grating 
(see  Figure  IV-1) . 

The  two  pump  beams  interfere  in  the  medium  and  optical  absorption  by 

the  active  ions  creates  a  spatial  distribution  of  excited  states  with  a 

sinusoidal  Dattern  of  wavevector  k  *  k.  -  k  .  Corresoondinq  to  the 

~g  ~b  -a  ' 

grating  wavevector  k  is  the  -grating  wavelength  A  given  by 


2  sin  6/2 

where  a  is  the  laser  wavelength  and  S  is  the  crossing  angle  of  the  two 
pump  beams . 

The  depth  of  the  grating  can  then  be  probed  by  Bragg  diffraction  of 

the  probe  beam  off  of  the  grating.  With  the  probe  beam  counterpropa- 

gating  against  the  second  write  beam  which  has  wavevector  k^,  the  Bragg 

condition  requires  the  Bragg  scattered  signal  beam  to  have  wave  vector 

k  =  k  +  k  =-k,  which  implies  that  the  signal  beam  counterpropagate 
-s  -p  -g  -a 

back  against  the  first  pump  beam. 

The  theory  of  FWM  has  recently  been  addressed  in  several  papers 
(26,40-43).  Two  fundamentally  different  approaches  have  been  used  to 
model  the  FWM  process. 

References  (40-42)  explicitly  consider  the  non-linear  wave  equations 


26 


Bragg  Condition  ks=kp+k 
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Figure  IV- 1.  Four  Wave  Mixing  Wavevector 
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where  the  electric  fields  are  coupled  by  the  non-linear  susceptibility  in 
the  material.  Furthermore,  (41)  and  (42)  consider  the  mechanisms 
creating  this  non-linear  susceptibility  by  modeling  the  system  as  an  en¬ 
semble  of  two  and  three  level  atoms,  respectively. 

References  (36,43)  model  the  system  in  a  very  different  way,  where 
the  probe  beam  Bragg  diffracts  off  of  a  sinusoidally  varying  complex 
index  of  refraction  grating.  In  this  type  of  development,  one  assumes  a 
spatially  varying  susceptibility  X(x,y,z)  which  forms  a  spatially  varying 
holographic  grating  in  the  material.  This  is  similar  to  the  approach 
ised  by  Kogelnik  (53).  The  emphasis  in  this  development  is-  m  under¬ 
standing  how  the  spatial  properties  (43)  or  temporal  properties  (36)  of 
the  grating  affect  the  Bragg  diffracted  signal.  However,  this  approach 
completely  ignores  the  mechanisms  creating  the  non-linear  susceptibility 
which  causes  the  spatially  varying  grating. 

The  approach  used  by  References  (41-42)  yields  important  information 
about  how  the  intrinsic  properties  of  the  material  affect  the  steady 
state  scattering  efficiency.  However,  References  (41-42)  assume  that 
the  pump  beams  are  counterpropagating  and  that  the  pump  beams  are  exactly 
phase  matched.  In  the  first  part  of  this  chapter,  the  FWM  scattering 
efficiency  will  be  derived  when  the  media  is  modeled  as  a  two-level  sys¬ 
tem  as  in  Reference  (41)  but  important  extensions  to  the  theory  will  be 
made.  Instead  of  assuming  counterpropagating  pump  beams,  the  assumption 
will  be  made  that  the  pump  beams  intersect  at  crossing  angle  6  ,a  more 
common  configuration  in  FWM  energy  migration  studies).  Furthermore , 
effects  arising  whenever  the  pump  beams  are  not  exactly  phase  matched 
will  be  explored.  Thus  the  first  part  of  this  chapter  will  give  infor¬ 
mation  ibout  how  the  pump  beam  properties  and  intrinsic  properties  of 
the  material  affect  the  steady  state  scattering  efficiency. 
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The  approach  used  in  Reference  (36)  is  much  more  useful  in  studying 

the  temporal  properties  of  the  grating  decay.  In  the  second  part  of  this 

chapter,  the  time  evolution  of  the  decay  of  the  signal  beam  intensity  is 

derived  in  a  manner  similar  to  Reference  (36). 

The  last  part  of  this  chapter  describes  FWM  data  taken  on 

Nd  La,  ?.0,  .  crystals, 
x  1-x  3  14 


Derivation  of  Scattering  Efficiency 


The  assumptions  will  be  made  that  all  beams  are  linearly  polarized 
in  the  same  direction,  with  the  pump  beam  electric  fields  given  by 
S.*r,t)  and  E,(r,t),  the  probe  beam  field  as  E  (r,t),  and  the  Brace  cif- 
fracted  signal  beam  field  as  2  (r,t). 

If  the  z  axis  is  taken  to  be  along  the  pump  beam  with  electric  field 
2.,  tnen  the  con  figuration  will  be  as  snown  in  Figure  17-2. 

If  one  makts  Che  "parametric  approximation"  that  the  pump  beams  ire 
undepie  ted  in  the  media,  then  the  four  electric  fields  are  given  as 

-u) t  -lk  .  •  r 

E,  (r,t)  =  A.  ( z)L  l  .."Jt  2  r  j 
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E  ,  i  r  ,  t !  =  A  i  l 
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Ir.  the  con  f iqura  t  ion 

used  , 

there  is  no 

incident  ;  icnai 

beam  i . e  . 

*  3)  and  tne  siunal 

•  ■) 

oeam 

Will  be  much 

weaker  than  the 

probe  bean 
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throughout  the  material.  Since  the  probe  beam  is  also  much  weaker  than 
the  oumo  beams. 


,  ,2  ,  .2  ,  .  ;  2 

I A_ ,  ,  , A  i  >>  I  '>  I A  ; 


<  p 

and  thus  j  —  |  <<  1  where  E  =  E2<r)  +  E  (r)  and  -E  =  E,  (r)  +  fr )  .  2 

“o' 

luit 

total  electric  tield  is  then  E(r,t)  =  l  (E  +AE)  . 

o 

The  wave  ecuation  that  these  fields  must  obey  m  the  material  is 


2  2  2 
7  S  -  s  m  3  E  =  u  3  ? 
o  o  t  at 


'  IV- 1) 


where  s  is  the  oermitivitv  constant,  u  is  the  oermeabilitv  constant, 
o*'  o 

and  ?  is  the  polarization.  The  polarization  ?  may  be  expressed  m  terms 
of  the  susceptibility  X  as 


?(E)  =  £  X<2>E 

o 


The  media  can  now  be  modeled  as  a  two-level  system.  As  is  she-/ 
Appendix  A,  the  susceptibility  for  a  two-level  system  is  given  by 


2a 

X<E)  -  -  -2  < 


1+0 


1  +  5"  +  ! E/E 


,  2 


where  5  is  the  normalized  detunina  from  line  center,  E  is  the  satur- 

s  1 

ation  intensify,  and  a  is  the  line  center  small- signal  field  attenuation 
'  o 

coefficient . 

IE 

Using  |— |  <<  1,  one  can  expand  tE)  and  P(E)  about  E^  to  first 


order  in  ; —  ■  to  on tain 
“o 


X (E  +AE) 
o 


X(E  )  U  - 


IE  IE*  +  E*AE) 
o  o 


1 


1  (l+o  )  +  :e 


„  i  ■ 


( I V-4 ) 
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The  polarization  is  given  by  (IV-2) ,  and  if  one  defines 

”>  *> 

E  [  “  ( 1+5*')  ,  then  to  first  order  in  |  AE/E  ;  / 

3  o 


i^t 

? (E  +  :1E)  =  i  z  X<E  )  IS  + 

o  0  o  o 


2  2 
(E  AE*  +  jE  I  AE) 

O  O' 


I  +  IE 


(17-5) 


t  2  2  2 

Since  =  lie.  1  =  Ik. I  ,  the  driving  term  y  3  ?  for  Eauation  (17-1) 

'1 1  ‘ -  3  o  t 

can  be  written  as 


.2 

j  3  ? 
o  t 


2  ,  2 

0  .  ( E  AE  *  +  i  E  j  IE) 

k~Z1Wt  X  (E  )  (E  +AE - 2 - - 

o  o  _  2 


The  left  hand  side  of  Equation  (IY-1)  now  must  be  calculated. 


?he  slowly  varying  envelope  approximation, 


id2Ai|  i,  dAi, 

I  «  tk  — i  i  =  1 

dz 


implies  that  the  field  amplitude  variation  due  to  nor.- linear  coup  line  is 
small  over  a  wavelength. 

.  ,  ,  ..ilkcosvz  ksir.cz) 

E.  and  E^  can  be  expressed  as  E,  =  A,  i2)t 


=  A 3  ( z)  i 


so  using  the  slowly  varying  envelope  approximation. 


iwt  _iV-  -iJv? 

E  -  e  y  3  E  =  2ikS.  (cos61  3  A  ‘.z!  *■  l  3  A.(z))  . 

o  O  t  z  1  z  : 


(17-7) 


lombming  ( 77-6)  and  ,17-7)  ,  the  wave  Equation  (IV-1.  Decones 


- ik . ■ r  -ik.-p 

2i  i cose 4.  1  3  A,  iz)  ♦  i  ■*  3  A.  (z)  )  *  -  kX(E 

z  1  z  '  o 
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One  mav  write  out  IE  ,  as 

o 


2  ‘  ■>  i(Ur!s^  •? 

!  A,  1  +  |A4!  +  A^i 


The  intensities  of  the  pump  beams  are  given  by 


A  i - ,  so  if  one  defines 


AC<ri<  )•-: 

=  A„A *1  1 


1  r  I  *• 

'  O  ’ 


Alone  the  z  axis,  ;  is  given  by  ;(z)  =  A, A  .  ■ 


oscillates  with  z  with  oscillation  wavelength 


>S  i-COS2 


has  an  interference  term  that  oscillates  with  o  with  oscillation 


wave .enc 


Tne  on^v  te 


rms  of  interest  on  the  right  hand  side  of  (17-3) 


these  that  satisfy  the  phase  matching  condition;  i.e.,  those  terms  that 


svnc.nronousiv  drive  .  E  -  z  u  j 
•  ■  o  o  t 


E  as  either 


Rewriting  ,'IV-3)  and  numbering  the  terms,  (IV-3)  becomes 


COS3?. 
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E  AE* 
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!eoI2  4e 


(I  +  |e  i2)2 

s  o' 


(I 


2  2' 


( IV-9) 


I'hJ 


term  (iii) 


term  (iv) 


Terms  (i) ,  (ii) ,  (iii) ,  and  (iv)  must  now  each  in  turn  be  examined 
to  determine  which  of  these  satisfy  the  phase  matching  condition. 


>ince  E  =  A, 2 
o 


term  !i)  is  clearly  net  synchronous 


with  the  left  hand  side  of  '.IV-9)  and  thus  term  (i)  is  not  phase  matened. 

-ik  •  r  - ik .  • r 

Since  IE  =  A  iz)(  A^(z)2  J  ,  term  (it)  clearly  is  phase 

matched  alt.nough  the  amplitude  of  this  tarn  will  be  modulated  by  tr.e 

oscillator1/  behavior  of  E  1  whicn  aooears  in  the  denominator . 

o' 

“* 

E *"  mv  be  written  exolicitly  as 
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■:  *  v  *  -a? 


i  .< ,  *  k  . •  r 


—<*A  rt  .  i. 
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and  IE  *  x.ay  be  written  out  exni  ;:i:iv  as  EE*  -  A.  - .  *. 


r\  .  •  -  /  A, 
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2 A  A-A.i  +  2A.A  A 
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Since  none  of  these  terms  are  synchronous  v/ith  the  lef.t  hand  side  of 
>,)  ,  term  (iii)  is  not  phase  matched. 

Term  iiv)  may  be  explicitly  written  out  as 
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E  j~  AE 
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i  (k  -k  ) • r  -ik  •  r 
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The  first  three  terms  in  this  expression  are  obviously  chase  match 

2  -2 

althouch  thev  are  modulated  by  (I  +  !E  ! 

s  ‘o' 

The  last  2  terms  in  square  brackets  must  be  looked  at  more  closely 
The  fourth  term  is 

-it*,-*-,)  *? 

(A*A  Z  "  )  -ik\ • r 

r  2  4  -i,  „  -I  ~ 

«.*  i*0:V  ' 1 

The  quantity  above  in  the  square  bracket  can  be  thought  of  as  mcdu 

-ik. 


iatmc  the  phase  matched  term  A,  i 

r 


■1 


hotice  that  since  the  exponent  in  the  square  bracket  above  has  the 
same  courier  comoonent  as  the  oscillatory  oart  of  'z  ,  this  term  will 

O' 

not  necessarily  average  to  zero. 

Similarly,  the  last  term  on  the  right  hand  side  of  the  expression 
far 


IE 


il  -  ;S  |2) 
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The  quantity  in  this  square  bracket  modulates  the  chase-  matched  term 
- lk ,  •  r 

A  i  J  ,  and  the  exponent  in  the  square  bracket  has  the  same 'four ter 

romoonent  as  tne  oscillatory  cortion  of  Z  ~ ,  so  tms  term,  aisc  will  net 

o 

necessarily  average  co  zero. 

Equation  (IV-9)  may  therefore  be  rewritten  as 


act  the  signal  beam  A_(z)  will  be  muen  weaker  than  the  prebe  beam 
,  and  since 
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v,  =  'A,U 
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The  phase  mismatch  of  the  two  pump  beams  u  is  -given  as  i,  =  ;n-5.  and 
,j|  and  (z)  become 


(2)  =  |A2A4j* 


i (k ( 1-cosc ) 2+1) 


(XV- 16) 


,*  ,  .  ■  ,  : „-i (k(l-cosd) z-1) 

(  2 )  —  ,  i  1 

1  h 


The  oscillatory  behavior  of  jE^j  can  be  shuwn  more  explicitly  by 


defining  the  parameters 


D  =  I  +  I  +  I 
o  s  2  4 


a  2  ,A2A4 ; 


Equations  (XV- 12)  and  ( IV— 13)  can  then  be  expressed  as 
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(XV- 13) 


These  eauations  can  be  simplified  somewhat  bv  defimnq  2  as 


P 
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3  =  k  { 1-cos vi  . 

Equations  (IV-18)  and  (17-19)  may  then  ce  expressed  as 

3  A,  (z)  +  ?  (z)  A,  (z)  =  0  (17-20) 

z  1  1 

and 

3  A,(z)  +  P iZ)A  iz)  =  2<z)  '17-21} 

z  3  3 


where 
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.17-24) 


The  solutions  to  (17-20)  and  (17-21)  fall  into  2  categories,  depending  on 
whether  or  not  ; (z)  and  ;*(z)  oscillate  very  rapidly  over  distances  where 
A,  z)  and  A^(z)  change  appreciably. 


Case  I  (9  is  Not  Small) 


is  not  small  then  A. (z)  and  A  (z!  do  not  change  appreciably 


over 
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as  1-cosci 


Since  for  Case  I  ;(z)  and  ;*(z)  oscillate  many  times  over  distances  where 
A  ( z)  and  A^(z)  change  significantly,  the  quantities  P(z),  P'zl  and  1‘z) 
can  be  averaged  over  an  oscillation  wavelength  a 
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o '  3  ,  o  s 


toss 


V-4aV/2'  ' 
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ie  complex  absorption  coefficient  i  can  be  defined  as 


cos0  <P  .  z) > 
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or 
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written  as 
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where  and  ;  are  given  explicitly  by  ( IV—  32)  and  (IV— 33) .  normally 
a  degenerate  FWM  experiment  the  two  pump  beams  and  the  probe  beam  are 
tbtamed  by  splitting  one  laser  beam  into  three  parts.  When  this  is  t 
case,  one  sees  from  '17-35)  that  the  output  Bragg  diffracted  signal  be 
intensity  will  vary  as  the  cube  of  the  laser  power. 

A  good  measure  of  the  "scattering  efficiency"  of  the  four  wave 
mixing  process  is  the  ratio  of  the  exiting  signal  beam  |A,(0) j  “  to  the 
exiting  probe  beam  in  the  absence  of  the  pump  beam  interactions. 

vising  (IV-25)  and  (IV-26)  ,  one  obtains 


lA1(0) 

o 

So  the  scattering  efficiency  n  is  given  by 


2, 
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-%  d+sec3)L 

cos [£  ( l-sec3 ) l] }  . 


The  dependence  of  che  scattering  efficiency  n  on  crossing  angle  3 
can  be  clarified  somewhat  whenever  j £L{l-sec5) •  is  small. 

This  will  usually  be  true  since  the  experiment  does  not  work  well 


ii)  j£L,l  is  large  since  if  |£l!  is  large  the  beams  are  essentially 
extinguished  in  the  crystal. 

(ii)  3  is  very  large  since  the  output  3raga  diffracted  signal  beam 
intensity  decreases  with  increasing  3. 

Thus  if  jlMl-secd)  '  is  small,  then  to  first  order  in  dbd-secl,  ;  , 
.17-36)  becomes 
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Equation  (IV-37)  shows  that  the  scattering  efficiency  n  decreases 
with  increasing  crossing  angle  4.* 

It  is  now  of  interest  to  calculate  the  scattering  efficiency  wr.er 
the  pump  intensity  is  well  below  saturation  intensity,  i..e.,  when 
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so  that  (IV-37)  becomes 
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( IV- 36) 


*In  thick  samples,  there  will  be  another  effect  that  will  decrease 
the  output  signal  intensity  with  increasing  crossing  ancle  t .  In  thick 
samples,  as  3  is  increased,  the  beam  overlap  volume  in  the  sample  may 
decrease,  which  is  not  taken  into  account  m  this  development. 
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To  lowest  order  in  - - - ,  (XV- 38)  becomes 
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'  IV- 3H ) 


weu  oe.ow 
saturation 


Experimentally,  one  usually  works  with  crossing  angles  large  sr.cucn 
that  Case  I  applies  and  Equation  ( I V — 3 6 )  ,  ;iV-37),  or  <IV-39)  will  be 

applicable.  However,  it  is  of  interest  tc  derive  the  scattering  effi¬ 
ciency  for  very  small  angles-  in  certain  limiting  cases  to  discover  what 
new  effects  should  be  expected  at  very  small  crossing  angles. 


Case  II  (3  is  Small; 


If  3  is  verv  small,  then  A, (z)  and  A_-,z;  do  chance  actreciablv  eve: 

i  u  —  -  -  - 

=  - — - —  and  thus  ? (z)  ,  P[z),  and  Q(z)  cannot  be  averacec  over  ' 

05  --COS"  "  OS 

Equations  (IV-20)  and  (IV-21)  cannot  be  solved  for  exactly  m  this 

case,  since  ?(z),  ? \z) ,  and  Q(z)  cannot  be  averaged. 

Equations  (IV-2C)  and  (IV-21)  will  be  approximately  so.ved  ftr  m 

the  case  where  the  pump  beams  are  well  below  saturation  intensity,  1.3., 
a 

tor  the  case  that  —  is  small. 

s 

One  can  now  define  a  as 
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ana  to  first  order  m  — ;  (IV-22)  becomes 
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P  (2) 


-  r{l  -  —  cos(oz+A)  -r  7—  sir.fiz+A); 


The  solution  to  Eauation  (IV-20)  is  then 
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Equation  (IV-41)  then  becomes 
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To  first  order  in  a/I_,  Equation  (IV-24)  becomes 
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and  since  I  A. 12)  '  <<  'A  z)  !,  ?(z)  may  be  approximated  as  ?• z) 
mat  Equation  '17-21)  may  be  approximated  as 
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w.iere  -•  given  by  ,17-42,'  and  2.  is  a  constant  :na:  remains  to  ce  i 

termined  using  boundary  conditions. 

To  perform  the  integration  indicated  m  Equation  117-46)  ,  one  r.ee 

to  expand  the  exponent  inside  the  integral  in  powers  of  a/Z  .  3ir.ce  - 

s 

factor  a/I_  in  the  exponent  is  multiplied  cv  j  factor  of  v  £ ,  it  is 

3 

important  to  briefly  examine  magnitude  of  a/c,  to  be  sure  that  it  is 
net  large.  This  parameter  can  be  rewritten  as 


uL 


-  Lr 

The  beams  are  approximately  attenuated  as  as  they  pas 

the  crystal,  so  j  iL|  must  be  of  order  unity  or  smaller  or  else 
are  so  severely  attenuated  through  the  crystal  tnat  the  Brace  i 
signal  beam  will  be  negligible. 

Furthermore,  ;k;l-cose)b|  will  be  >  1  for  typical  sample  i 
unless  than  angle  is  too  small  to  be  experimentally  reaiizucle. 

Thus  taking  |a/6!  <  1  will  be  a  good  assumption  for  almost 
meaningful  experimental  configuration. 

Expanding  the  exponent  inside  the  integral  from  Equation 

seep  m  a  term,s  first  order  in  a,  I  ,  '17-46)  becomes 

s  i 


one  re 


ire  ns : 


any 
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U) 


„ idaz 
a l  I  aC- 

„  ctz  1  .  .  :aa 

=  C  H  -  - ; -  >  3z\  1  -  —  sin 

~s  Jt5 


-  —  (4  +  -T-)  COS  (  3z-ri)  }  . 


>er forming  the  integration,  this  becomes 


,ic,iz  - 

JL  X.  k  31 


1  3aa  - 

—  i  z  f  — —  cos  (rz+i/  - 


3“  I 
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(ne  car.  use  the  boundary  condition  A  ,T~)  =  0  to  cfc 


id 

^  ct  , 


-jaa 


+  — r —  CCS ( £L+u) 


-j  SI 


I 


te  can  then  substitute  <.17-49)  ir.to  '.17-43)  to  ob 


aaC, L~  i' 


\  <L-z)  -  ~ 


oaa  - 

_s  ir.  i. — — 


.a  .  :i,  r  .3'u-z) 

— —  ',4  -  — r)  Lsin' — t - cos  — - 


Cne  :an  expand  .17-42)  to  first  order  m  a, I  to 


- 1  cos  ,  . 


tne  car.  cow  substitute  (17-51)  into  (17-50)  to  obt 
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odd  r  ,3vl< 

— —  lsm( — 7- 


+  2.)  sin  (- 


^Lj  1 


2a  ,  iuN  r  ,d(L-z)  x  ,^Uz) 
— —  +  — )Lsin( — r )  cos  ( — 

Z>  X  2  Z  Z 


At  the  output  face  z  =  o. 


-iaA. !L) 


3s in !  ib-i;  -  .cm  tL-1; 


fL.  -i  2a  ia,  .  fL, 

ca>— )J  '  —  '4  '  —'sir.  ,—  ccs  r 


The  output  intensity  of  the  3ragg  scattered  signal  bear,  will  be 
*■» 

equal  to  A ^ 1 0 )  | “ ,  and  to  first  order  in  a7 I  ,  this  is  given  by 


2.r2:4 


-2a  L 

1  R 

•>  ia1(l)  ri 


'r  2aI-  r  '  .  . 

'L  *  rt—  Lid.  Sin  1  zL~u  -  l  cos 'in- 

i  x  ri  1 


— —  sin\—  +  2, )  3 m  v  — ; 


■4a  ^  f»  4  t  2^ 

__  ,  4  _  -j-)  l_Sm  (■— ^)  cos  (-p  -  2.)]} 


I  +  I 


=  -  vx 


50 
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Equation  (17-54)  is  che  expression  for  the  output  2ragg  scattered 


jigr.al  team  to  first  order  in 


:A2A4 


- .~\b L.  _  *  iS na i  i  ir.C»SSj  . 


Examining  this  result  more  closely,  one  car.  see  from  Equation 

rt  .*1 

17-54;  that  ail  of  the  terms  in  (17-54)  that  are  first  order  m  — g- 

C 

ill  vary  sinusoidally  with  the  phase  mismatch  1  of  the  two  pump  beams 
If  the  2  oumc  beams  are  within  the  coherence  let  nth  cf  each  ctr.e 


tr.en  tr.e  or.ase  mismatcn  „  will  oe  riven  ov 


■v r.ti  - w  - 


the  difference  in  oath  lencth  between  the  twc  turns  beams. 


Thus  the  ccm.ccnenc  cf  A. (C)  whicn  is  first  order  m 


exceptionally  sensitive  tc  vibration  and  to  tr.e  precise  alignment  of  the 
comp  beams.  Thus  at  verv  small  ancles,  Equation  17-54)  srodiscs  the 


observed  3rana  scacts 


;::r.a.  seam  to  have  a  comconent  tr.at  osci_- 


-c-tes  very  rapiu-y  cue  to  visr  atior.s  or  miner  ac;  istner.es  superimposed 
on  an  "envelope  signal"  that  is  insensitive  to  vibrations  and  precise 
alignments. 

Remembering  that  tr.e  "scattering  efficiency"  g  cf  the  four  wave 
mixing  process  is  the  ratio  of  the  exiting  signal  beam  A.  )0V  *"  to  the 
exiting  probe  beam  m  the  acsence  of  the  pump  interactions  A,  Z, 

„  ,  one  obtains  from  (IV-43j  that 


A,  '0) :  =n  =  :a. ,b) 


17-34),  the  scattering  efficiency  -  ftr  the  3ir.al2 


:r:er  m 


’’(3 


2a  L 


r2:4  1+62  "2ClRL 

- — )  4  l 


R  .  ,3L  , .  .  ,3L, 

— —  sm  (—  +  A)  sin  (— ) 

SI 

s 


2iiL  r  ^ 

{L  -  —  [  2  sin  l  sLf  A)  -  u . cos ( ub+_) J 

■2 1  rf  l 

3 
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Transient  Sehavior 


When  the  pump  beams  are  chopped,  two  processes  tan  contribute  to  trie 
decay  of  the  sinusoidal  excited  state  population  prating.  The  grating 
can  decay  because  of  the  decrease  in  the  excited  state  population  cy 
normal  fluorescence  decay  and  because  of  exciton  migration  from  the  peak 
to  the  valley  regions  of  the  grating,  whenever  the  exciton  motion  is 
diffusive,  the  density  of  excited  states  is  given  by  (26) 


3  n(x,t) 

r 


n 


.  4L 

o  n'x,t) 
:< 


where  x  is  along  the  direction  of  the  grating  wavevectc-r,  2  is  the  dif¬ 
fusion  coefficient,  and  t  is  the  fluorescence  decay  tine.  If  it  is 
assumed  that  n  { x ,  t )  initially  has  a  sinusoidal  spatiai  distribution,  the 
solution  to  (IV- 56) 


n  v  x ,  t ) 


-  k  ^  0 1 

i  -  cos : kcx) } 


where  k  is  the  magnitude  of  the  crating  wavevector.  The  depth  of  the 
grating  in  is  given  by 

K 

In  =  n ix=0, t!  -  n(x  =  f.  t)  =  i  ~  [ZV-5S, 
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where  K  is  the  decay  constant  given  by 


2  1 
K  =  2  :  k  D  *  — ) 


rhe  magnitude  of  the  grating  wavevector  k  is  given  by 


T/sin(e/2) 


for  small  3,  !<  may  be  approximated  as 


Since  the  scattering  efficiency  is  proportional  to  In  (43),  the 
race  scattered  signal  beam  intensity  I_(t)  may  be  'written  as 


I  (t)  *  I  Ll  In/ 
s  o  w 


I  (t;  *  :  :2i"'Nt 
s  o  w 


:;v-dO) 


where  I  is  the  incident  orobe  beam  intensitv  and  1  is  the  cumc  beam 
o  "  w 

intensity.  Thus  the  3ragg  diffracted  signal  beam,  should  decay  exponen- 
tiallv  with  iecav  constant 


sin" 


which  for  small  angles  becomes 


_  2  _  3? 
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III.  PROPERTIES  OF  Nd  Y,  P,0,  .  AND  Nd  La,  PrO, .  CRYSTALS 

x  1'X  )  14  x  1-x  5  14 

A  significant  amount  of  time  was  spent  characterizing  the  optical  proper¬ 
ties  of  stochiometric  laser  materials  for  possible  use  in  minilaser  applica¬ 
tions.  The  most  significant  result  of  these  investigations  are  summarized 
in  the  following  three  manuscripts.  Additional  information  on  the  quantum 
efficiency  of  these  materials  is  presented  in  Section  V.I. 
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III.l 


LIFETIME  MEASUREMENTS,  INFRARED  AND 
PHOTOACOUSTIC  SPECTROSCOPY  OF  NdP5014 
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\bstract  We  have  malic  a  series  of  spectroscopic  measurements  to  heller  unilerstano  Ihe  optical  pro|serties 
and  concentration  quenching  characteristics  of  the  stoichiometric  laser  maieri.u  NdP.O,..  7  he  icsuils  ■  f 
fluorescence  lifetime  measurements  in  the  presence  of  JiHerent  surface  environments  mdicaies  that  the 
surface  condition  alfeeis  the  concentration  quenching  whereas  lifetime  measurements  under  applied  uniaxial 
stress  show  that  the  presence  of  internal  strains  in  the  crystal  is  not  significantly  effective  in  iluorescence 
quenching.  A  comparison  of  ihe  internal  rellection  spectrum  with  the  normal  infrared  spectrum  did  -tot 
reveal  any  significant  differences  in  ihe  Nil1'  energy  levels.  Photoacoustic  spectroscopy  results  proved 
difficult  to  interpret  but  appear  to  he  consistent  with  the  presence  of  surface  quenching  ol  cveitons  and 
indicate  that  the  intrinsic  quantum  efficiency  of  NdJ  *  ions  in  ihe  penluphosph.ue  host  in  the  absence  of 
concentration  quenching  is  approximately  0.9(1. 


I.  INTRODUCTION 

Understanding  [lie optical  properties  of  NdP,Ol4  is  of 
ugnilieani  importance  because  ol  ihe  use  of  this 
material  in  mmtlaser  applications  [1  \  One  of  the  most 
interesting  properties  of  this  material  is  the  small 
amount  of  concentration  quenching  of  the  fluo- 
reseence  compared  to  that  iff'  NJ3  *  in  other  types  of 
hosts  There  appear  to  he  two  possibilities  Tor 
explaining  concentration  quenching  in  NdP,Ol4.  The 
tirst  is  energy  migiation  In  "sinks"  which  are  delect 
sites  al  which  the  energy  is  dissipated  radiation- 
iesslv  2  7  1  he  second  is  ihe  quenching  of  exeitons  al 
the  surface  of  the  crystal  '37  We  describe  here  the 
results  of  several  different  types  of  spectroscopic 
measurements  which  provide  lurther  information 
concerning  the  characteristics  ol  fluorescence  quen¬ 
ching  IP  Sdl’.U.x. 

We  recently  published  the  results  of  an  investigation 
of  Sd.V,  ,Pt(J.l4  crystals  using  time-resolved  site- 
seieciion  spectroscopy  techniques  [3. 4].  One  of  the 
conclusions  drawn  from  this  work  was  that  spatial 
energy  migration  was  taking  place  in  ihis  material 
without  spectral  diffusion.  More  recent  fluorescence 
line  narrowing  experiments  corroborate  this  conclu¬ 
sion  at  least  at  low  temperatures  [57.  The  observation 
,ff  no  energy  transfer  between  ions  in  nonequivalent 
.  rvstal  field  sites  casts  some  doubt  on  the  possibility  of 
a  quenching  transfer  occurring  fo  ions  in  randomly 
distributed  sink  sites  However,  a  consistent  interpret¬ 
ation  of  the  results  can  be  provided  by  a  model  of 
trapping  and  quenching  of  exeitons  at  a  macroscopic 
delect  region  such  as  the  sample  surface. 

We  report  here  the  results  of  measuring  the 
fluorescence  lifetime  of  Ndl\Ou  crystals  in  the 
presence  of  different  surface  environments  and  ihe 


results  of  laser  pholoacoustic  measurements.  Al¬ 
lhough  ihe  results  do  not  conclusively  prove  which 
model  of  concentration  quenching  is  correct  for  this 
system,  they  are  consistent  with  an  interpretation 
based  on  surface  quenching  and  cannot  be  easily 
explained  by  migration  to  randomly  distributed  sinks. 
In  addition  measurements  ol  the  lluoteseenee  lifetime 
under  conditions  of  uniaxial  stress  indicate  that 
internal  strains  are  not  effective  m  producing  increased 
quenching  as  would  be  expected  !pr  sink  ales  IR  ..nil 
internal  rcffeclion  spectra  were  obtained  m  an  attempt 
to  characterize  quenching  centers  hut  tl  e  results  were 
inconclusive. 

:.  LIFETIME  MUM  HE  MEATS 
The  Iluorescence  lifetime  of  NdJ*  emission  can  be 
used  as  a  measure  of  the  rale  of  concentration 
quenching  with  the  quenching  rate  defined  as 


where  is  the  intrinsic  decay  time  and  :  is  Ihe  decay 
time  measured  in  the  presence  ol  quenching.  In  lightly 
doped  crystals  of  Nd,Y.  J\Ou  Jl  room  icmpcra- 
lurc  the  value  of  is  measured  to  be  about  32(1  psec 
whereas  for  concentrated  \dP,0M  (he  measured 
values  of  :  range  from  about  100-1 20  usee  at  room 
temperature  depending  on  ihe  exact  wavelength  ol 
excitation  [3  7  We  performed  a  series  of  measurements 
of  the  fluorescence  lifetime  of  NdP.  0, 4  crystals  with 
different  types  of  perturbations  affecting  either  the 
surface  or  interior  environments  of  the  NdJ*  ions  to 
determine  how  such  environmental  changes  affect  the 
concert  (radon  q  ucnching 
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(a)  Ejjec:  uj  surface  environments 

'he  fir  s'  of  these  experiments  consisted  ol 
measuring  llie  ihioresccnce  lifetime  of  the  crystal 
immersed  in  sarious  types  ol  solutions  which  should 
perturb  the  environment  of  NdJ  ’  ions  at  the  surface  ol 
the  crystal.  Measurements  were  made  at  room 
temperature  using  the  same  experimental  equipment 
and  procedure  described  previously  .3’  with  the 
excitation  laser  tuned  to  give  a  value  of  r  x.  12U nsec 
for  an  air  or  vacuum  environment.  We  found  that  in  a 
concentrated  solution  of  hydrofluoric  acid  the  decay- 
time  was  decreased  to  about  104  nsec,  (n  concentrated 
solutions  of  carbon  tetrachloride  or  carbon  disulfide 
the  decrease  in  lifetime  was  about  half  this  amount 
w  hereas  when  the  sample  was  immersed  in  acetone,  no 
change  in  lifetime  was  detectable. 

To  further  check  the  effects  of  the  surface 
concentration  quenching  a  sample  was  grown  which 
had  an  \dP<014  substrate  and  an  epitaxial  layer  of 
Gd„  .4La„  40PxO,4.  !ror  the  preparation  of  the 
epitaxial  layer,  the  methods  used  for  the  growth  of 
crystals  were  employed.  The  seed  was  lowered  into  the 
solution  which  was  kept  at  a  temperature  in  the  range 
of  500  -  5  50J  C  Epitaxial  growth  proceeded  for  three 
days  until  a  1  mm  thick  layer  was  grown.  This 
corresponds  to  a  relatively  fast  growth  rate  of 
14  urn  hr  A  smail  region  of  the  substrate  at  one  end  of 
the  sample  was  left  uncovered.  To  qualitatively 
evaluate  the  perfection  of  the  overgrowth.  Lauc  X-ray 
buck  reflection  patterns  were  obtained.  Sharp  spot 
patterns  with  the  same  symmetry  confirmations  were 
obtained  from  both  the  substrate  and  overgrowth 
indicating  (hat  both  regions  were  single  cry  stalhne  and 
i  here  is  no  misalignment  between  them.  No  evidence  of 
mosaic  structure  or  lattice  constant  tnismuicn  could  be 
found. 

A  comparison  was  made  of  the  fluorescence  lifetimes 
measured  when  the  sample  was  excited  through  the 
expitaxial  layer  and  in  the  region  of  uncovered 
substrate  material.  At  room  temperature  with  the 
particular  laser  excitation  wavelength  which  was  used 
the  lifetime  measured  with  substrate  excitation  was 
1 10  nsec  whereas  the  lifetime  measured  after  excitation 
through  the  epitaxial  layer  was  11‘Jirsec. 

The  fact  that  the  fluorescence  lifetime  can  be  altered 
bv  the  surface  environment  of  the  crystal  indicates  that 
the  surface  plays  some  role  in  the  concentration 
quenching  of  fluorescence  in  NdP50,j.  However,  it 
has  not  yet  been  possible  to  find  the  environment 
necessary  for  significantly  reducing  the  quenching. 

lb i  t  mu i\iai  stress  measurements 

Fluorescence  lifetime  measurements  were  made 
while  samples  were  under  uniaxial  xtiess  to  determine 
the  contribution  to  the  fluorescence  quenching  due  lo 
energy  migration  to  randomly  distributed  sinks.  A 
sample  of  Nd„  4Y„  ,P,0,4  was  mounted  m  a  dewar 
containing  an  hydraulically  driven  piston  and 
immersed  in  liquid  nitrogen  The  sample  was 
compressed  along  the  x-axis  and  excited  by  the 
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Fig.  I  Fluorescence  .aetanc  ol  lc  A  eve:  .!  "6  -  .. 

lunctionoi  umjxi.n  tre  o  .pph-.-u  m-iic  o:.-  .  n; 

nitrogen  laser-pumped  lye  ,a.e:  at  '•■/  \  ! -v: 

fluorescence  liietime  was  monitored  it  v  ..rn  ip  e.s  ; 
applied  sticvs  The*  results  ate  how:,  .;  I  .g  I 
application  I  exiremeiv  high  levels  ol  si;.--  .  up  ' 
’OOOpsi.  results  in  only  as;  addition  a!  in"  ,  quetidia.g 
of  the  fluorescence  lifetime. 

Sinks  lor  radiationless  quencr.ing  o;  energy  ire 
genet  ally  associateq  with  V;-1'  urns  ioc.iteo  near 
crystal  imperfections  where  iocai  -trams  .e 
differences  ,n  the  energy  ievel  structure  allowing  .;o--- 
relaxation  to  occur  rhus.  the  implication  e\te:  -s.u 
stress  would  be  cxj'oeted  i<  introduce  addition.:, 
strains  in  the  sample  and.  therefore,  ni.res.c  the 
number  of  sink  -ites  resulting  m  ennanced  uueii.  hu.c 
It  is  well  known  that  nendvmium  pentaphosp,.  ,-e  i. 
lerroelaslic  crystal  .mo  twin  boumi.oic  .a,,  v 
introduced  into  the  crystal  ">>  very  sitia.i  iorce  1 
The  fact  that  only  a  sm.nl  .iiange  in  quench  ng  -s 
observed  even  w  uh  a  large  magnitude  oi  ippiiexi  .tre 
indicates  that  randomly  distributed  sink .  are  probanly 
not  the  sites  responsible  for  radiationless  quenching  m 
this  case. 

.1.  NPhCTRONt  OIM<  \n  VNI  RI-MtNfs 
la l  IK  spectra 

If  radiationless  quenching  occurs  at  the  -ur:,..e  of 
NdP.O,4  crcstals  nut  not  in  the  hulk,  .t  might  be 
expected  that  the  energy  leveisofNd3'  ions  ut  ..I'.ave 
sites  are  shifted  in  position  with  respect  to  those  ol  ions 
m  normal  ales  m  such  a  way  that  cross  relaxation 
transitions  are  more  favorable  One  way  toi  tin.  to 
occur  is  the  shifting  of  ihe  ■*! . «  ,  level  to  lower  energy 
An  attempt  was  made  to  detect  such  a  shill  S\ 
comparing  the  internal  reflection  spectrum  to  the 
ordinary  transmis.ix'n  ivctrum 

The  internal  reflection  spectrum  of  NdP  O  .  wa. 
obtained  m  the  region  of  the  ‘I.,  ..  'I.  ,  .  md  ‘  I , .  . 
levels  using  a  I’erkin  l  Inter  modci  -I  .p..lropiioto 
meter  with  the  Wilks  reflectance  attachment  pie 
cleaved  samples  were  mounted  on  iha.i  um  ho  mo 
iodide  KRS-5|  prisms  with  index  -nai.lune  'aim 
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Fie  1  IR  la  1  ,.nd  internal  reflection  (hi  spectra  of  NdP,0M 
crystals  at  room  temperature 

in  =  I.mOi.  The  b  direction  was  chosen  because 
NdP,0|4  can  be  cleat ed  perpendicular  to  the 
n  axis.  The  flatness  01  the  (010)  cleaved  surface  is 
t>  picallv  HI  11m  over  the  total  area  of  the  sample.  This 
eliminates  polishing  procedures  and  yields  a  surface 
free  of  contaminants  The  Wilks  micro-sampling  trans¬ 
mission  attachment  was  used  to  obtain  a  transmission 
spectrum  under  similar  conditions  to  those  used  in  the 
internal  rejection  measurement.  The  two  spectra  ate 
shown  in  big.  2  Wide  slits  |4o0,uni|  were  used  to 
improve  the  signal  10  noise  ratio.  Although  this 
broadens  the  spectra,  no  shifts  in  the  centroids  of  the  JI 
mulnpieis  can  lie  detected. 

The  penetration  depth  of  the  evanescent  wave  in  the 
rarer  medium  in  internal  reflection  is  tv  picallv  on  the 
order  of  0  2/  In  the  spectral  region  investigated  the 
penetration  depth  would  be  approx.  U  3  n  V  uni  which 
is  .lboul  t(H)  500  times  the  latnce  spacing  in  the  h 
direction.  The  lack  of  observ  ation  of  a  spectral  shift  in 
the  internal  reflection  results  does  not  necessarily 
prove  that  the  energy  levels  of  Nd3'  ions  in  surface 
sites  are  not  perturbed.  It  may  be  that  the  region 
sampled 's  so  large  that  the  mammy  of  ions  involved  in 
transitions  are  located  m  bulk  sites  instead  ol  surface 
-ties 

Another  point  to  note  about  the  infrared  spectrum 
m  Fig.  2  is  that  no  impurity  bands  can  be  observed  in 
the  region  between  2S00  and  3800cm  '  li  is  well 
known  that  phosphate  materials  can  contain  hydro¬ 
gen  impurities  which  can  act  as  quenching  centers  lor 
luminescence.  Tins  does  not  appear  to  be  the  case  here 
since  OH  vibrationai  transitions  would  be  easily 
visible  in  the  near  IR  region  of  the  spectrum 
investigated. 

fb)  Photoacousttc  spectroscopy 

Photoacoustic  (PAl  spectroscopy  is  a  technique 
whidi  directly  monitors  radiationless  relaxation 


processes  by  detecting  the  .ouiul  produced  in  an 
acoustical  cell  when  j  sample  is  excited  isv  a 
periodically  varying  light  source^7  We  previously 
reported  on  the  PA  spectra  of  Nd,Y,  .  ,P5Om  crystals 
and  the  basic  experimental  apparatus  is  described  in 
Rel.  [3],  Two  additional  types  of  PA  measurements 
have  now  been  performed  on  this  set  of  samples  using 
the  individual  tines  of  an  argon  ion  laser  for  excitation. 
This  light  was  chopped  at  frequencies  varying  from 
HXI-2700  Hz  and  the  laser  power  was  continuously 
monitored  and  stabilized  at  a  level  of 0.15  W  The  PA 
signal  intensities  from  the  samples  were  more  than  two 
orders  ol  magnitude  greater  than  noise  and  einptv  cell 
signals. 

The  lirst  experimeni  was  to  record  the  maximum  PA 
signal  from  the  various  samples  at  twelve  or  more 
chopping  frequencies  v  alter  excitation  with  the 
5 145  A  argon  laser  line.  The  results  are  shown  in  Fig.  3 
for  the  NdP,014  and  Nd,,  'tlJ.,P.O.j  samples.  The 
lightly  doped  sample  shows  a  PA  signal  intensitv 
which  varies  as  y  *  throughout  the  whole  range  of 
frequencies  investigated  whereas  the  concentrated 
sample  signal  varies  as  v.  1  at  frequencies  up  to  about 
450  Hz  and  then  varies  as  »v  3  ;  at  higher  frequencies. 
Samples  with  intermediate  concentrations  exhibit  PA 
signal  intensities  which  vary  with  chopping  frequency 
between  a  y‘ 1  and  v~ 3 

The  thickness  of  these  samples  is  of  the  order  of  one 
millimeter  and  the  optical  penetration  depth  (I 
=  x~  ' )  varies  between  about  700/mi  lor  the 
concentration  sample  and  '(XX) /tm  for  the  I0"„ 
sample  Jl  this  excitation  wavelength.  According  to  the 
theory  of  Rosencwaig  and  Gerslto  (RGi  the  thermal 
diffusion  length  varies  wall  chopping  frequency  as  8] 


3  PA  Mcn.il  'nlcnsil)  \s  choppmc  Ircqucncv  at  room 
temperature  lor  A  c\i:uuon 
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where  ;i  is  the  thermal  dtffusivity.  For  NdP,0’.i  the 
largest  value  of/J  is  [9]  ~  10“ '  em:  >ec  1  and.  thus.  f.h 
ranges  from  about  54  gm  at  v.  =  1 10  Hz  1  I  jun  at  r 
=  3700  Hz.  It  is  approximately  37  ;uil  at  450 11/  where 
the  change  in  dope  occurs  in  the  curve  lor  the  100",, 
-.ample  as  seen  in  Fig.  3. 

From  the  above  considerations  the  results  of  these 
ex periments  should  be  classified  in  terms  ol  the  special 
cases  of  the  RG  theory  as  an  optically  thin-lhermally 
thick  case  for  the  10",,  sample  and  an  optically  thick  - 
thermally  thick  case  for  the  I00‘‘„  sample.  Since 
/,h  s  f in  all  cases  the  PA  signal  intensity  is  expected 
to  vary  as  v*  ;  The  fact  that  this  is  not  observed  is 
probably  due  to  the  fact  that  the  RG  theory  is  based  on 
a  one-dimensional  model  where  the  thermal  gradient 
at  the  excited  surface  causes  uniform  heat  tlow  back 
into  the  cell.  In  these  experiments  a  small  area  ol  the 
total  sample  surface  is  excited  by  the  laser  beam  and 
the  heat  is  generated  in  a  cylindrical  volume  within  the 
sample.  The  heat  generated  in  a  layer  close  to  the  front 
surface  still  is  dominated  by  the  temperature  gradient 
at  the  surface  and  contributes  to  the  PA  signal  as 
predicted  by  the  one-dimensional  theory.  However, 
further  into  the  sample,  the  heal  Row  is  dominated  by 
the  temperature  gradient  radially  out  from  the  excited 
cylinder  and  most  -if  it  will  never  reach  the  surface 
within  a  duty  cycle  to  contribute  to  the  detected  P  \ 
signal.  Comparing  these  considerations  w  ith  the  results 
obtained  on  the  It)'1,,  sample  indicate  that  the  effective 
optical  layer  depth  lor  heat  contributing  to  the  PA 
signal  /  tl  is  CIO  urn.  the  smallest  value  of  f,h.  For 
....  lh  the  RG  theory  predicts  a  i\.  1  variation  of  the 
PA  signal  intensity  as  observed.  For  the  ItX)"  ,  sample 
applying  the  above  considerations  to  the  data  shown 
in  Fig.  5  indicate  that  the  effective  optical  layer  is  two 
and  a  half  times  greater  than  that  in  the  !0"„  sample  A 
possible  explanation  for  (his  is  the  added  contribution 
to  the  PA  signal  from  the  hcjl  generated  by  the 
diffusion  and  surface  quenching  of  excitons.  Although 
these  PA  results  are  obviously  too  complicated  to  be 
considered  proof  of  exciton  migration  and  surface 
quenching,  this  interpretation  is  at  least  consistent 
with  previous  results  indicating  exciton  migration 
lengths  of  the  order  of  tens  of  microns  and  quenching 
at  the  surface. 

The  second  PA  experiment  performed  on  these 
wimples  was  an  attempt  to  measure  their  absolute 
quantum  etliciency  The  maxima  PA  intensity  and  the 
phase  at  signal  maximum  was  measured  for  two 
ditferenl  wavelengths  of  excitation  at  two  different 
chopping  frequencies.  The  4765  and  5I45A  lines  of 
the  argon  laser  were  used  for  excitation  and 
measurements  were  made  with  =  312  and  KXXIHz. 

The  Nd '  *  energy  level  model  used  to  interpret  these 
data  is  shown  in  Fig.  4  After  absorption  into  levels  F- 
or  F-  heal  is  generated  by  radiationless  relaxation  in 
the  metastable  level  F,.  This  level  has  radiative 
transitions  to  the  various  multiplets  of  the  il,  term 
with  branching  ratios  There  arc  also  ihe 
possibilities  of  nonradialive  decay  and  concentration 


I  ig  4  Inert!)  levels  .mil  transition i  :<>r  Nu*'  i*r.> 


quenching  in  which  all  of  the  energy  in  is  lost  to 
heat.  After  transitions  to  the  higher  !>  ing  ground  term 
levels,  heat  is  generated  hy  radiationless  relaxation  ti¬ 
the  iowest  ground  level.  These  levels  .md  transition 
rales  arc*  labeled  in  Fig.  4.  The  quantum  ctlic.ency. 
probability  ofmniradiative  relaxation,  and  probability 
of  concentration  quenching  of  the  tiielastaole  state  are 
given  by 

Qt  —  U ...  1 11  /<  ~  H-,#  *  H.  I 

^\K  ~  f  ~  1 1  l  I 

P,  =  IV,.  (IF,  -  U'%)(  -  H\  i 
with  tlie  condition 

y/i  ~  - /\  -  I  14. 

file  expression  for  the  PA  signal  at  ptu.se  angle  •• 
alter  excitation  in  energy  level  E_,  is 

IjU I  =  ClP^  -  tan  !G-v  r  i  -  " 

1 5 1 

where  l\  is  the  power  absorbed.  (  :s  a  factor 
accounting  for  properties  of  the  -c*!!  and  detection 
system,  and  the  summation  is  over  ail  transitions  in 
which  heal  is  produced,  o”  is  the  probabiiilv  tor  the 
nonradialive  transition  to  occur  between  levels  i  and  -. 

is  the  energy  of  this  transition  tan  1 1 2 rt >  _  t ,  i  is  the 
phase  shift  of  the  signal  due  to  the  lifetime  of  the  initial 
level  and  yr  is  the  phase  shift  due  to  the  detection 
procedure.  For  \dP.O,4  til  of  ;he  energy  level 
positions  and  branching  ratios  are  known  Id  t  >:il\ 
the  mctastable  levei  has  j  lung  enough  lifetime  to 
introduce  a  measurable  phase  shut  into  the  Mgn.il. 
Evaluating  the  sum  in  eqn  i5i  exphctiK  leads  to 

/  ||l|  =  C{P  E  j  /  .  Co v I y/  -  n. 

■*  1,1  K>| 

~  (FM  —  (.(JF icosiy.  •  i.m  i3.ni-:,i  -d 
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Mhcrc  is  the  lifetime  of  ihe  meustable  level  and 

f  =  -  h}iEu  —  />«,£,,  (7) 

In  order  lo  eliminate  the  unknown  experimental  factor 
we  took  the  ratio  of  PA  signals  at  two  different 
excitation  wavelengths.  Solving  for  QE  gives 

(JE  = 

AE-,  cos  i to  -  i/-l  -  .•<£,,  cos  tan  't2itvt:,i  -  il.  i 
k  -I  cos  li//  -  tan  1  (2.tv  r,|  —  n.) 


a  acre  the  h  and  '  subscripts  refer  to  excitation  at  5145 
aim  4~t>5  A.  respectively  and  the  factor  1  is 

■i  - /. 1 1 /•;„ ,  E-,  nf.  i\i.  i»i 

SnisC  me  values  to  be  used  for  m  are  still  in  question,  a 
sowoiki  equation  was  generated  by  maximizing  cqn  ioi 
with  rcsjx-st  to  •’  nut  solving  for  QE  which  gives 


liquation  v S >  and  the  twoequations  obtained  front  eqn 
■  o  lor  two  ihilereni  excitation  wavelengths  can  now 
l,e  so.vea  imult.meousiy  tn  an  iterative  wav  to  obtain 
a  .miotic  value  lor  QE.  To  check  the  uniqueness  of  the 
result,  the  procedure  w-us  repeated  lor  both  chopping 
avquenvies 

I  or  the  Sdl’.O.j  sample  a  quantum  elliciencv  ol 
i)  45  was  found  b;.  the  method  described  above  In  this 
,ase  the  power  adsorbed  in  eqn  i9i  was  assumed  to  be 
me  total  incident  power  for  both  wavelengths  of 
excitation.  The  value  of  the  quantum  elliciencv  ;n  the 
ihscnce  of  concentration  quenching  can  be  found  from 
the  expression  OEtO 1  =  ijEi  c:.  The  combined 
lifetime  and  Qt tvi  results  indicate  that  QElOl  rs 
approximately  n  )i|  for  NdP.O)4.  Similar  analysis  of 
the  v  =  o20  sample  gives  ijt't.vi  =  ').b2  ind  QElO  i 
=  Odd 

The  above  procedure  for  determining  quantum 
elliciencv  b>  taking  the  ratio  of  two  PA  signals  was  lirst 
suggested  by  Rocklev  I!  and  applied  to  a  simpler 
case  where  phase  angles  for  different  transitions  did 
not  piesenl  a  problem.  Although  the  results  obtained 
on  SdP.t)14  appear  to  he  reasonable,  it  should  be 
pointed  out  that  the  value  obtained  for  QE  is  extremely 
■ensitive  to  the  parameter  I  and  even  very  stnail 
changes  ;n  factors  .uch  as  the  ratio  of  power  absorbed 
can  cause  sigmhcaiit  changes  in  die  predicted  QE.  This 
is  especially  a  problem  when  we  are  dealing  with  a 
situation  having  involvement  with  an  effective  thermal 
lavet  as  discussed  previously  and.  thus,  the  value 
obtained  for  Q E 1 0 1  mud  be  considered  as  only  a  very 


rough  estimate.  Quimhy  and  Veil  12  have  reported 
measurements  of  the  quantum  efficiency  m  Nq  '  ’  .ons 
in  a  glass  host  using  PA  measurements  I  he  model 
used  in  the  interpretation  of  their  data  is  similar  to  the 
one  used  here  except  that  phase  shifts  for  dnleient 
transitions  were  not  considered.  Their  results  indicate 
a  Qfc(O)  of  between  o.()5  and  0.75  lor  \d3  ’  in  glass. 


-  ,  cos i ip  -  0b)  -  cosly  t.m  'i2t:>  t,i  -  *'„/ 

-cosly  -  tan  1  i2.t.  :.i 

n  i 


4.  MV1MVHV  tSDlOMl.tslOSS 

The  various  results  described  here  provide  further 
mloimation  on  : lie  optical  properties  of  NdP.O  .  „ 
Holh  ihe  lifetime  and  P\  measurements  indicate  thai 
the  surface  pla's  some  role  :n  the  concentration 
quenching  of  this  malcnal  bui  ihe  "Csulls  arc  m  l 
.  ''idavive  n  identifying  '.he  quenching  mechanism 
I’hoio.icousl.c  measurements  made  r>>  \u/e;  n 

powdered  samples  indicate  that  radiatiomess  quench- 
mi!  is  create:  <  m  powders  of  jo  uni  grain  size  than  lor 
'urge  single  srvsi.ds  which  is  ais..  cons. stent  with  'he 
piedictioiis  ol  suit.ivc  quenching  I  •  More  direct 
measurements  ol  evenon  diffusion  U'lne  iiaitsienl 
gratmg  techminies  .ire  presently  heme  used  in 
leiermme  the  possible  contribution  of  tone  u.ce 
even  >n  imgrudon  m  :i:e  quenenme 

•U  (inmlci/ip  rm  ula  III.,  work  was  supported  bv  i.hc  L  S 
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\b.s  tract.  \  afiations  m  the  MuoreseetKc  line  positions,  hnvrwidlhs.  and  lifetimes  *iih  hydio 
static  pressure  up  •«>  aroul  sO  kb  ir  ire  reported  tor  samples  ut  NdP  O  .  and  s'il  V ,  ,P.t )  , 
Ihe  icsults  for  the  hentlv  doped  sample  are  consistent  with  changes  in  the  radiative  dcv.iv 
rate  while  the  vor\  different  results  obtained  on '.he  heavily  doped  sample  are  inter,  ted  in 
teimsol  increased  radiationless  uiienehmu  due  to  energy  migration  to  trapping  sites  Results 
are  tlso  reported  on  ’he  variations  m  the  Muorescence  lifetimes  oi  the  R  lines,  ano  the 
positions  •  *f  the  N  lines  and  v ihromc  peaJw  of  ruby  as  a  Junction  «>!  In  arosutic  pressure 


I.  Introduction 

The  study  of  the  effects  of  pressure  on  the  optical  properties  of  impurity  ions  in  solids 
can  be  useful  in  many  regards.  Changes  in  ihe  widths,  positions,  and  lifetimes  of 
electronic  transitions  and  their  vibrational  sidebands  as  a  function  of  applied  pressure 
can  give  important  information  concerning  the  interaction  of  the  impurity  iou  with  its 
surroundings  (Drotning  and  Drickamer  ISC6.  Tyner  and  Drickamer  1^77.  Klick  et  ul 
llC7,  Webster  and  Drickamer  1980).  Diamond  anv  il  high-pressure  cells  are  convenient 
tor  such  investigations,  requiring  the  study  of  very  small  samples  In  this  regard,  the 
light  intensity  and  small  size  of  the  illuminated  area  available  with  lasers  are  valuable. 
The  work  reported  here  represents  an  exploratory  application  of  the  techniques  of 
pulsed  laser  excitation  and  gated  electronics  for  time-resolved  fluorescence  studies  at 
high  pressures  in  a  diamond  anvil  cell 

The  measurements  reported  here  include  studies  of  the  effects  of  pressure  on  the 
absorption  and  fluorescence  spectra  and  the  lifetimes  of  Nd,  Yi  -  ,ILOu  for  x  =  (1. 1  and 
1.0.  in  recent  years  there  has  been  significant  interest  in  understanding  the  optical 
characteristics  of  NdPOu  because  of  its  properties  as  a  stoichiometric  laser  material  for 
minilaser  applications  (Weber  1975).  One  of  the  most  important  properties  of  this 
material,  vv  hich  is  still  not  understood,  is  the  nature  of  the  mechanism  giving  rise  to  the 
weak  concentration  quenching  of  the  fluoresence.  These  pressure-dependent  results 
give  new  information  on  the  fluorescence  quenching  in  this  material  which  is  consistent 
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with  an  energy  migration  model.  In  addition  results  are  presented  on  the  pressure 
dependence  of  the  tluorescence  spectra  and  lifetimes  of  heavily  doped  ruby  crystals, 
see  section  V.2. 

2.  Experimental  details 

The  diamond  anvil  high-pressure  cells  used  for  this  study  were  designed  by  one  of  us 
(ILS)  and  are  very  similar  to  a  cell  described  previously  (Yu  er  al  1979)  The  major 
modification  made  in  the  cells  used  here  is  the  opening  of  a  conical  hole  of  16°  included 
angle  in  the  piston  and  driving  screw,  allowing  optical  access  with  approximately  f  3.5 
aperture.  A  fluid  of  approximately  four  parts  methanol  to  one  part  ethanol  served  to 
maintain  hydrostatic  and  homogeneous  pressures. 

Pressure  was  monitored  by  observation  o!  the  wavelength  shift  ot  the  R,  fluorescence 
line  of  ruby .  For  pressures  in  the  range  of  interest  to  this  work,  this  shift  has  been  found 
to  be  +  0.3b5  A  kbar  1  (Piermarini  et  al  1975).  Also,  a  correction  was  made  tor  the 
temperature  shift  of  the  line,  which  is  -1-0.007  A  K  1  near  room  temperature  (Paetzold 
1951 )  due  to  small  drifts  in  room  temperature,  at  which  all  measurements  were  made, 
t'he  ruby  fluorescence  was  excited  by  about  10  mVV  of  power  from  the  5145  A  line  of  an 
argon  laser,  which  has  been  found  empirically  to  be  considerably  less  than  the  power 
required  to  measurably  raise  the  temperature  ol  the  small  chips  of  rubv  typically  used 
( 10-30  pm  on  a  side).  The  R-line  spectra  were  measured  with  sufficient  resolution  that 
the  accuracy  of  the  pressure  measurements  was  limited  by  noise  and  wavelength  cali¬ 
bration  to  about  a:  I  kbar. 

The  lifetime  and  high-resolution  fluorescence  measurements  of  the  JF  ;  —  JI , :  tran¬ 
sitions  in  theNd.Yi-  ,PsOu  samples  and  of  the  ruby  transitions  were  made  with  excitation 
by  the  5145  A  line  of  the  argon  laser  focused  onto  the  sample  through  the  driving  screw 
end  of  the  pressure  cell.  Fluorescence  w'us  collected  Irom  the  opposite  end  of  the  cell, 
focused  onto  the  entrance  slit  of  a  one-meter  monochromator  equipped  with  a  grating 
blamed  at  1  tun;  and  detected  by  a  cooled  RCA  C3I034  photomultiplier  tube.  For 
analysing  the  ruby  vibronic  spectra,  a  resolution  of  about  0.4  A  was  used  and  for  the 
pentaphosphate  spectra  the  resolution  was  about  1 .2  A.  To  record  steady-state  spectra, 
the  excitation  power  was  about  It)  mW.  the  beam  was  mechanically  chopped  at  100  Hz. 
and  the  signal  analysed  by  a  lock-in  amplifier.  For  lifetime  measurements  an  electro- 
optic  shutter  was  used  to  generate  50  us  pulses  with  pulse  heights  chosen  to  keep  the 
average  power  well  below  10  mW.  The  decay  times  were  analysed  by  a  boxcar  integrator. 

Optica)  absorption  measurements  were  made  using  a  UK)  W  quartz-halogen  lamp 
through  a  1.4  M  monochromator  for  the  illumination  source.  The  light  was  detected  hv 
a  cooled  RCA  C31034  photomultiplier  tube  whose  output  was  measured  by  a  pico- 
ammeter.  The  observed  shift  in  the  ruby  R,  fluorescence  line  position  indicates  that  the 
light  source  caused  sample  temperatures  to  be  raised  by  about  5  to  Kl  E  in  these 
experiments. 

For  infrared  fluorescence  measurements,  the  -.ample  was  illuminated  by  up  to 
200  m\V  of  power  of  the  5145  A  line  of  the  argon  laser,  which  caused  less  temperature 
rise  than  the  transmission  experiment,  and  the  signal  was  chopped  at  500  Hz.  The 
fluorescence  was  detected  by  a  cooled  PbS  detector  whose  output  was  analysed  by  a 
lock-in  amplifier. 

The  samples  used  in  these  measurements  were  chips  broken  from  single  crvstals. 
The  rubv  was  analysed  by  atomic  absorption  and  found  to  have  about  ().03ri  Cr  *  ions 
replacing  the  aluminium.  The  Nd.i  ;Y  nP:0,j  and  NdP<Ou  crystals  cleave  easily  and  it 
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was  thus  possible  to  choose  flat  pieces  with  parallel  faces  lor  the  optical  transmission 
experiments.  However,  it  was  difficult  to  measure  the  thicknesses  of  these  small  pieces 
accurate!) .  so  only  relative  absorption  coefficients  of  a  given  piece  are  reported  at  the 
various  pressures 


3.  Nd,Y|  .  ,P<Ou  results  and  interpretation 

The  fluorescence  decay  times  of  the  Nd1'  emission  in  the  two  pentaphosphate  samples 
were  measured  at  several  pressures  up  to  about  ‘>0  kbar  The  fluorescence  decay  rates 
are  plotted  in  figure  I  The  decav  curves  for  NdP.(_).4  were  simple  exponentials  at  all 
pressures.  At  low  pressures,  the  decay  curves  lor  Nd„  j  were  also  pure  expo¬ 

nentials  but  at  higher  pressure  a  taster  initial  decav  was  observed  to  grow  m  I:  was 
difficult  to  get  a  quantitative  value  for  this  initial  decay  time,  since  it  was  not  greatlv 
different  from  the  asy  mptotic  time  and  lasted  for  only  about  the  tiist  m)  us  of  the  decav 


Figure  I.  I  ota!  fluorescence  decay  rate  of  *F.:  level  in  \d,Y. .  ,iM  >  4  .or  c  =  >>  j  jnd  l  U.  .is 
a  function  of  pressure 


i  r  i  i  i  r  i - 1  i  i 


P  krc.' 


Figure  2  shows  the  fluorescence  spectra  of  the  JF-. ;  — •'‘R^  Nd  "  transitions  at  low 
and  high  pressures  for  both  pentaphosphate  samples.  The  shifts  in  positions,  changes  :n 
relative  intensities,  and  variations  in  lineshapes  with  pressure  are  different  for  the  two 
samples.  Figure  3  shows  plots  of  the  line  positions  versus  pressure  and  the  pressure- 
dependent  w  idths  of  two  of  the  transitions  are  shown  in  figure  4.  These  are  the  only  tw  o 
lines  well  enough  resolved  to  allow  their  widths  to  be  measured  accurately.  The  labels 
for  these  transitions  are  R|  and  R;  for  the  lower-  and  higher-energy  crystal-tield  com¬ 
ponents  of  the  *F-. ;  metastable  state  and  Z,  refers  to  the  ground  state.  The  shape  of  the 
Rr-Z|  transition  in  NdP«Ou  was  analysed  to  determine  the  homogeneous  (lorentzian  I 
contribution  to  the  linewidth  and  this  is  also  shown  in  figure  4.  All  of  the  lines  in  the 
spectra  of  NdP<Ol4  appeared  to  narrow  at  high  pressures  similarly  to  the  data  show  n  in 
figure  4.  The  equivalent  transitions  in  Nd„  . Y0 .jP.Ou  all  broadened  considerably  with 
pressure,  but  they  were  not  well  enough  resolved  to  make  quantitative  measurements 
of  the  widths  as  can  be  seen  in  figure  2. 
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K:;;ure  4.  Pressure  JcnenJciici:  >>!  the  R.  •  /.  . . .  i.  X.  —  '  lot  M  l  jnd  X  —  /. ;  I hoin  i 

1  •  '  l.ni/'s’illlv,  ui  Ndl'.O  . 


The  line  positions  tor  SdP«Ou  .»re  shown  in  table  !  at  low  and  high  pressures  These 
aie  compiled  from  \  isthle  absorption  and  fluorescence  spectra  and  from  infrared  emis¬ 
sion  spectra.  The  fluorescences  to  the '!.? :  ano  *1;, :  manifolds  were  measured  with  \ ery 
w  ide  slits  tin  the  ntonochromator  oue  to  the  poor  sensitivity  of  the  I’bS  detector  and  thus 
oniy  the  total  breadths  of  these  manifolds  are  given.  Transitions  to  the  JIi< :  levels  were 
too  wear.  to  observe  either  in  fluorescence  or  absorption  tor  a  sample  small  enough  for 
diamond  cell  studies  and  the  levels  listed  in  table  1  are  taken  from  Blatte  etui 1 1973;  and 
st'.v'wn  tor  completeness.  Note  that  most  of  the  levels  move  to  lower  energies  with 
increasing  pressure  and  some  appear  to  be  more  sensitive  to  pressure  than  others. 

The  line  broadening  with  pressure  observed  in  the  fluorescence  spectrum  ot  the 
unfitly  doped  . ample  and  the  line  narrowing  reported  in  tic. ire  c  tor  the  neaviiy  doped 
sample  are  ditficult  to  interpret.  In  other  materials,  such  as  Y.Ai-Q,;.NJ'‘ .  the  homo¬ 
geneous  contribution  to  the  iinewidths  has  neen  attributed  to  phonon  absorption  tr.d 
emission  transitions  among  the  various  levels  of  each  tnamfoid  t  Kushida  i'/fv-'i  Those 
processes  may  change  with  pressure  as  the  positions  of  the  energy  levels  change,  but 
exact  calculation  of  such  changes  is  not  possible  without  data  on  the  pressure  v  ariation 
ot  the  vihrntu.nai  properties  of  these  crystals. 

The  ohserv  ed  increase  in  the  liuorescence  deea;.  rate  at  high  pressures  oouid  possibly 
oe  due  to  one  ot  three  effects:  an  increase  in  tile  radiative  decay  rate  ot  each  Nd"  ion: 
in  increase  in  the  non-radiative  decay  rate  of  each  Nd"  ion:  or  an  increase  in  the 
interaction  between  Nd  ’  ions  leading  to  liuorescence  quenching.  The  third  possibility 
appears  to  be  the  most  likely .  since  single-ion  effects  .should  be  independent  of  concen¬ 
tration  and  die  observed  pressure-dependent  quenching  properties  are  different  for 
lightly  doped  and  he.iviiy  doped  samples,  in  addition,  quantitative  estimates  ol  the 
non-iadiative  decay  rate  of  the  T  :  level  give  values  of  cU.2  of  the  total  decay  rate 
l  Weber  a  al  l‘C4.  Powell  a  ai  ll'M)bj  and  it  is  doubtful  that  pressure  would  cause  this 
to  merea-e  to  the  extent  necessary  to  explain  the  NdP  Ou  data. 

idle  radiative  decay  rates  can  be  estimated  irom  the  absorption  strength  of  the 
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4I., .  s .  transitions  and  the  relative  fluorescence  intensities.  If  the  transition  proba¬ 
bilities  HA  Kl  and  lV,rR.  are  known,  the  emission  rate  is  found  from  Kinston's  relations 
to  be 

If  it-/.  ^  hr.VVr,./,  +  HrTVh-  /,  *  tijt.VV'/.,  r.  +  hr, IV... .«  l  ' ) 

where  hr,  and  nR  are  the  tltermal  population  factors  for  the  two  crystal-field  components 
of  the  JF ; :  level.  By  working  with  the  ratios  of  transition  rates  for  high  and  low  pressures, 
problems  with  accurately  measuring  sample  thickness  can  be  avoided  in  determining 
the  changes  in  absorption  cross  sections.  Correction  must  be  made  to  account  for  the 
change  in  the  number  of  ions  in  the  path  of  the  light  beam  This  can  be  accounted  for  by 
a  tactor  of  inb.S  kbar !  a~(()  kbar)  where  a  represents  the  lattice  constant.  Table  2  lists 
the  integrated  absorption  coefficients  for  the  two  transitions  and  pressures  of  interest. 


Fable  2.  Relative  radiative  transition  rate  parameter*  lor  )■.*  at  iinm-plienc  pro* 'arc 
and  hS  kbar. 
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It  also  gtves  the  relatis e  thermal  populations  of  two  crystal -field  components  ot  the  ’F- ; 
!e\el  and  the  computed  radiative  decay  rate  to  the  Z,  level.  Comparison  ol  the  relative 
fluorescence  intensities  of  transitions  to  other  components  of  the  '!,•  term  and  correction 
for  the  known  branching  ratios  at  1  bar  i  Lomheim  and  DeSha/.er  I97.S).  noting  that  the 
transitions  to  the  levels  make  a  negligible  contribution  to  the  decay  rate,  gives  the 
results 
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A  recent  paper  by  Asaumi  cr  ai  1 10, SO)  gives  the  compressibility  of  NdP-O.j  along  the  a 
and  c  axes  but  not  the  b  axis.  Their  data  indicate  that  tile  compression  factors  in  the 
above  equations  must  be  no  smaller  than  about  0.95.  The  radiance  decay  rate  it 
atmospheric  pressure  has  been  estimated  to  be  about  3225  s~‘  and  thus  the  increase  vith 
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pressure  predicted  by  equation  (31  is  far  too  small  to  account  tor  the  NdP--0:4  data  shown 
in  figure  I. 

It  has  been  shown  that  the  oscillator  strengths  tor  Nd1'  transitions  in  Nd.La, 
P.Ou  are  independent  of  Nd5’  concentration  (Auzel  197b).  Thus  the  change  in  the 
radiative  decay  rate  with  pressure  should  also  depend  little  on  concentration  and  tile- 
prediction  of  equation  (3 )  should  be  applicable  to  the  Nd,, ,  V,  ..ITOn  sample.  Comparing 
equation  (3)  with  the  data  for  the  lightly  doped  sample  shown  in  figure  1  shows  relatively 
good  agreement. 

The  above  considerations  indicate  that  the  pressure  dependence  of  the  fluorescence 
decav  rate  in  NdP-Ou  is  not  due  to  changes  in  single-ion  decay  rates.  Two  multi-ion 
quenching  processes  are  possible,  cross-relaxation  between  pairs  of  Nd’’  ions  (Strek  el 
ai  !l.C7 1  and  energy  migration  to  sinks!  I.empicki  W77.  Hubert  v  .md  Powell  1479.  Powell 
el  al  19,X0a).  for  electric  dipole-dipole  interaction,  both  ot  these  processes  have  an 
R-*  dependence  on  the  separation  between  Nd"  ions  and  this  will  vary  only  weakly 
with  pressure  because  of  the  small  co  r.pressibility  ot  the  material  ( Asatnni  ct  al  IPSO). 
The  efficiency  ot  cross-relaxation  between  pairs  ot  Nd"  ions  is  dependent  upon  the 
degree  of  resonance  between  the  'F> ;  —  'l|- ;  transitions  on  one  ion  ami  the  ;  —  ; 

transitions  on  the  other  ion.  The  branching  ratio  for  fluorescence  to  the  *1 :5  -  manifold 
is  quite  small  at  atmospheric  pressure  and  the  fact  that  these  transitions  are  still  not 
i  ibserv  ed  at  high  pressure  suggests  that  the  matrix  elements  for  these  transitions  remain 
smail.  Since  transitions  to  the  41|<;  levels  could  not  be  observed,  the  change  in  their 
positions  with  pressure  is  unknown.  However,  since  the  'F-. ;  levels  move  to  lower  energy 
with  increased  pressure  and  the  spectral  lines  become  narrower,  in  order  to  have 
increased  resonance  for  cross-relaxation  transitions  at  high  pressures  the  'I,?;  levels 
would  have  to  move  very  strongly  to  lower  energies  and  reach  very  specific'  positions 
favourable  for  energy  matching.  Thus  the  observed  spectral  properties  do  not  appear  to 
favour  strongly  increased  cross-relaxation  at  the  pressures  attained  in  this  study. 

The  linal  possibility  is  that  the  increased  decay  rate  at  high  pressures  is  due  to  the 
increased  energy  migration  to  traps.  This  process  involves  transitions  between  the  JF;,; 
and  JL :  states  for  which  the  matrix  elements  were  shown  above  to  increase  with  pressure. 
Also,  since  the  transfer  involves  ions  whose  relevant  energy  levels  ;fre  in  resonance,  the 
narrowing  of  the  fluorescence  lines  increases  the  interaction  strength  at  high  pressures. 
A  simple  diffusion  mode!  may  be  used  to  check  the  plausibility  of  this  hypothesis 
quantitatively .  In  this  ease  the  rate  of  energy  transfer  to  traps  is  proportional  to  the 
diffusion  coefficient.  D.  and  for  resonant  electric  dipole-dipole  interaction  this  is  given 
by  (Trlifaj  1  V5s ) 

\VmK  *  D  x  a  A;V.!m  (4! 

where  a  is  an  average  lattice  spacing.  H'lon  is  the  single-ion  transition  rate  and  A;\Jra  is 
the  homogeneous  linewidth  which  is  used  to  approximate  the  -  pectral  overlap  integral 
between  the  two  interacting  ions.  Using  equal  ion  i  2 1.  and  the  homogeneous  line  widths 
measured  for  the  R -Z,  line,  the  ratio  of  the  rates  of  energy  migration  and  trapping  at 
high  and  low  pressures  is 

VVr"B(ft,X  kbar)  _  aJ(i)  kbar)  ,a:(b.X  kbar)  ,  ___  ()  ,;;35\  =  ,4ll 

VVn"*(0  kbar)  a4(68  kbar)  '  <r(()  kbar)  ~  ""  -21' 

This  prediction  can  be  compared  with  the  ratio  of  the  radiationless  quenching  rates 
found  from  the  experimental  data  shown  in  ligure  1  and  corrected  for  the  change  in 
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radiative  decay  rate  using  equation  (3) 

lUVSkbar)  .  tr^bSkbar)  -  IV^bS  kbar)  .  ,  _  ,  .rfftXkbar) 

H’^lOkbar)  U^"f()  kbar)  -  U'-'iO  kbar)  ~  (<U  *  ,U)  7, -(77k bar) 

where  the  radiative  decay  rate  at  atmospheric  prcssute  lias  been  estimated  to  be 
.-22h  s  As  noted  earlier,  the  compression  laclor  is  nearly  unity,  so  that  the  observed 
ratio  of  the  radiationless  quenching  rates  is  about  2.2  r  0.3,  which  is  in  close  agreement 
with  the  prediction  ot  equation  [5). 


5.  Discussion  and  conclusions 

The  physical  mechanisms  leading  to  the  observed  changes  in  the  energy  level  positions 
and  transition  strengths  with  increased  pressure  in  the  neodymium  pcniaphosphate 
samples  are  probanly  due  to  changes  in  the  symmetry  of  the  crystal  held  and  the  amount 
ot '  configuration  mixing  due  to  the -strength  ot  the  odd  crystal -Held  components,  Further 
studies  are  necessary  to  clarity  the  details  ot  these  changes.  Di  tie  rent  transition  energies 
have  been  associated  with  several  types  ot  sue  symmetries  in  Nd,V.  J'.U;;  crvstais 
i  Kruhlci  ft  at  Id"? )  but  no  consistent  trend  toward  one  tv  pe  ot  symmetry  can  ne  seer,  in 
our  nigh-pressure  results.  Some  work  on  crystal-held  changes  with  pressure  has  been 
reported  by  i  iuber  ct  al  t  W)  on  europium  pentapnospiiate. 

For  the  lightly  doped  mixed  pcniaphosphate  sample  the  decreased  Inctime  with 
pressure  is  consistent  with  the  increased  radiative  decay  rate  whereas  the  -trongcr 
variation  observed  in  the  heavily  doped  sample  can  only  be  explained  by  increased 
concentration  quenching.  The  v  ariationsof  the  quenching  rate  and  transition  linewiutns 
with  hydrostatic  pressure  are  more  consistent  with  a  inodei  oi energy  migration  quench¬ 
ing  than  one  ot  cross-relaxation.  However,  they  do  not  help  to  distinguish  between 
randomly  distributed  quenching  sites  and  surface  quenching. 
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Transient  four-wave  mixing  techniques  .'.a---  teen  me™.-,  .r  ■  ...-  r  r.c.e  •  -1 

•  -...-  excited-state  copulation  grating  in  N-i,  La: P.U:<  .-  as  a  :  gr  ;  .g 

s[acini;.  The  results  show  tl-.at  at  room  ten,;. .-rat  a  re  there  till™, i-  ■  -.-.igr.ition 

jser  distances  of  the  order  of  0.3(1  pm  :n  the  highit  •  •■'•  c-r.t.- ..  ,e:.. 


PACS  lumbers:  .  s.PO.-e 

D'-ger.erate  four-wave  mixing  (FWM)  .spectros¬ 
copy  ha.-:  beers  shown  recently  to  be  a  powerful 
.method  for  studying  spatial  migration  without 
spectral  transfer  of  electronic  excitation  energy 
in  solids. •'*  This  technique  has  been  used  to  de¬ 
termine  the  diffusion  coefficient  of  molecular  ex- 
•  ::;•  ns  in  organic  solids,'  but  attempts  to  make 
similar  measurements  in  Frenkel  excitons  in  in¬ 
organic  materials  have  resulted  mly  in  placing  an 
upper  bound  on  the  diffusion  coefficient  because 
the  migration  is  too  short  to  ibserve  in  the  sam¬ 
ples  which  have  been  studied.-"’  We  report  here 
the  results  of  FWM  measurements  on  single  crys¬ 
tals  of  Nd,  Li..,  P.0,,  at  room  temperature.  The 
presence  of  energy  migration  is  easily  observed 
and  :t  is  found  to  be  diffusive  with  a  migration 
iis'ance  of  the  order  of  0.3^  ,-m.  These  results 
..re  of  significant  importance  since  they  represent 
an  unambiguous  determination  of  rhe  range  of  en¬ 
ergy  migration  in  this  class  of  stoichiometric 
materials  used  as  “minilasers”  for  low-thresh¬ 
old,  high-gain  applications. The  mechanism 
causin'  .■ur.gor.tratiun  quenching  of  the  fluores¬ 
cence  in  these  materials  is  .not  understood  but  it 
is  -mown  to  !...vc-  quite  different  properties  from 
other  neodymium  laser  materials.5  This  is  not 
only  an  interesting  fundamental  physics  question, 
but  it  :s  aiso  .rr.portant  in  material  development 
f  r  I;.  .-r  applications.  The  com  ribution  of  ener¬ 


gy  migration  to  :•  nwn!  ration  quenching  .:.  :ie  : - 
(iymium  ncr.tapiiusphate  has  beer,  the  subject  of 
significant  .merest  and  some  controversy,  with 
estimates  of  migration  lengths  ranging  from  a 
few  angstroms  to  a  few  microns.7'1'  These  FWM 
results  provide  an  answer  to  one  of  the  f™nda.-.,en- 
tal  questions  surrounding  this  controversy. 

The  samples  used,  for  tins  investigation  were 
cleaved  from  high-quality  single  crystals  grown 
by  the  technique  described  recently.15  Several 
different  experimental  arrangements  have  been 
used  to  establish  and  probe  population  gratings  ;f 
excited  states.  A  schematic  diagram  -,f  the  setup 
used  for  this  work  is  shown  in  Fig.  1.  The  51-45 -A 
line  of  an  argon  laser  was  used  since  it  falls  on 
one  edge  of  a  Nd5*  absorption  line.  The  absorp¬ 
tion  of  this  light  by  the  Nd1*  ions  in  the  -'ample 
creates  a  spatial  distribution  of  excited  states 
with  a  sinusoidal  pattern  of  wavelength  -  *  u. 
where  :•  is  the  crossing  angle  of  the  write  beams 
and  A  is  the  wavelength  of  the  light  in  the  .-.ample. 
This  "population  grating”  can  ''ause  the  probe 
beam  to  be  scattered  with  the  maximum  -c.rti  r- 
ing  efficiency  occurring  when  the  P.racg  condition 
is  satisfied. 

Figure  2  shows  the  results  from  m.t  .is  .ring  the 
decay  rate  K  of  the  scattered  probe  beam  ,.s  a 
function  of  the  square  of  the  crossing  angle  >f  the 
write  beams  for  single  crystals  of  N.iP  0_,  and 
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KIG.  1.  Ex pen menial  stitup  for  transient  four-wave 
mixing  measurements.  1  and  2,  the  write  beams;  3. 
the  incident  probe  beam;  and  4,  the  scattered  probe 
beam . 

Nd^  Xa^ jP.,0,,,  at  room  temperature.  The  decay 
curves  were  found  to  be  exponential  and  the  data 
points  shown  in  the  figure  represent  the  average 
obtained  from  repeating  the  measurements  nu¬ 
merous  times  at  the  same  crossing  angle.  The 
error  bars  represent  the  spread  in  the  data  ac¬ 
quired  at  each  point.  The  increase  in  background 
“noise"  due  to  extraneous  light  scattering  and 
the  decrease  in  probe  beam  transmission  caused 
the  data  obtained  on  the  heavily  concentrated 
sample  to  be  less  accurate  than  that  obtained  on 
the  20*  N'd  sample. 

For  both  of  the  samples  investigated,  the  probe 
beam  decay  rate  increases  approximately  linearly 
with  increasing  values  of  This  can  be  related 
directly  to  the  mechanisms  responsible  for  the 
destruction  of  the  population  grating.  This  grating 
can  be  destroyed  both  by  the  decrease  in  the  ex¬ 
cited-state  population  by  normal  fluorescence 
decay  and  by  the  migration  of  excited  states  (ex- 
citons)  from  the  peak  to  the  valley  regions  of  the 
grating.  If  the  exciton  motion  is  diffusive  the 
time  evolution  of  the  spatial  distribution  of  the 
excited-state  population  rc(x,  /)  is  given  by  the 
one-dimensional  diffusion  equation1 

3n(x,  l  )/!>t  =  D32n(x,  t  )/'d  x3  -  T''n(x,  t),  (1) 

where  D  is  the  diffusion  coefficient  and  ~  is  the 
fluorescence  decay  time.  For  the  initial  condition 
of  a  sinusoidal  spatial  distribution  of  n(x,  0)  in  the 
x  direction,  the  solution  of  Eq.  (1)  is 

ufx,  t )-  i  e' ,/r  {l  ♦  exp{-  (2x  /.\)2D/jcos(2rrx <A)}. 

(2) 


A(Prn) 


FIG.  2.  Measured  decay  rales  of  the  scattered  probe 
beam  for  NdP<014  and  Nd0- ; L.a/j,  sP^O:i  crystals  at  room 
temperature  for  various  values  of  the  crossing  angle  of 
the  write  beams.  The  filled  points  are  twice  '.he  meas¬ 
ured  fluorescence  decay  rales  and  the  straight  lines 
represent  the  predictions  of  Eq.  i-tl . 

The  scattering  efficiency  is  proportional  to  the 
depth  of  the  grating. I7,1“  This  leads  to  the  ex¬ 
pression  for  the  diffracted  probe-beam  intensity 

=  (3) 

where 

K  =  2[(2rt/  A)2D  *  r ' 1  j  =(8jr£>/A2}t?:  -  2t"‘.  (4) 

Here  !f  is  the  probe  beam  intensity  in  the  absence 
of  scattering  and  Jt  is  the  write  beam  intensity. 
Both  X  and  0  have  been  corrected  for  the  index  of 
refraction  of  the  sample. 

The  theory  outlined  above  predicts  that  in  the 
presence  of  diffusive  e.xciton  migration  the  meas¬ 
ured  scattered  probe-beam  intensity  will  decay 
exponentially  with  a  decay  rate  which  increases 
linearly  with  o2.  These  predictions  are  consis¬ 
tent  with  the  experimental  results  obtained  on  the 
NdjLa^PjO,,  crystals.  The  solid  lines  in  Fig.  2 
represent  the  fit  of  the  theoretical  expression  in 
Eq.  (“5)  to  the  experimental  data.  The  intercepts 
of  these  lines  should  be  2/r  and  their  slopes  can 
be  used  to  determine  values  of  1).  Table  I  lists 
the  measured  fluorescence  lifetimes  for  these 
samples  and  they  are  used  to  plot  the  solid  points 
shown  in  Fig.  2  at  £i2  =  0.  The  theoretical  lines 
extrapolate  nicely  to  these  points.  Table  1  also 
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table  i 

La  1  -  *  P iOu 

.  Exciton  diffusion 
crystals. 

parameters  for  Nd,  - 

Sample  U) 

i  (>jsec) 

O  (cm3/ sec) 

lx  tpm) 

0.2 

300 

6.85*  10"’ 

0.20 

1 .0 

125.3 

5.09*  10"s 

0.3C 

lists  the  calculated  diffusion  coefficients  for  these 
samples.  The  average  displacement  of  an  exciton 
along  the  grating  directions  in  its  lifetime  is  giv¬ 
en  by  Ii  =  [2DtY/2.  The  values  for  this  one-dimen- 
siona!  exciton  diffusion  length  for  the  two  penta- 
phosphate  samples  are  listed  in  Table  I. 

The  results  reported  here  were  all  obtained 
with  the  laser  beams  incident  on  the  6  plane 
(cleavage  plane)  of  the  crystal  and  with  the  sam¬ 
ple  oriented  to  measure  diffusion  in  approximate¬ 
ly  the  a  direction.  Although  some  changes  in  the 
results  were  observed  when  the  sample  was  ro¬ 
tated  to  different  orientations,  this  experiment 
is  so  sensitive  to  precise  alignment  that  it  is  not 
clear  at  this  point  whether  the  observed  changes 
are  due  to  anisotropy  in  the  exciton  diffusion  or 
simply  a  slight  misaligning  of  the  experiment. 

The  question  of  orientation  dependence  of  D  will 
be  thoroughly  investigated  in  the  future.  At  high 
"probe-beam"  intensities  a  double  exponential 
decay  was  observed  with  a  faster  initial  decay 
rate.  At  angles  larger  than  those  reported  in 
Fig.  2  this  fast-decay  component  dominated  the 
delected  signal  even  at  smaller  probe-beam  in¬ 
tensities  and  thus  limited  the  range  of  grating 
spacings  available  for  study.  The  decay  rate  of 
the  fast  signal  is  independent  of  the  crossing  an¬ 
gle  of  the  "write  beams"  and  does  not  extrapo¬ 
late  to  2/t  and  thus  is  not  directly  associated 
with  the  "population  grating."  There  are  several 
other  sources  of  refractive  index  gratings  in 
crystals  due  to  different  types  of  nonlinear  opti¬ 
cal  .nteractions  and  it  is  not  surprising  that  other 
interference  effects  are  present  in  a  birefringent, 
ferroelastic  crystal  such  as  Nd,  La^P^O,,.  In 
this  paper  we  focus  our  attention  on  the  proper¬ 
ties  of  the  population  grating  and  plan  to  study 
the  characteristics  of  other  types  of  gratings  in 
future  investigations.  It  should  be  pointed  out 
that  thermal  gratings  can  also  be  studied  by  this 
technique  and  that  the  thermal  diffusion  coeffi¬ 
cients  are  of  the  same  order  of  magnitude  as 
those  reported  here.  We  can  be  sure  that  the  dif¬ 
fusion  coefficients  reported  are  not  due  to  ther¬ 
mal  diffusion  since  the  observed  decay  rates  ex¬ 
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trapolate  to  twice  the  fluorescence  decay  rates 
and  the  value  of  D  varies  significantly  with  Nd 
concentration.  Neither  of  these  characteristics 
is  true  for  thermal  gratings. 

One  important  aspect  of  characterizing  energy 
migration  is  identifying  the  microscopic  nature 
of  the  interaction  causing  the  motion.  A  rough 
approximation  for  the  diffusion  coefficient  of  an 
exciton  undergoing  an  incoherent  random  walk  re¬ 
sulting  from  electric  dipole-dipole  interaction  is1* 

0=t(t"'VJ4/3*oVto.  (5) 

where  .V,  is  the  sensitizer  concentration,  t.  is 
the  intrinsic  fluorescence  lifetime,  ar.d  /<„  is  the 
"critical  interaction  distance."  For  NdP.O,.,, 

=  4*1031  cm"3  and,  if  concentration  quenching  is 
attributed  to  energy  migration,  ro  =  3  50  Use 
of  these  values  and  the  measured  value  of  1J  in 
Eq.  (5)  gives  ft0=4 5  A.  Although  this  is  extremely 
high,  compared  with  typical  values  for  Ra  in  light¬ 
ly  doped  solids,  it  is  consistent  with  the  theoreti¬ 
cally  predicted  value  fl0  in  a  stoichiometric  crys¬ 
tal  such  as  NdPsOM  if  resonant  interaction  be¬ 
tween  ions  in  identical  sites  with  small  homoge¬ 
neous  linewidlhs  is  assumed.11  Also,  Eq.  (5)  pre¬ 
dicts  that  the  diffusion  coefficient  should  vary  as 
the  y  power  of  the  concentration,  which  is  con¬ 
sistent  with  the  results  given  in  Table  I.  In  addi¬ 
tion,  it  should  be  noted  that  the  electric  dipcle- 
dipole  interaction  lias  a  long  enough  range  to  over¬ 
come  any  structural  anis  nropyU 

These  FWM  results  show  that  at  room,  tempera¬ 
ture  in  NdP^O,,,  excitor.s  diffuse  over  an  average 
distance  of  0.36  ,-m  in  a  specific  direction.  The 
relatively  long  range  of  the  exciton  diffusion 
measured  here  disagrees  strongly  with  se'.  ■  ral  of 
the  previously  published  results.9,  “■  15  However, 
it  should  be  noted  that  the  authors  of  Ref.  15  have 
obtained  more  recent  results  in  agreement  with 
those  reported  here.5  It  is  not  clear  why  the  re¬ 
sults  of  Weber  and  co-workers9’  11  are  in  contra¬ 
diction  to  the  present  results,  but  one  source  of 
difference  may  be  the  greatly  improved  quality  of 
the  large-size  crystals  now  available  for  study.19 
Although  the  results  reported  here  do  not  give 
direct  information  concerning  the  mechanism  of 
concentration  quenching,  they  do  show  the  pres¬ 
ence  of  long-range  energy  migration  and  thus  it 
is  quitepossible  for  exciton  diffusion  and  trapping 
to  play  an  important  role  in  the  quenching  pro¬ 
cess.  The  complete  theoretical  characterization 
of  energy  transfer  in  Nd,La1.1P50M  must  await 
further  experimental  results  and  this  investica- 
tion  is  presently  being  extended  to  several  other 
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samples  with  different  Nd3’  concentrations  in  or¬ 
der  to  establish  the  concentration  dependence  of 
the  diffusion  parameters,  which  can  then  be  com¬ 
pared  to  the  variation  of  the  quenching  rate.  AJso, 
these  experiments  will  be  extended  to  low  tem¬ 
peratures  to  see  whether  any  coherent  contribu¬ 
tion  to  the  exciton  migration  can  be  detected.20 

This  work  was  supported  by  the  U.  S.  Army  Re¬ 
search  Office. 
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5/-Eiectron  Localization  in  Uranium  Compounds 

Wolf-Dieter  Schneider  and  Clemens  Laubschat 
Institut  IVY  Atom  -  und  Fc  stkorperphysiM .  Freic  I'uivcrsitdt  Berlin,  D-1000  Berlin  .73,  Cermany 

(Received  17  December  1980) 

Observed  7-eV  satellites  In  the  x-ray-photoemission  spectra  of  UGa,  and  iher  B- 
group  compounds  are  shown  to  be  due  to  two-bole  final  states  as  cor.fir med  by  existing 
Auger  data.  The  presence  of  Uiese  satellites  is  an  indication  for  a  weak  id  hybridization 
and,  when  compared  to  uranium— transition-metal  compounds,  increased  5/  localization. 


PACS  numbers.-  71. 70. Ms,  79. GO. On 

The  nature  of  the  5/  electrons  in  the  actinides 
his  attracted  much  attention  both  from  the  exper¬ 
imental  and  the  theoretical  point  of  view.  For  ar¬ 
il  metal  it  has  been  shown  recently  that  the  5/ 
electrons  have  primarily  itinerant  character.1 
The  degree  of  5/  localization  depends  on  the  over¬ 
lap  of  the  corresponding  5/  wave  functions  on 
neighboring  atoms  and  on  the  5/ -6 d  hybridization. 
Since  in  uranium  compounds  the  interatomic  U-U 
distance  and  the  bonding  properties  are  changed 
us  compared  to  or-U,  a  change  in  the  degree  of 
the  localization  can  be  expected. 

X-ray  photoemission  spectroscopic  (XPS)  in¬ 
n-stir. Hums  a)  uranium  compounds  with  77-group 


elements  have  shown  a  characteristic  satellite 
structure  at  7  eV  higher  binding  energy  of  the  U 
4/  core-level  spectra  (see  Table  If.  The  only 
serious  attempt  to  explain  these  structures  has 
been  made  for  UO,  where  this  satellite  has  been 
attributed  to  a  shakeup  process  from  the  oxuie- 
derived  p  band  to  unoccupied  5/  states.1,3 

In  this  Letter  we  show  that,  in  the  XPS  spectra 
of  the  intermetallic  compount  UGa^,  a  similar 
strurture  exists  for  the  U  4/  as  well  as  for  the 
U  5/  levels  (we  note  that  even  the  valence  bands 
of  UO,  show  such  a  structure,  which  has  been 
interpreted  as  belonging  to  the  O  Ip  band3).  In 
UGa,  this  structure  cannot  be  explained  by  a 
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IV.  ENERGY  TRANSFER  IN  SEVERAL  NEODYMIUM  DOPED  LASER  MATERIALS 

The  following  three  sections  describe  the  results  of  studies  of  energy 

A 

transfer  in  three  different  types  of  NdJ+  doped  laser  materials.  Section 
IV. 1  describes  the  work  done  on  YVO^  crystals  including  both  host  sensitized 
energy  transfer  and  site-selection  spectroscopy  studies.  The  manuscript  in 
Section  IV. 2  describes  energy  migration  amoung  NdJ"  ions  in  neodymium  aluminum 
borate  crystals  using  time-resolved  site-selection  spectroscopy  techniques. 

The  manuscript  in  Section  IV. 3  discusses  the  results  obtained  on  energy  migra¬ 
tion  among  neodymium  ions  in  four  different  types  of  glass  hosts.  Of  special 
interest  are  the  results  obtained  as  a  function  of  temperature  in  this  latter 
study. 


IV. 1 

Energy  transfer  processes  in  YV04:Nd3  + 

Dhiraj  Sardar  and  Richard  C.  Powell 

Physics  Department.  Oklahoma  Stale  University,  Stillwater.  Oklahoma  74074 
(Received  1 1  January  1980;  accepted  for  publication  31  January  1980) 

Laser-excited,  time-resolved  spectroscopy  techniques  were  used  to  investigate  the  characteristics 
of  energy  transfer  in  YV04:Nd  J  *  crystals.  The  temperatuare  and  time  dependencies  of  the 
energy  transfer  rates  were  determined  for  both  host-sensitized  and  impurity  ion  interaction 
processes.  The  former  process  is  found  to  have  characteristics  consistent  with  the  migration  and 
trapping  of  localized  host  excitons  whereas  the  latter  process  is  consistent  with  multistep  energy 
migration  via  a  two-phonon-assisted  electric  dipole-dipole  interaction  mechanism. 

PACS  numbers:  78.20.  -  e,  78.40.  -  q,  78.50.Ec 


1.  INTRODUCTION 

Yttrium  orthovanadate  (YVOJ  is  an  important  materi¬ 
al  for  applications  in  various  ty  pes  of  optical  systems.  It  is 
useful  as  an  infrared  polarizer,  and  a  host  for  rare-earth 
phosphors  and  lasers.1'5  In  fact,  YV04:Nd3  *  laser  crystals 
have  been  found  to  have  some  performance  characteristics 
superior  to  those  of  standard  YAlG:Nd3  +  lasers.5  Because 
of  the  significant  potential  of  this  material  it  is  important  to 
characterize  and  understand  the  fundamental  physical  pro¬ 
cesses  underlying  its  optical  properties.  We  describe  here  the 
results  of  an  investigation  of  processes  of  energy  transfer  in 
YV04:Nd‘*  crystals. 

The  general  spectroscopic  properties  of  undoped 
YVOj,  and  YV04  doped  with  Eu  3  * ,  Er3  *  ,  and  Nd3  '  ions 
have  been  investigated.1"10  Host-sensitized  energy  transfer 
has  been  studied  in  samples  doped  with  Eu3  "  and  Eri  ' 
ions.'’  This  process  is  especially  important  in  applications  as 
a  rare-earth  phosphor  material  but  also  can  be  useful  in  en¬ 
hancing  the  pumping  of  a  laser  material.  The  process  of  ener¬ 
gy  transfer  between  impurity  ions  in  YV04  has  been  investi¬ 
gated  for  Eu  '  This  process  can  be  important  in 
concentration  quenching  in  luminescent  materials  and  in  the 
operational  characteristics  of  lasers.  We  have  characterized 
both  host-sensitized  energy  transfer  and  ion-ion  interaction 
in  YV04:Nd  '  "  and  the  results  are  compared  to  those  ob¬ 
tained  previously  on  yttrium  vanadate  samples  containing 
other  tnvalent  rare-earth  ions. 

In  addition  to  the  practical  importance  of  the  informa¬ 
tion  obtained  in  this  work,  the  results  enhance  our  general 
understanding  of  the  fundamental  proceses  of  energy  trans¬ 
fer  in  optical  materials.  The  model  used  to  interpret  the  data 
obtained  on  host-sensitized  energy  transfer  is  one  involving 
exciton  migration  and  trapping  which  has  recently  been  ap¬ 
plied6  to  inorganic  “molecular”  crystals  such  as  YV04.  The 
model  employed  to  explain  energy  transfer  between  Nd3  * 
ions  in  this  host  is  based  on  one  of  the  new  two-photon- 
assisted  interaction  processes  developed  recently  to  interpret 
high-resolution  laser  spectroscopy  data.11  In  both  of  these 
areas  only  a  few  experimental  investigations  have  been  car¬ 
ried  out  and  further  work  is  necessary  to  develop  a  complete 
understanding  of  these  energy  transfer  processes. 

I  he  single  crystal  samples  of  YV04  Nd 1  ‘  were  investi¬ 
gated.  One  had  a  doping  concentration  of  2%  (2.55  \  I0;" 
cm  )  and  the  other  3 mc  (3.82  y  10:”  cm  ’). 


The  experimental  technique  of  laser  time-resolved  spec¬ 
troscopy  was  employed  for  investigating  the  properties  of 
energy  transfer.  A  block  diagram  of  the  experimental  appa¬ 
ratus  is  shown  in  Fig.  1.  The  samples  were  mounted  in  a 
cryogenic  refrigerator  for  controlling  the  temperature.  For 
studying  host-sensitized  energy  transfer  a  pulsed  nitrogen 
laser  was  used  to  excite  the  sample  whereas  for  direct  excita¬ 
tion  of  the  Nd3  *  impurity  ions  a  tunable  dye  laser  contain¬ 
ing  rhodamine  6G  and  pumped  by  the  nitrogen  laser  was 
used.  In  both  cases  the  excitation  pulses  were  less  than  10 
nsec  in  duration.  The  half  width  of  the  dye  laser  pulse  was 
less  than  0.4  A.  The  sample  fluorescence  was  analyzed  by  a 
1-m  spectrometer,  detected  by  a  cooled  RCA  C31034  photo- 
multipler  tube,  averaged  by  a  boxcar  integrator  triggered  by 
the  laser,  and  displayed  on  a  strip-chart  recorder. 

II.  HOST-SENSITIZED  ENERGY  TRANSFER 

The  absorption  edge  of  Y'VOj  crystals  occurs  at  about 
3400  A  and  the  fluorescence  emission  appears  as  a  broad 
band  peaked  near  4500  A.  Although  some  thermal  quench- 
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f  '  Then  the  magnitude  of  the  energy  transfer  rate  and 
the  proportionality  constant  K  were  treated  as  adjustable 
parameters  and  the  theoretical  solutions  were  fitted  to  the 
FKS  data. 

It  was  found  that  at  low  temperatures  the  theoretical 
predictions  for  a  constant  energy  transfer  rate  could  not  give 
a  cond  tit  :o  the  Jala  whereas  the  solid  line  giving  the  good  til 
to  the  Jala  shown  in  Fig.  5  was  obtained  using  the  time- 
•dependent  energy  transfer  rate.  In  this  ease  the  theoretical 
expression  :s  given  by 

=  K  [exp  (  t B i,  -  B ,  )t  —  la)!  1  ;  —  - — -  ) 

[  \  (B„  -B,  )t  1  -  -rto  J\ 

(4) 

w  here 

ul  =  (5) 

This  expression  desenbes  the  rate  of  energy  transfer  by  a 
single-step  electric  dipole-dipole  interaction  where  C,  is  the 
concentration  of  Nd1  '  ions  and  R„  is  the  critical  interaction 
distance  Using  the  measured  value  of  the  host  fluorescence 
decay  time  in  the  undoped  sample  r1/,  and  the  value  of  a> 
obtained  from  fitting  the  experimental  TRS  data  at  1 1  K 
gives  a  value  of  14.6  A  for  R„.  A  theoretical  prediction  for  R„ 
can  be  obtained  from  the  expression 

R  =  [5.36  •  10  (6) 

where  6  ’  is  the  quantum  efficiency  of  the  sensitizer,  SI  is  the 
overlap  .ntegral  of  the  absorption  spectrum  of  the  neody¬ 
mium  tons  and  the  emission  spectrum  of  the  host,  n  is  the 


index  of  refraction.  vw  is  the  average  wave  number  in  the 
region  of  spectral  overlap,  and  the  numerical  factor  is  lor 
unit  consistency  and  includes  a  (actor  of  *  for  the  average 
angular  dependence  of  randomly  oriented  dipoles.  The  over¬ 
lap  integral  is  found  from  spectral  data  to  be 
SI  —  7.3  x  10  ‘  l/mol  cm  Using  this  value  along  with 
v,a  =  2.0  x  10  ‘cm  1 .  ii  —  0.54,  and  n  =  1 .53  predicts  a 
value  of  5.5  A  for  R,t. 

At  high  temperatures  (about  100  K  and  above/  it  was 
found  that  Eq.  (4)  no  longer  gave  a  good  fit  to  the  time- 
resolved  spectroscopy  data.  Instead  the  best  fit  to  the  data 
was  found  by  assuming  a  time-independent  energy  trnasfer 
rate.  In  this  ease  the  solutions  to  Eqs,  (2)  and  (31  lead  to 


/,//„  =  A" 


Ul, 


l>  1  t 


•  I  exp  [US,,  -  u  -  B,  u  \  -  \\  .  c  • 

The  dashed  line  in  Fig  5  represents  the  best  fit  to  the  data 
obtained  by  treating  u'  and  K  '  as  adjustable  parameters. 

A  constant  energy  transfer  rate  is  indicate,  e  ufa  mules- 
tep  energy  migration  process.  This  can  be  treated  through 
either  the  mathematics  of  diffusion  or  random  walk  with  the 
resulting  energy  transfer  rates  expressed  as 

<u'  =  4  ~DR.C:p  I  Sal 


or 


u'  -  y„,C,  C  I  1  ’  1  '/».  (kb: 

respectively.  Here  D  is  the  diffusion  coefficient,  R  ,  is  the 
trapping  radius  at  the  activator  site.  is  the  volume  "er 
molecule,  I,,  is  the  excitation  energy  hopping  time.  C  .4  :s 
the  capacity  of  the  random  walk  which  reflects  the  size  and 
shape  of  the  trapping  region,  an  dp  represents  t  he  ;>r<  >•  .jiiii;; 
of  the  host  excitation  energy  being  transferred  to  an  act.v  •  • 
tor  after  migration  has  occurred  tv,  u  liu.;  site  next  ’o  an 
activator.  Sicnconly  one  primary  experimental  parameter  o 
obtained.  («./),  it  is  very  difficult  to  separately  determine  the 
characteristics  of  the  migration  process  and  the  trapping 
processes.  The  best  w  av  to  do  this  is  to  compare  the  results 
obtained  from  similar  experiments  on  the  same  type  of  host 
crystal  with  different  types  of  activ  ator  traps  keeping  in 
mind  that  the  migration  properties  of  energy  transfer  should 
be  the  same  for  all  samples  with  any  differences  in  the  data 
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TABLE  I  Comparison  of  the  results  obtained  on  YV04  samples  containing 
Eu  '  ‘  .  Er  '  .  and  Nd*  ‘  ions  as  activators. 


Y  VO,.Nd 1  ‘ 

YVO.fcu'  ' 

YVO.  Er  ' 

Parameters 

(is,  lO-'Vm  ') 

<l.27x  lO’-cm 

‘)  ( 1 .27  x  ’1 

R.,  (A) 

14 

4 

2 

JA'  -m  *) 

250 

156 

1000 

{ see  ‘ )  - 

4  1  »  10" 

I.OX  I04 

2.1  10' 

c |i .  p 

:  40 ...  i  o' 

1.75  10“ 

3.68X10“ 

UK,p 

O 

0.27\  10  11 

1.32.x  10  '“ 

1  0  v  to 1 

7.90  x  10  1 

1 

'7'=  297  K 


being  attributable  to  different  trapping  characteristics  of  the 
various  activators.  Table  I  compares  the  results  obtained  on 
YVO,  samples  containing  EuJ  * ,  ErJ  '  ,  and  NdJ  '  ions  as 
activators.* 

In  order  to  compare  the  fundamental  energy  transfer 
properties  in  these  three  samples  it  ts  first  necessary  to  adjust 
for  the  different  thermal  properties  of  the  energy  transfer 
rates.  The  same  qualitative  dependence  as  shown  in  Fig.  3 
was  observed  in  the  other  two  samples  but  the  measured 
activation  energies  are  quite  different  as  shown  in  Table  1. 
This  can  be  attributed  to  thermal  effects  being  present  in 
both  the  migration  on  trapping  processes.  Previous  work  on 
undoped  YVQ,  crystals  indicates  that  host  energy  migration 
occurs  w  ith  the  activation  energy  of  about  t>25  cm  indi¬ 
cating  that  the  differences  between  this  number  and  the  val¬ 
ues  of  AE '  appearing  in  Table  1  are  due  to  the  need  for  ab¬ 
sorbing  or  emitting  phonons  in  the  trapping  process  at  the 
activator  sites  which  can  be  accounted  for  as 
-  -  p  .,  expt  •*-  AE i  ,'k  T ). 

Alter  adjusting  the  data  to  account  for  the  different 
(hernial  characteristics  of  the  trapping  process  there  are  two 
possible  ways  to  interpret  the  results.  The  first  ts  to  assume 
that  the  trapping  regions  for  the  three  types  of  activators  are 
all  the  same  size  and  the  difference  in  the  trapping  character¬ 
istics  is  Cue  to  the  factor pu.  Assuming  a  point  trapping  re¬ 
gion1'’  C  ‘  .4  i  =  0.059  or  assuming  a  trapping  radius  equal  to 
the  nearest  neighbor  V  distance  of  R ,  =  3.15  X  10  ' 
cm.  the  values  obtained  for  p,,  are  listed  in  Table  1.  In  this 
.use  the  random  walk  model  predicts  a  hopping  time  of  tH 
2.  1  •  ID  "  sec  and  the  diffusion  model  predicts  that 
.0-2  o  •  10  K  enr  sec  1  In  the  limit  of  many  steps  in  the 
random  walk  the  two  models  should  predict  the  same  re¬ 
sults.  This  can  be  checked  with  the  relationship 

t,  =  a:/[bD  ),  (9) 

w here  a  is  the  average  stepping  distance  for  the  excitation 
energy.  If  this  is  taken  to  be  the  shortest  V-  F  separation  of 
4  ~5  A  and  the  value  of  D  obtained  above  is  used,  the  hop¬ 
ping  time  :s  predicted  to  be  1.5  X  10  '  sec  which  is  consis¬ 
tent  with  the  value  obtained  from  the  random  walk  model. 
The  number  of  steps  in  the  random  walk  is  given  by 

i  -  r, ,/t,.  t  10) 

which  for  (he  Nd-doped  sample  investigated  here  is  755 

!  h-  second  possible  interpretation  of  the  data  is  that 
n ,  -  1  tor  ail  cases  but  the  size  of  the  trapping  regions  arc 
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different  for  the  three  activators.  However  pursuing  this  line 
of  reasoning  leads  to  the  conclusion  that  the  trapping  radius 
for  ErJ  '  is  of  the  order  of  1830  A  which  is  unphysically 
large.  It  also  predicts  such  a  slow  hopping  time  that  only  line 
step  could  take  place  in  a  random  walk  w  hich  is  not  consis¬ 
tent  with  the  observed  time-independent  energy  transfer 
rate.  Thus  we  conclude  that  the  first  approach  to  interpret¬ 
ing  the  data  is  the  correct  one 

From  the  analysis  of  the  data  described  above  the  diffu¬ 
sion  length  at  room  temperature  is  determined  from 

I  =  r.Dr  .\\) 

to  be  8. S  '  10  ■  cm.  The  diffusion  coefficient  can  be  ex¬ 

pressed  as  D-D,  ex  pi  -  JA’, .  kT  •  where  JA\  is  the  acti¬ 
vation  energy  for  hopping  The  value  of  D,  is  found  to  be 
5.4  .,  10  enr  sec 

III.  SITE-SELECTION  SPECTROSCOPY 

The  Nd  '  '  ions  can  lie  directly  excited  with  a  dye  laser 
using  rhodamine  t>G  dy  e.  When  this  is  done  structure  can  be 
observed  in  each  of  the  spectral  transitions  due  to  the  excita¬ 
tion  of  ions  in  different  types  of  crystal  field  sites.  This  -.true- 
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111.)  •  Rjiios  i>f  th<  integrated  fluorescence  intensities  of  sensitizer  and 
activator  Nvl  '  ions  as  a  .’’unction  of  time  after  the  laser  pulse  for 
t  SO.  Nd  ’  ,3fci  i Sec  text  for  explanation  ofthe  theoretical  lines. i 


cure  changes  with  laser  wavelength  due  to  the  selective  exci¬ 
tation  of  ions  in  different  sites  as  shown  in  Fig.  7.  Energy 
transfer  takes  place  between  Nd  1  *  ions  and  this  can  aga.n  be 
characterized  through  time-resolved  spectroscopy  tech¬ 
niques.  Figure  8  shows  a  typical  emission  spectrum  near  890 
nm  at  two  times  after  the  excitation  pulse.  For  both  transi¬ 
tions  show  n,  the  high  energy  structure  decreases  and  the  low- 
energy  structure  increases  as  a  function  of  time  after  the  laser 
pulse.  Designing  the  high-energy  structure  as  transitions  due 
to  ions  in  selectively  excited  "sensitizer”  sites  and  the  low- 
energy  structure  as  transitions  due  to  ions  in  "activator" 
sites,  energy  transfer  from  sensitizers  to  activators  can  be 
characterized  by  plotting  the  ratio  of  the  integrated  fluores¬ 
cence  intensities  as  a  function  of  time  after  the  laser  pulse. 
The  results  of  doing  thts  at  low  and  high  temperatures  are 
shown  :n  Fig.  9. 

The  solid  lines  in  Fig  resent  the  best  theoreticai 

fits  to  tlie  data  based  on  the  r..  ,ndel  shown  in  Fig.  10.  This 

is  similar  to  the  model  used  m  Sec.  II  to  analyze  the  host- 
sensitized  energy  transfer  results  except  that  direct  excita¬ 
tion  of  the  activators  is  allowed  and  no  significant  difference 
could  be  detected  in  the  fluorescence  lifetimes  of  the  ions  in 
the  sensitizer  and  activator  sites.  The  rate  equations  for  the 
excited  stale  populations  in  this  model  are  given  by 

dn  /dt  =  W,  — Bn ,  —  w..  n,,  (12) 

dn  ,,'dt  =  Wj  -*-oq,  n.  —  BiJj.  (13) 

Tuese  can  be  solved  assuming  delta-function  excitation  and 
j  specific  time  dependence  for  the  energy  transfer  rate.  It  was 
found  that  the  best  fits  to  the  data  at  all  temperatures  were 
blamed  with  time-independent  energy  transfer  rales.  For 


s  d 

FIG  10.  Rate  equation  model  for  explaining  energy  ;  ransfer  between  Nd 
ions  in  V  V()4 

this  case  the  solutions  to  Eqs.  (12)  and  ( !  3 ;  give 

4,//,  =  K"  \  [nj0)/n,(0)  -  1  ]  expltq.  t )  -  1  ,  '( 14| 
where  the  radiative  decay  rates  are  assumed  to  be  me  same 
for  sensitizers  and  activators.  The  solid  lines  in  Fig  9  repre¬ 
sent  the  best  fit  to  the  experimental  data  obtained  using  Eq. 
(14)  and  treating  cu,4  and  K  "  as  adjustable  parameters.  The 
temperature  dependence  of  the  energy  transfer  rate  can  be 
obtained  using  Eq.  (14)  and  measured  fluorescence  intensity- 
ratios  at  short  and  long  times  after  the  excitation  pulse.  The 
results  are  shown  in  Fig.  1 1 .  The  rate  appears  to  be  approxi¬ 
mately  constant  at  low  temperatures  and  then  increases  sig¬ 
nificantly  above  about  25  K.  The  activation  energy  for  this 
increase  is  found  to  be  about  15  cm  :.  This  is  essentially  the 
same  as  the  crystal  field  splitting  of  the  JF,  -  level  in  this 
system. 

The  temperature  dependence  described  above  is  indica¬ 
tive  ofa  two-phonon-assisted  energy  transfer  process  involv¬ 
ing  a  real  intermediate  level. 1 1  In  this  case  the  energy  trans- 
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FIG.  I  1  Temperature  dependence  of  the  r;ite  of  energy  iranMer  between 
Nd  ions m  ^  VO,  Nd  <30  >  tSoelc.il  lor  explanation  .  :  the  tluoteue.it 
line. ) 
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TABLE  II  \d  Nd  '  '  energy  ir»n*/er  pirmwiers. 


Parameters 

YVO.Nd’- 

(3.1t  \  I0:u  cm  ’) 

C  (CfTI  J 

3ft  ■ 

10" 

t'_  (cm 

i  8  . 

nr 

j£„  \em  ) 

1  fa 

AE  .  tern  ') 

7  9 

$  (cm  1 » 

15 

r  icmi 

4  4 

J  :  R  N,  icrtd) 

40,5 

.  io-,,j 

J  i  R  * tcrtl4) 

4,41  . 

A  10 

/N„  (v no 

3  56  ■ 

10*- 

D  I  cm  sec  *) 

1,3  • 

lu 

.J  (jUS 

*  3  « 

H)  ’ 

•  r  -  :u5  K  *5  K. 


fer  rate  is  given  by : ' 

=  7 !  bC, C, /? fi)[  i  -  expU E.JkT) ]w4 


1  4 

ex p(  —  5/kT ),  (15) 

.v here  V,  and  J.  are  the  ion-ion  interaction  strengths  involv¬ 
ing  the  initial  and  intermediate  states,  respectively,  r  is  the 
width  of  the  intermediate  state,  and  6  is  the  energy  ofthe 
phonon  needed  to  change  the  ion  from  the  initial  to  the  inter¬ 
mediate  state  In  this  case  the  intermediate  state  is  the  upper 
crystal  held  component  ofthe  4/\,%  level,  so  <>  is  15  cm  : 
The  width  of  this  level  is  measured  from  the  transition  be¬ 
tween  it  and  the  ground  state  to  be  4. 0  cm  ’  JA',„  is  mea¬ 
sured  to  be  about  '.9  cm  1  .  d  E.s  is  taken  to  be  about  one- 
thtrd  ofthe  inhomogeneous  linevvidth  ofthe  sensiu/.er  tran¬ 
sition  which  is  measured  to  be  1.6  cm  1  .  This  is  consistent 
with  previous  results  on  garnet  samples, 1 '  The  squared  ma- 
ti.v  elements  can  be  expressed  in  terms  of  Judd-Ofelt  param¬ 
eters  as  4 


J  : 
j£ ; 


J] 

7r 7 


u- 


At':. 


<■)(£-,  R  '): 


CJ  -  1H27,  -  1)  Vr 


( !6l 


I'he  summations  over  reduced  matrix  elements  have  been 
evaluated  for  the  4/',.:  to  ‘/w;  transitions  of  Nd'  *  in  YVO, 
and  are  given  bv ' 

J  '  —415  ■  10  '"/?  ■  '’em1, 

(17) 

J  =  d.J !  -  10  '  R  ~  '’em4, 

where  ihe  population  differences  between  the  "a"  and  "b  " 
a  mponents  of  (he  V,  ,  level  have  been  taken  into  account. 
W , ; h  these  considerations.  Eq.  (151  predicts  values  for  the 
energy  transfer  rate  at  different  temperatures  shown  as  a 
v  >liu  line  in  Fig.  1 1 .  This  is  in  good  agreement  with  the  val¬ 
ue.  obtained  experimentally  at  high  temperatures  but  pre¬ 
dict  smaller  values  tor  the  energy  transfer  rates  than  those 
■  ibser'-'-d  at  low  temperatures.  The  theoretical  parameters 
for  this  model  are  summarized  tn  Table  II  q. 


IV.  DISCUSSION  AND  CONCLUSIONS 

The  results  of  the  investigation  of  host -sensitized  ener¬ 
gy  transfer  in  YVO,  Nd'  '  are  consistent  with  the  model 
proposed  previously  from  studies  of  yttnum  vanadate  with 
other  types  of  trivalent  rare-earth  impurities.''  In  this  mod¬ 
el  the  host  excited  state  considered  to  be  a  charge  transfer 
transition  localized  on  a  V04  molecular  ion.  This  can  be 
treated  as  a  Frenkle  exciton  which  is  self-trapped  at  low  tem¬ 
perature  and  undergoes  thermally  activated  hopping  among 
the  other  vanadate  molecular  ions  at  high  temperatures.  Pre¬ 
vious  work  indicates  that  the  mechanism  responsible  for  the 
exciton  hopping  is  exchange  interaction  and  the  large  Stokes 
shift  due  to  lattice  relaxation  in  ihe  excited  state  is  ihe  cause 
of  ihe  self-lrappmg  at  low  temperatures.  Tabic  I  compares 
some  of  the  relevant  model  parameters  obtained  from  stud¬ 
ies  of  YV04  samples  doped  with  Eu  *  ,  Er  '  ,  and  Nd’  ' 
The  trapping  probability  per  visit  to  an  activator  site  and  the 
thermal  activation  energy  for  trapping  are  found  to  be  sig¬ 
nificantly  different  for  the  samples  with  the  three  different 
impurity  ions.  The  reason  for  the  different  trapping  prob¬ 
abilities  is  not  known  but  it  should  be  noted  that  they  in¬ 
crease  with  increasing  number  of  4 /  electrons  ana  decreas¬ 
ing  ionic  radius.  Additional  activators  in  this  host  should  be 
investigated  to  determine  the  relevance  of  this  correlation. 

At  low  temperatures  where  the  exettons  are  seif- 
trapped  the  results  indicate  that  energy  transfer  to  Nd 
ions  lakes  place  by  a  single-step  electric  dipole-dipole  inter¬ 
action  process.  This  is  consistent  w  ith  the  results  found  from 
investigating  Eu’  '  and  Er’  *  doped  samples  *  [n  the  pre¬ 
vious  cases  the  values  of  R„  which  characterize  this  process 
were  found  to  be  of  the  order  of  a  few  angstroms  and  in 
agreement  with  the  predictions  of  spectral  overlap  measure¬ 
ments.  For  the  data  described  here  for  Nd-doped  samples  R 
was  found  to  be  significantly  larger  than  for  other  actr.  alors 
and  even  larger  than  'dial  predicted  by  spectral  overlap  c-  n- 
siderattons.  The  reason  for  this  is  not  presently  understood. 

A  previous  site-selection  spectroscopy  investigation  has 
been  carried  out  on  Eu  ’  ions  tn  YVO,  crystals.’  tn  that 
case  it  was  found  that  the  results  w  ere  consistent  with  energy 
transfer  by  a  single-step  resonant  electric  dipole-dip. >ie  inter¬ 
action  process  but  the  data  analysis  was  complicated  by  hav¬ 
ing  three  distinct  'vpes  of  sites.  The  results  obtained  here  are 
consistent  with  a  mulusiep  diffusion  type  of  energy  transfer 
with  the  mechanisms  for  diffutions  an-1  "'•’oping  both  in¬ 
volving  two-phonon-assisted  processes,  one  reason  for  the 
significant  difference  in  energy  transfer  characteristics  for 
different  impurity  ions  may  be  due  to  the  much  higher  con¬ 
centration  of  Nd'  ’  ions  in  the  YVO,:Nd' ‘  samples.  The 
model  used  to  explatruhese  results  is  similar  to  that  used  for 
interpreting  the  results  of  similar  investigation  of  Nd-doped 
garnet  crystals." "  The  main  difference  in  the  models  used 
for  these  two  sy  stems  is  that  in  the  garnet  hosts  the  interme¬ 
diate  state  lor  the  phonon-assisted  energy  transfer  is  ihe  first 
level  of  the  ground  state  manifold  whereas  for  the  yttnum 
vanadate  host  it  is  the  upper  component  ofthe  nietastable 
stale.  The  reason  for  this  may  be  that  the  splitting  of4/-'  . 
level  is  much  smaller  in  YVO,  crystals  than  tt  is  in  garnets. 

It  should  be  emphasized  that  the  agreement  between 


theory  and  experiment  at  high  temperatures,  shown  in  Fig. 

1 1 ,  is  quite  good  considering  that  there  are  no  adjustable 
parameters  in  the  theory.  The  poorer  agreement  at  low  tem¬ 
peratures  indicates  that  other  energy  transfer  processes  are 
becoming  important  in  this  temperature  range. 

The  origin  of  the  different  types  of  crystal  field  sites  for 
the  Nd3  ’  ions  is  not  known  but  such  site  multiplicity  has 
been  observed  in  this  system  by  previous  investigators. 16  It  is 
well  known  that  a  major  problem  with  yttrium  orthovana¬ 
date  crystals  is  the  presence  of  imperfections  introduced  dur¬ 
ing  crystal  growth.  Chemical  or  structural  defects  located 
near  to  Nd3  "  ions  could  result  in  nonequivalent  sites.  The 
sites  which  are  selectively  excited  by  the  laser  are  the  domi¬ 
nant  types  of  sites  for  Nd3  +  ions  in  the  activator  type  of  sites. 

The  value  of  D  found  here  for  energy  migration  among 
Nd1  '  ions  in  YV04  is  within  the  range  of  values  between 
10  3  and  10  13  found  for  other  rare-earth  ions  in  solids. 
The  diffusion  coefficient  for  energy  migration  among  the 
Nd3  ’  ions  is  given  by 


Z>=  3.4C4 


■c'/3ri  (j\  w:  -  2(J-)' 

\  •  (j  £j-  +  &r2J 

:r/*J£j2][l  +  exp(d£„A/)]exp(  —  S/kT) 

(18) 


and  is  found  to  be  1.3  V  10  -  IJ  cm 2  sec  at  295  K.  This  is 

somewhat  smaller  in  the  vanadate  host  than  in  the  garnet 
host  due  to  the  increased  value  inhomogeneous  broadening 
of  the  transitions. 

An  attempt  was  also  made  to  fit  the  data  using  a  Monte 
Carlo  computer  simulation  procedure17  to  account  for  the 
ihe  fact  that  the  Nd3  '  “exciton"  is  hopping  on  a  lattice  of 
randomly  distribution  sites.  A  good  fit  to  the  data  is  shown  in 
Fig.  9  for  the  100  K  data  assuming  a  critical  interacation 


distance  of  Ru  =  23  A  which  is  similar  to  the  results  obtained 
for  Nd3  v  ions  in  gamet  hosts. 

In  conclusions,  laser  time-resolved  spectroscopy  tech¬ 
niques  have  been  used  to  determine  the  characteristics  of 
both  host-sensitized  energy  transfer  and  energy  migration 
among  activator  ions  in  YNO<:Nd3  *  crystals.  Understand¬ 
ing  these  processes  is  of  fundamental  importance  in  the  use 
of  this  material  in  laser  and  phosphor  applications. 
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A  pulsed,  tunable  dye  laser  was  used  to  selectively  excite  Nd  '  ions  in  nunequiv  jler.  .- 
ervstal  field  sites  in  "NdAijtBO,  i4  crystals  and  energy  transfer  between  ions  in  different 
types  of  sites  was  studied  bv  monitoring  ihe  tune  evolution  of  ihe  fluorescence  spectrum 
The  results  show  that  the  energy  transfer  rate  vanes  as  r  1  -  and  increases  with  tempera¬ 
ture.  The  predictions  of  vanous  models  of  phonon- assisted  energs  transfer  are  compared  to 
the  results. 


1.  Introduction 

Neodymium  aluminum  borate,  NdAljfBOj).,,  is  one  of  the  several 
"stoichiometric  laser  materials”  which  have  gained  recent  interest  due  to  their 
use  in  mimlaser  applications  (1,2).  In  these  materials  a  high  concentration  of 
Nd3*  ions  can  be  present  (5.43  x  10:t  em~3  for  NdAlj(BOj)4j  without  strong 
concentration  quenching  of  the  fluorescence  emission.  Since  these  materials  are 
relatively  new,  there  are  still  many  aspects  of  their  physical  properties  that  are 
not  well  characterized  and  understood. 

NdAljfBOj  )i  (NAB)  has  a  complex  crystal  structure  with  marked  deviation 
from  inversion  symmetry  and  significant  isolation  of  the  Nd  polyhedra  [3.4], 
The  strong  odd-parity  crystal  field  admixture  of  the  f  and  d  configurations 
results  in  a  short  radiative  lifetime  and  the  highest  cross  section  of  any  Nd 
material.  Although  laser  action  has  been  observed  in  NAB.  the  spectroscopic 
properties  of  the  material  have  not  been  well  characterized  [5,6],  One  property 
of  special  interest  which  has  not  been  investigated  is  the  interaction  between 
Nd3*  ions  which  may  lead  to  fluorescence  quenching. 

We  report  here  the  results  of  a  study  of  energy  transfer  among  Nd3*  ions  in 
NAB  crystals.  Time-resolved  site-selection  spectroscopy  techniques  were  em¬ 
ployed  to  determine  the  magnitude,  time  dependence,  and  temperature  depen¬ 
dence  of  the  rate  of  energy  transfer  between  Nd’*  ions  in  nonequivalent  types 
of  crystal  field  sues.  Several  possible  interpretations  of  the  results  are  dis¬ 
cussed. 
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2.  Experimental  techniques 

The  samples  consisted  of  small  single  crystal  pieces  of  NdAI,(BO,)4.  These 
were  mounted  on  the  cold  finger  of  a  cryogenic  refrigerator  capable  of  varying 
the  temperature  from  about  10  K.  up  to  room  temperature.  A  block  diagram  of 
the  experimental  apparatus  is  shown  in  fig.  I.  A  nitrogen  laser-pumped  tunable 
dye  laser  was  used  to  excite  the  sample.  Using  rhodamine  6G,  a  pulse  of  less 
than  10  ns  in  duration  and  about  0.4  A  halfwidth  was  obtained.  The  sample 
fluorescence  was  focused  onto  the  entrance  slit  of  a  one-meter  monchromator 
and  delected  bv  a  cooled  RCA  C3I034  photomultiplier  tube.  The  signal  was 
averaged  by  a  boxcar  integrator  triggered  by  the  laser  so  that  either  the 
lifetimes  or  the  fluorescence  spectra  at  different  times  after  the  excitation  pulse 
could  be  recorded. 

The  wavelength  of  the  dve  laser  excitation  light  was  approximately  5920  a 
which  is  absorbed  in  one  of  the  higher  excited  states  of  the  Nd'1 '  ions.  After 
relaxation  of  the  4Fj  ,  metastable  state,  the  emission  to  the  various  compo¬ 
nents  of  the  4lt)/2  ground  state  manifold  was  monitored.  The  lifetime  of 
fluorescence  was  found  to  be  20  ns  and  independent  of  temperature  between 
10  and  295  K. 

The  fluorescence  spectrum  of  NdAI  at  1 1  K  in  the  region  from  8780 

to  8940  A  is  shown  in  fig.  2  for  two  different  excitation  wavelengths.  The  lines 
from  each  transition  exhibit  a  significant  amount  of  structure  and  the  domi¬ 
nant  structural  components  vary  significantly  with  excitation  wavelength.  The 
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Fig  2.  Fluorescence  spectra  of  NdAl  ,(1)0,  )4  ji  !  1  K.  for  two  different  excitation  wavelengths 

relative  intensities  of  these  peaks  vary  in  a  uniform,  continuous  wav  between 
the  two  extremes  shown  in  fig.  2  as  the  laser  wavelength  is  scanned  from  5919 
to  5925  A.  This  structure  can  be  attributed  to  the  presence  of  ions  in  slightly 
different  types  of  crystal  field  sues.  By  tuning  the  dve  laser  wavelength,  ions  in 
specific  types  of  sites  are  selectively  excited.  Energy  transfer  between  ions  in 
spectrally  different  sites  can  be  studied  by  monitoring  the  time  evolution  of  the 
fluorescence  spectra  after  the  laser  pulse.  Fig.  3  shows  the  spectrum  at  two 
times  after  the  excitation  pulse  for  a  specific  excitation  wavelength  at  low 
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i ;!£.  .*  fluorescence  spectra  of  NdAl  5(B()i  i4  ai  1 1  K.  ji  iv*o  times  jficr  ihe  cxuution  pulse 
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Fi^.  4  Ratios  of  the  integrated  fluorescence  intensities  of  activator  and  -cnsitizer  transitions  as  a 
function  of  tune  after  the  excitation  pulse  at  1 1  K.  and  4U  K  See  text  for  explanation  of  theoretical 
lines 


temperature.  At  short  times,  the  high  energy  components  of  the  transitions  are 
dominant;  as  time  increases  the  low  energy  components  grow  in  intensity  with 
respect  to  the  high  energy  lines. 

The  quantitative  analysis  of  energy  transfer  from  Nd3'  ions  in  sites 
producing  the  high  energy  transitions  (sensitizers)  to  ions  in  sues  resulting  in 
the  low  energy  transitions  (activators)  can  be  obtained  by  plotting  the  time- 
dependence  of  the  ratios  of  the  integrated  fluorescence  intensities  of  the 
activator  and  sensitizer  transitions.  The  results  are  shown  in  fig.  4  for  1 1  and 
40  K.  At  higher  temperatures  the  lines  are  broadened  to  the  extent  that  the 
transitions  from  ions  in  nonequivalent  sites  cannot  be  esolved  accurately.  The 
solid  lines  in  the  figure  represent  theoretical  fits  to  the  data  discussed  in  the 
following  section. 


3.  Interpretation 

The  time-resolved  spectroscopy  results  can  be  interpreted  using  the  phe¬ 
nomenological  rale  parameter  mode!  shown  in  fig.  5.  The  "sensitizer"  ions  are 
those  in  the  site  preferentially  excited  by  the  laser  at  a  rate  IT,  whereas  the 
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"activators'1  ire  ions  in  Mies  which  receive  the  energy  through  energy  transfer 
as  vvell  as  being  excited  by  the  laser  at  a  rate  IT.  n  and  /tj  are  the 
concentrations  of  ions  in  the  excited  states,  is  the  rate  of  energy  transfer, 
and  3  is  the  fluorescence  decay  rate  which  was  found  to  be  approximately  the 
same  for  ions  :n  both  types  of  sites. 

The  rate  equations  describing  the  time  evolution  of  the  populations  of 
excited  states  are 


0  n  k  d !  =  \V\  -  Bn  (  -  iovJ  n  v .  ( 1  I 

Jrtj  df  =  ll\  ~  B>i  j  -r  u.\an ( 2 ) 


These  equations  can  be  solved  assuming  a  delta-function  excitation  pulse  and 
an  explicit  time  dependence  for  the  energy  transfer  rate.  A  v  ariety  of  different 
theoretical  models  were  tested,  but  it  was  found  ihat  the  best  fits  to  the  data  at 
both  low  and  high  temperatures  were  obtained  with  an  energy  transfer  rate 
which  varied  as  /  :  *.  The  solutions  of  eqs.  ( 1 1  and  (2)  are 

n  ( -  n  <0l  e  a‘ ' ) 

/MM  -«,(0)[e  -  e  tu  -  -•/,('))  e  (4, 

where  the  time  dependence  of  the  energy  and  transfer  rate  is  written  explicitly 
as  w.  =u :i  '  :.  The  ratio  of  integrated  fluorescence  intensities  is  then  given 
by 


'M01 


(5' 


where  K  is  a  constant  containing  factors  including  the  radiative  emission  rates 
The  solid  lines  in  fig.  4  represent  the  predictions  of  this  expression  treating  K. 
u  and  rtj 0)/;js(0)  as  adjustable  parameters. 

The  temperature  dependence  of  the  energy  transfer  rate  was  determined  bv 
measuring  the  fluorescence  intensities  at  different  times  after  the  excitation 
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pulse  .is  ,i  function  of  temperature  and  fitting  die  data  using  eq.  (5i.  1  he  re-  a.':- 
ire  shown  in  fig.  5.  I  he  transfer  rate  was  found  to  he  approximate!'.  ->  lotum 
helow  25  K.  and  to  increase  at  higher  temperatures. 


4.  Discussion  and  conclusions 

i  he  tune  dependence  of  the  vnvrus  tr.mstei  i.ce  ,  ,  cn-itr.e  nea  it-  r  : 
the  mechanism  of  clients  transier.  I  lie  :  san.it;. -n  .  r-ser-c  :  ::er-  s  in 

due  to  several  tspes  >>i  processes,  '.ngie-step  cievtKe  dipole  m:  -ic  nter.n'.i.  a 
between  random!}  distributed  -eri'itr/.ers  and  act r s a t.  r - ;  muitotep  ::ie:e. 
migration  on  a  one-Jtmensional  lattice:  .md  :;ap-niodaia:ed  .-tier..-.  -marati...-: 
:n  three  dimensions  Pile  first  if  these  possibilities  ,s  genera:!-,  the 
common  and  can  he  serified  tluougn  theoretic  a!  piectctioiis  a  the  :iu.'m,|..ue 
of  the  energy  transfer  .ate.  At  I  I  K  the  measured  salue  of  w-  .s  ■>.'.<  us 
The  theoretical  expression  for  the  r.ite  of  simtie-'iep  electric  dipo.e  -hp-'ie 
energs  transfer  among  randomls  distributed  ions  is  " 

w here  c'4  is  the  concemratton  of  xtivutois.  -  n  tlte  minnstc  fan  rc  -cence  deca-. 
time  of  the  sensiti/ers  and  the  critical  interaction  d. stance  o  e:\en  bv 


where  o,  is  the  quantum  etficienc;.  -  i  the  en.stti/.cr.  •  .s  t,c  se:'...uor  arenctit 
of  the  activator.  i> j  is  the  average  'vase  number  m  die  r-titon  .f  spear..) 
overlap,  and  '.l  the  spectral  oserlap  integral.  If  sing  tlte  measured  value  t  tt.c 
energy  transfer  rate  m  eq.  |h|  along  with  the  v allies  of  (  \  -  pa  •  ] tj-'  cm 
determined  by  comparing  relative  -.peetra!  intensities  at  different  excitation 
wavelengths  and  r  -  5"  u>  which  is  the  value  measured  m  iightlv  doped  norate 
samples  jdj  gises  a  measured  value  of  H  ,  -  5.0  A.  1  his  can  he  compared  to  the 
theoretical  salue  predicted  by  eq.  <?).  where  /,  is  calculated  from  .-peetrui 
measurements  to  be  2..'  -  10  *  and  the  overlap  integral  .s  approximated  bv 
the  overlap  of  iwo  l.oremt-ian  lines 


2.1-  -  2o"  )'  '  I  v‘  -  r\’  ' 


At  1  I  K  the  predicted  vaiue  is  R,  -  7.6  A.  I  Ire  difference  m  live  predicted  and 
measured  values  of  R  ,  is  significant.  Although  some  of  'lie  values  of  the 
parameters  used  m  obtaining  these  results  are  rough  estimates,  -uch  as  o.  and 
12.  any  reasonable  choice  ol  values  given  at  least  a  factor  of  2  difference 
between  the  two  determinations  of 


12  o 


1  /  j 
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/J  !\  c.  /'  w.ti  ‘-.'v..".-  .  ;n  '  \  l4  I,  ,  lit  ■  •  .  r  ii , .  ; ' 

i  ho  inconsistency  described  above  nisiujics  that  ihc  single-step  mechanism 
in  iK. i  responsible  lor  the  observed  encrgv  transfer  in  this  s'. stem.  Hie  econo 
t\pe  of  possible  energy  transfer  process  bavinc  the  observed  time  dependence 
:s  .i  random  '.sails  in  one  dimension.  In  this  case  the  Nd'  '  u  n  m  the  essited 
motjslable  stale  is  treated  as  a  Frenkel  evsitoii  which  .an  migrate  to  'Hkt 
-•ensiti.-er  ions  before  becoming  trapped  .it  an  .utivjtoi  site.  Site  'i.cicc. 
■.rar.ster  rate  for  a  random  '.talk  process  c.m  re  e\p:c.M.'d  as  ::.e  :ate  at  .slttcti 
'tie  ."...ton  sample  new  -ate  :miit. plied  t.ve  frusta  n  f  ‘ot...  Hies  .  e.  aptesi 

o.  itivatois  , o  :|  l  -mg  the  expression  for  the  -umplma  functii  n  ’.r  a 
iniiensioiiai  lattice  give  ;  I  !  j- 

~  -  l2  '  ;  V  |  -  «  !<,']'■  'i 

si  ere  :t  has  been  assumed  that  each  step  in  the  random  w  .A  .ikes  place  ns  an 
-lecttte  shrole  dipole  interaction  median. sin  between  ions  -eparated  by  a 
.int.ince  H.  I  he  closest  Nd  ion  separation  in  NdAI  ,i  B();  cr.  st.Js  is  r.J!  V 
I  ..ng  this  and  the  measured  value  of  the  energv  transfer  rate  at  ea  <'■)  etses  a 
.  dae  ot  H  "  '.d  A  which  i>  consistent  with  the  the-  •rettcai  predtctions  of  ou. 

I  lit*  third  possibility  that  the  ihser-.eii  energy  transfer  ctt-trasicrislics  are 
.tee  o  trap-'iioduiateO  energy  migration  :n  three  ditnensums  can  not  be  easih. 
proven  or  dtsproven  since  any  type  of  irap  distribution  necessary  to  lit  the  data 
can  be  postulated  [Id],  An  important  question  is  whether  or  not  exciton  motion 
restricted  to  .me  dimension  is  feasible  m  this  type  -f  ervstai.  Flic  published 
crystal  structure  of  NJAl-tliO,  j4  indicates  that  a  Nd  ion  site  has  ti  nearest 
neighbor  Nd  tons  in  various  crystallographic  directions  'o'.  However,  more 
recent  structural  measurements  describe  this  crystal  as  having  a  two- 
dimensional  layered  structure  |lc).  I  he  feasibiutx  of  a  preferred  migration 
atrcstton  must  await  publication  (  the  details  of  this  structure.  It  should  be 
emphu.'i.’.esi  that  theoretical  models  lor  fitting  the  I  RS  data  .ire  •. cry  sensitive 
to  the  'tine  siependence  ot  the  proposed  energy  transfer  mechanism  and 
.itts'inpis  to  fit  the  data  with  two-dimensional  or  three-dimensional  random 
talk  models  or  with  single- step  higher  order  muitipole  processes  were  ail 
unsuccessful. 

The  observed  temperature  dependence  of  the  energy  Tansfcr  rate  .an 
provide  further  information  about  the  nature  of  the  energy  transfer  mechanism 
arts!  the  role  placed  by  phonons  at  high  temperatures.  Several  different 
processes  can  contribute  to  the  increase  ;n  the  transfer  rate  with  increasing 
emperature.  These  include  the  increase  in  the  spectral  overlap  integral,  the 
increase  in  the  number  of  phonons  available  to  make  up  the  mismatch  between 
sensitizer  and  activator  transitions  ( A £sa  1 2  cm  !,  and  the  increase  in 
population  of  higher  energy  sensitizer  and  activator  levels  which  can  take  part 
in  the  transfer  process  [14j.  For  this  case  the  first  two  types  of  processes  were 
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found  to  be  inconsistent  with  the  data  siiown  in  fig.  ft.  The  third  mechanism 
van  be  applied  to  either  .-angle-step  or  ntultistep  transfer.  F-<r  a  single-step 
process  the  expression  for  the  energy  transfer  rate  in  eq.  (6.  would  be  written 
as  the  sum  of  transfer  rates  for  each  set  of  transitions  weighted  by  the 
populations  of  the  initial  states  of  the  transitions.  The  observed  temperature 
dependence  ts  not  consistent  with  populating  ’he  first  excited  level  of  the 
ground  state  manifold  of  the  activator.  £  ,  *s  bO  cn.  ‘ !  A  qualitative  fit  to  the 
data  shown  ;n  fig.  ft  can  be  obtained  assuming  the  important  activation  energy 
is  A  a  ib  ~  25  cm  “ 1  which  is  needed  to  populate  the  higher  crystal  field 
component  of  the  metaslable  state  of  the  sensitizer.  However,  the  room 
temperature  spectrum  shows  that  the  b-  I  transition  is  less  intense  than  the 
u  -  I  transition  and  thus,  energy  transfer  stiould  not  increase  as  the  h  level 
becomes  populated  from  the  a  level  (91. 

The  most  consistent  fit  to  the  data  is  found  from  the  one-dimensionai 
random  walk  theory  described  by  expressing  eq.  (9>  j>  a  sum  over  the 
thermally  activated  transitions  involved  in  the  migration  If  the  -urn  is  am¬ 
plified  to  include  only  the  terms  for  the  ground  -tale  and  tne  first  excited 
component  the  energy  transfer  rate  is 
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where  rM.  .inJ  .ire  the  hupping  times  involving  the  two  tvpes  of  tr.in.sition •>. 
File  good  fit  •.*>  the  data  shown  hv  die  solid  line  in  fig.  6  was  found  with 
A t'._  =  f’(l  cm  .  =s  o.trf  is  •  and  S:  6.7?  us  .  Both  the  v.iiue  M 

the  activation  energy  and  the  fact  that  are  consistent  with  spectral 

oriservations.  The  magnitude  of  i,{  is  consistent  with  the  theoretical  predict¬ 
ion  =  r  I  H.,  R 

In  summers,  tne  energy  transfer  ehat  ictenstics  of  NA1J  crystals  are  con¬ 
sistent  with  a  one-dimensional  random  walk  in  which  the  hopping  rate 
increases  as  the  first  evened  component  of  the  ground  state  multiple’,  becomes 
thermally  activated.  The  diffusion  coefficient  at  low  temperatures  is  I) 
-  R~  rM,  =  1.3  •  id  ’  ear  sec  These  characteristics  are  quite  different 
from  those  found  :n  other  Nd-  '  -doped  crystals  ( 1  5 j .  The  differences  may  he 
associated  with  the  different  crystal  structure  of  NAB  which  til;  must  ->e 
clarified.  It  should  be  pointed  out  that  the  conventional  approach  to  the  stuitv 
of  energy  transfer  lias  been  used  in  this  investigation  in  that  a  phenomenologi¬ 
cal  rate  parameter  model  was  utilized  to  obtain  the  primary  characteristics  of 
the  transfer  rate  from  the  experimental  data  :ad  seeondarv  information  was 
obtained  from  comparison  with  the  predictions  of  -tandard  theoric-  Mote 
sophisticated  theories  are  now  being  developed  which  raise  question'  about 
some  aspects  of  this  conventional  approach  [16]  but  they  are  not  ;  e;  estah- 
iisned  to  the  pomi  where  they  ear.  be  used  to  analyze  new  data  itch  as  that 
obtained  in  this  studs. 


Acknowledgements 

This  research  was  supported  by  the  LS  Arms  Research  Office.  The  audit  rs 
are  grateful  to  Professor  O.  Blasse  for  stipplving  the  samples.  .mu  apnreeiate 
the  helpful  discussions  with  Dr.  W.  Zwicker  concerning  the  crvdnl  structure  of 

NAB. 


References 


;lj  S  R  Chinn.  HA  P  Hong  and  J  A  Pierce.  Laser  Focus  id  i  !<5Au  (4 
’d]  H  (j  Damelmcvcr,  Festkorperproblcme  (Advances  n  Sot.d  State  Phvxs'  .  '  ti  1 

Ouetsscr  ' Pcreamon  Wwei,  Braunschweig,  .  s  "  3  p  d53 
31  HY-P  Hone  and  K  Dwight.  Mat  Res  Bull.  1  .  1  ma,  'sd 

.4)  O  iarchow.  F  Lmz  .eiJ  K  1!  Kiaska.  Dokussionslacune  dcr  ArKmsceuKnoou.a  k::s:...ioL; 
rafinc  t'J  1 1  rdl  '  hZ. 


91 


J 


^  -  a  .Vu/Jo/\  .V  t  /*"*»*•</'  S;>t'(  •<,'*»  i/Ut/;n  ./ff*  -j.M-.'un 

V  1  liilxk.  1  l  K^jraies,  N  1  l.ci>n\uk.  V  ,\  PaMikov  A  V  Pa.shk' ’•  j.  1  i  l.nivhc 
A  V  Shcsukov.  So*  PhvN  Dok.1  23  ( 2in 
A;  ! i  I >  Haitendorlf.  (j  I lubcr  and  H  <. i  DaniclnwaT  J  Ph.N  C‘||  i19~m23*^ 

■’]  M  liuAuti  and  P  Hira\ama.  J  Chem  Ph%>  15  <  P^-Oi  .  v7j< 

Sj  D  L.  Dcxicr.  J  (‘hem  Phv>  21  (1953)  a3(» 

w!  SR  Chinn  and  HYP  H<>n»;.  Opuo  Comm  !5i!'C5»345 

Oj  /.  (.j  S*h>n  and  R  C  Powell.  Phvs  Ro  lit>  i  1^"'2j  4»'35 

t.W  Moniroil.  P:  v  S-.mp  Appl  Maih  Am  Math  Soc  5  i !  *>*>3 1  i-0 


•  -  ; 

R  C 

P.MACiJ 

and  - 

'  i  / 

Soon.  Ph\>  Rc\  !i- 

•  .  ;‘^2)  1 54" 

i  ?: 

il  C 

i'ametmeNei 

r.  prt 

j'.  the  intc;: 

tat.onjj  l.a.NCf  SstnpOMUl 

m  ; 

orvr. 

Ph\  > 

R,- 

.  \U  . 

‘•“h 

J  N1 

ii;err 

.  ..no 

.1 

i\  V* 

I'..  Phv.  Rr. 

ii !  j  [tj-o»  '*2. 

1.  i.) 

'•U"kie 

u»10 

R  v 

Pi  -VtOli 

.  .  *ro  RC\ 

!i2< 1  t  "5. 

M  / 

.I'fca!,  R 

('  W 

*vtr;: 

:  u  ! 

l.nmj.v.1;  .nn 

1  H  ! ):ii.:rtok>.  J  \p;  i  .■ 

C  K. 

Sardar 

and 

P  as,-:i 

>  \ppi  PS: 

a  >  '  ;  : ,  ;:s2'> 

U>j 

i) 

Si  j  her. 

Ph\  ■ 

R; 

n  !J2n  • 

2}ir. 

K.  N 

y  »h«  •>:! 

anu 

is  i. 

iiu-'er. 

,  j  !  .miiii  2 

.  i  :«MJ|  225 

K.  K. 

*  n, 

aiar 

:\  and 

:  )  L.  Huiv;. 

Phw  Ke\  Ml  :  ‘Mil  2 

IV.  3 


ENERGY  TRANSFER  AMONG  NEODYMIUM  IONS  IN  GLASS  HOSTS 

Larry  L).  .'VrlN  ,  Richard  0.  ?r  well  and  Edwin  E .  rr<a-u 
Physics  Department,  Ok  1 auoma  State  L'niv.  *  ,  Stillwater,  I.  ~  j  7  A  , 

and 

Marvin  J.  Leber 
Lawrence  Livermore  National  Laboratories , 

Energv  transfer  between  NdJ+  ions  in  silicate,  y ema¬ 
nate,  phosphate,  and  pentaohosniiate  glasses  has 
Leer,  studied  between  12  and  300  K.  The  results  at 
12  K  acree  with  relative  values  predicted  : or  electric 
i  :>e  le-d  ipole  energy  transfer.  Yerv  di;  fer-.-nt  :  c.rera 
t  ire-  dependences  are  observed  for  the  energy  transfer 
in  these  samples. 


Tunable  dye  lasers  have  proven  to  be  an  important  tool 
Properties  of  ions  in  disordered  host  materials  su-.h  as  -La 
me  results  -  :  us  ins:  laser  excited  tine-resolved  ••• :  te-selec 
imiques  t ..  study  the  energy  t  roaster  between  Nd  ions  i 
sites  in  four  different  class  hosts . 

The  compositions  of  the  glasses  investigated  are  liste 
experimental  apparatus  used  in  this  work  has  been  described  rrevieusl 
Rhodaraine  bG  was  used  in  the  dve  laser  to  obtain  an  excitation  pulse 
bandwidth.  This  excited  the  Ley  -  -G--y  "Gj  y  transition  and  the 
emission  from  the  “F  :  ■>  -  - 10  .  .  transition  was  observed. 


Table  1:  Summary 
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Figure  1  shows  examples  of  the  fluorescence.  For  all  of  the  samples  the 
selectively  excited  spectra  evolves  with  time  toward  the  broadband  excited 
spectra.  The  fluorescence  decay  times  vary  across  the  emission  band  and  the 
amount  of  variation  is  different  for  each  sample  and  becomes  less  at  high 
temperatures.  The  lifetimes  measured  at  the  peaks  of  the  emission  bands  at 
12  K  are  listed  in  Table  1. 

The  results  described  above  show  that  laser  excitation  selectivity  excites 
Nd3+  ions  in  different  types  of  local  site  environments  and  that  energy  transfer 
takes  place  between  ions  in  these  different  types  of  local  sites.  The  function 
describing  the  energy  transfer  can  be  written  as  [2]. 

a(t)  =  /*Ui2  [l(w,t)-I(w,°°)  W2  [l(w,o)-I(u),°°)]d<jj. 

J  ^1  U>1 
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Figure  1 

Normalized  emission  spectra  at  12  K  -  broad  band  excitation;  ....  laser 

excitation  at  short  times  (10  us  for  ED2,  LG,  and  LP  and  1  us  for  NP) ;  - 

laser  excitation  at  long  times  (400  us  for  ED2,  600  us  for  LG,  300  us  for  LP  and 
50  us  for  NP) .  94 


a(t)  can  be  approximated  by  an  exponential  curve  with  the  e  time  given 
by  the  values  of  a  listed  in  Table  1. 

3+ 

The  rate  of  energy  transfer  between  two  Nd  ions  separated  by  a  distance 
R  can  theoretically  be  estimated  by: 


(2Ja  +1)  (2Jfi  +1) 


(2/3)  (2n/fi)  e*/R6 


L  £  Q 

t»2,4,6 


x  £  nBt  |<  jb|  |u  (t)|  I  jb>|2. 
t=2,4,6 


This  assumes  a  resonant  electric  dipole-dipole  interaction.  L  is  the  integrated 
spectral  overlap  and  the  line  strengths  are  expressed  in  terms  of  the  Judd-Ofelt 
parameters  and  the  reduced  matrix  elements  of  the  unit  tensor  operator.  The 
matrix  elements  and  Judd-Ofelt  parameters  were  calculated  from  spectral  data  and 
the  results  were  used  in  Eq.  (2)  along  with  estimates  of  average  values  of  R  ob¬ 
tained  from  the  known  concentrations  assuming  uniform  distributions  of  Nd  ions. 
The  spectral  overlap  integral  in  this  case  involves  all  of  the  transitions 
between  the  various  Stark  levels  of  the  excited  state  and  ground  state  m  nifolds. 
This  is  proportional  to  the  extent  of  splitting  of  the  J-manifolds  which  can  be 
approximated  by  the  effective  linewidth  6Xeff  3  and  these  are  listed  in  Table 
1.  Note  that  6Aeff  is  actually  measured  for  the  ^3/2  ^*11/2  transition  but 

this  is  proportional  to  the  effective  width  of  the  4v  ^  4t  transition  3. 

*3/2  *4/2 


Figure  2 

Measured  energy  transfer  and  quenching  rates  at  12  K  plotted  versus  theoretically 
predicted  electric  dipole-dipole  energy  transfer  rates. 


As  a  rough  approximation  we  assume  the  variation  in  6v  from  sample  to  sample  is 
proportional  to  the  effective  line  widths  (3).  In  Fig.  2  the  measured  values  of 
at  12  K  are  plotted  versus  the  predictions  for  P  expressed  in  arbitrary  units. 
a**3  varies  uniformly  with  P  from  sample  to  sample  which  shows  that  the  glass  com¬ 
position  effects  the  strength  of  the  energy  transfer  through  L  and  Q. 


Another  parameter  of  interest  is  tne  radiationless  quenching  rate  Q  = 

-  Tra(j-1  where  Tfi  is  the  fluorescence  decay  time  and  Tra{j  is  the  radiative 
decay  time.  The  values  of  Q-^-  at  12  K  are  plotted  versus  P  in  Fig.  2.  The  vari¬ 
ation  of  Q-3  with  P  is  the  same  as  that  of  which  indicates  that  the  physical 
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processes  causing  the  differences  in  from  sample  to  sample  have  similar 
effects  on  Q"l. 

The  temperature  dependences  of  a  are  quite  different  as  shown  in  Fig.  3. 
Since  there  are  several  different  physical  processes  which  can  contribute  to  the 
temperature  dependence  of  the  energy  transfer  rate  it  is  difficult  to  give  a  def¬ 
initive  interpretation  to  the  observed  results.  One  possible  interpretation 
which  is  also  consistent  with  the  different  spectral  characteristics  seen  in  Fig. 
1  is  the  following.  Energy  transfer  in  the  NP  and  LP  glasses  may  be  associated 
with  single-phonon  assisted  process  involving  stimulated  phonon  emission  whereas 
the  transfer  mechanism  in  the  LG  glass  may  be  a  single-phonon  assisted  process 
involving  spontaneous  phonon  emission.  The  behavior  of  the  ED2  glass  can  be  ex¬ 
plained  either  by  two-phonon  assisted  processes  or  enhanced  absorption  transition 
strengths  due  to  thermally  populating  higher  Stark  components  of  the  ground  state 
manifold.  Further  studies  are  currently  underway  to  correlate  energy  transfer 
properties  with  properties  of  the  glass  hosts  with  different  compositions. 


Figure  3 

Temperature  dependences  of  average  ion-ion  energy  transfer  times. 
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V.  MISCELLANEOUS  STUDIES 

The  manuscript  in  Section  V.l  describes  the  results  obtained  using  photo¬ 
acoustic  spectroscopy  techniques  to  determine  the  quantum  efficiencies  of 
several  different  Nd°^  doped  samples.  In  Section  V.2  results  are  presented 
on  the  spectroscopic  property  of  ruby  crystals  under  hydrostatic  pressure. 
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Radiationless  decay  processes  of  Nd 3  +  ions  in  solids 
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Laser  photoacoustic  spectroscopy  measurements  were  made  on  Nd  * '  ions  in  garnet,  vanadate,  and 
pentaphosphate  host  crystals.  The  variations  of  signal  intensities  with  chopping  frequency  of  the 
incident  light  are  not  in  agreement  with  the  predictions  of  standard  photoacoustic  signal  generation 
theory.  The  results  are  distinctly  different  for  concentrated  and  dilute  Nd-doped  crystals,  indicating 
that  the  mechanisms  for  generating  heat  have  different  characteristics  in  these  two  types  of  samples. 

The  determination  of  radiative  quantum  efficiencies  of  these  materials  by  photoacoustic  spectroscopy 
techniques  is  also  described.  These  results  are  compared  with  those  obtained  by  other  methods. 


INTRODUCTION 

Photoacoustic  spectroscopy  (PAS)  techniques  have  recently 
received  a  great  deal  of  interest  as  a  method  for  characterizing 
radiationless  relaxation  processes  of  ions  in  solids.1"’  We 
describe  here  the  results  of  PAS  investigations  of  Nd:,+  ions 
in  several  different  types  of  host  crystals.  The  data  is  inter¬ 
preted  in  terms  of  radiationless  decay,  concentration 
quenching,  and  energy  migration  processes  and  a  method  for 
determining  the  quantum  efficiencies  of  the  samples  is  de¬ 
scribed.  Problems  with  understanding  the  signal  generation 
process  with  laser  excitation,  and  problems  with  the  accuracy 
of  quantitative  measurements  are  discussed. 

Our  experimental  setup  for  PAS  measurements  has  been 
described  previously.1'4  The  only  difference  for  the  work 
described  here  is  that  the  individual  lines  of  an  argon  ion  laser 
were  used  as  the  excitation  source.  The  laser  power  was 
continuously  monitored  and  stabilized  at  a  level  of  0.15  W. 
The  exciting  light  was  chopped  at  frequencies  varying  from 
110  to  '2700  Hz  and  focused  onto  the  samples  which  were 
placed  on  the  quartz  exit  window  of  the  PA  cell.  The  micro¬ 
phone  was  mounted  at  90°  to  the  exciting  light  and  behind  a 
baffle  to  prevent  scattered  laser  light  from  reaching  it.  The 
signal  was  sent  through  a  preamplifier  to  a  lock-in  amplifier 
and  read  out  on  a  digital  voltmeter  after  adjusting  the  lock-in 
phase  lor  signal  maximum.  The  signal  to  noise  ratio  was 
(tetter  than  100  to  1  in  all  cases  and  the  measured  background 
cell  signal  was  at  least  200  times  smaller  than  the  sample  PA 
signal. 

The  samples  investigated  were  all  single  crystals  between 
1  and  2  mm  thick  and  between  5  and  10  mm  on  a  side  with 
polished  faces.  Three  types  of  samples  were  studied:  the 
garnets  Y.fAI50,g  (0.85%  Nd),  Y.)Ga50,2  (0.25%  Nd), 
Y:i(Alo.sGai>s).',Ois  (0.85%  Nd).  and  Nd.iGa^Oiz,  the  vanadate 
YV04  (2.0%  and  3.0%  Nd);  and  the  pentaphosphates 
Nd,  Yi_,  P'.On  with  x  ranging  from  0.1  to  1.0. 

Two  types  of  experiments  were  performed;  the  measure¬ 
ment  of  PA  signal  intensity  versus  chopping  frequency  and 
the  determination  of  quantum  efficiencies. 

I.  PA  SIGNAL  VERSUS  CHOPPING  FREQUENCY 

The  PA  signals  at  maximum  phase  were  recorded  at  12  or 
more  chopping  frequencies  cc  after  excitation  with  the  5145-A 
line  of  the  argon  laser.  Examples  of  the  results  are  shown  in 
Figs.  1-3.  For  both  samples  of  YVOr-Nd3*,  the  PA  signal 
varies  as  e.71  throughout  the  entire  range  of  frequencies.  For 


alt  three  of  the  lightly  doped  garnet  samples,  a  r~'  dependence 
was  also  observed  as  shown  for  Y:tAI>)0|2:Nd;,+  in  Fig  2. 
However,  for  the  concentrated  neodymium  gallium  garnet 
sample  the  PA  signal  was  observed  to  vary  as  vi  'n  throughout 
the  entire  frequency  range.  For  the  pentaphosphate  crystals 
the  10%  Nd  sample  exhibits  a  PA  signal  intensity  which  varies 
as  v~'  over  the  entire  frequency  range  while  the  PA  signal  for 
the  100%  Nd  sample  varies  as  ef1  up  to  about  450  Hz  and  then 
varies  as  at  higher  frequencies.  Samples  with  interme¬ 
diate  Nd  concentrations  have  PA  signals  which  van'  as  r~" 
where  n  is  intermediate  between  1.0  and  1.5. 

Understanding  the  results  described  above  poses  a  difficult 
problem  in  the  light  of  current  theories  for  photoacoustic 
signal  generation.  Three  parameters  are  necessary  for  the¬ 
oretical  analysis  of  PA  measurements:  the  sample  thickness 
l.  the  optical  penetration  depth  l„p,  and  the  thermal  diffusion 
length  Iq,.  For  the  samples  investigated  /  is  between  1  and 
2  mm.  The  optical  penetration  depth  is  characterized  us  l/rr 
where  or  is  the  absorption  coefficient  at  the  wavelength  of  the 
exciting  light.  For  the  laser  line  used  for  excitation,  :::  of 
the  order  of  80  mm  for  the  lightly  doped  samples  and  ap¬ 
proximately  0.7  mm  for  the  concentrated  Nd  sample*.  The 
former  (l„p  >  I)  is  the  “optically  thin"  case  while  the  latter  l 
<  I)  is  the  “optically  thick"  case. 

The  thermal  diffusion  length  is  given  by  llh  =  v' 3/ xe .  . 
where  is  the  thermal  diffusivity.  For  garnet  crystals’’ d  is 
about  5.0  X  10“2  cm2s-1,  whereas  in  NdPd)i  t,  |3  is  anisotropic 
with  the  largest  value  being  about’  1.0  x  10~2  cm2s _l.  The 
thermal  diffusivity  of  YVO«  should  be  of  the  same  magnitude 
as  that  for  the  garnets  and  pentaphosphates.  Thus  for  the 
range  of  chopping  frequencies  investigated,  Iq,  varies  from 
about  55  down  to  10  pm.  For  all  samples  lo,  «  1,  which  is  the 
"thermally  thick”  case.  Also,  Ip,  <  l„p  for  all  samples  inves¬ 
tigated. 

For  a  thermally  thick  case  with  lo,  <  l„p  as  investigated  here 
the  Rosencwaig-Gersho  (RG)8  theory  predicts  that  the  PA 
signal  intensity  will  vary  as  e~:M.  Only  if  Iq,  >  lop  should  a 
dependence  be  observed.  Obviously  there  is  a  striking  dif¬ 
ference  between  these  theoretical  predictions  and  the  data 
shown  in  Figs.  1-3  for  the  lightly  doped  samples.  A  possible 
explanation  for  this  discrepancy  is  the  difference  in  dimen¬ 
sionality  between  theory  and  experiment  The  RG  theory  is 
based  on  a  one-dimensional  model  where  the  thermal  gradient 
at  the  boundary  between  the  excited  surface  and  the  gas 
causes  uniform  heat  flow  back  into  the  cell.  In  the  experi¬ 
ments  described  here  only  a  small  area  of  the  total  sample 
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FIG  1  Photoacoustic  signal  intensity  as  a  (unction  of  chopping  frequency 
tor  Nd-doped  YVO,  crystals. 


surface  is  excited  by  the  laser  and  the  heat  is  generated  in  a 
cylindrical  volume  within  the  crystal.  Near  the  front  surface 
the  dominant  temperature  gradient  is  still  that  between  the 
sample  and  the  gas,  and  the  heat  generated  in  this  region 
diffuses  to  the  surface  and  contributes  to  the  PA  signal  as 
predicted  in  the  one-dimensional  model.  However,  in  the 
interior  region  of  the  sample  the  dominant  temperature  gra¬ 
dient  is  radially  outward  from  the  heated  cylinder.  A  three- 
dimensional  photoacoustic  theory  which  accurately  models 
the  details  of  our  experiment  has  recently  been  developed.9 
The  predictions  of  this  theory  have  been  shown  to  be  in  ex¬ 
cellent  agreement  with  our  data  with  no  adjustable  parame¬ 
ters. 

II.  DETERMINATION  OF  QUANTUM  EFFICIENCY 

The  accurate  determination  of  absolute  quantum  ef¬ 
ficiencies  of  ions  in  solids  has  been  an  experimental  problem 
of  interest  for  many  years.  It  has  been  hoped  that  PAS 
techniques  will  provide  a  method  for  doing  this  and  we  de¬ 
scribe  there  the  results  of  measuring  the  quantum  efficiencies 
of  the  neodymium-doped  crystals.  The  maximum  PA  signal 
intensity  and  the  phase  at  signal  maximum  were  recorded  for 
excitation  wavelengths  of  4765  A  and  5145  A  at  chopping 
frequencies  of  312  and  1000  Hz. 

The  PA  signal  at  phase  angle  0  can  be  described  by 
/„((!)  =  CiP„/Ea)  £  cos[^  +  tan_1(2iri'cr.)  -  U],  (1) 

‘.j 

where  P„  and  Ea  are  the  power  absorbed  and  the  energy  of  the 
level  where  absorption  occurs,  C  is  a  factor  accounting  for  the 
properties  of  the  cell  and  detection  system,  and  </»"/  is  the 
probability  for  a  nonradiative  transition  between  levels  i  and 
j  separated  by  energy  E,;  where  the  initial  state  has  a  lifetime 
r,.  The  summation  is  over  all  relaxation  processes  that  occur 
after  absorption.  4*  is  the  phase  shift  due  to  the  detection 
process. 

The  summation  in  Eq.  ( 1 )  can  be  evaluated  by  considering 
the  energy-level  scheme  and  transitions  for  Nd;,+  ions  shown 
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FIG.  2.  Photoacoustic  signal  intensity  as  a  (unction  ol  chopping  frequency 
lor  Nd3*  ions  in  garnet  crystals. 

in  Fig.  4.  For  the  lines  of  the  argon  laser  used  for  excitation, 
absorption  occurs  in  either  the  4f7  n  '2  or  ‘Ggn  levels  at  energies 
and  Eg.  This  is  followed  by  a  cascade  of  nonradiative 
decay  processes  to  the  iF-xn  metastable  state.  From  this  level 
radiative  decay  can  occur  with  a  rate  VV/<  to  the  various  4/j 
multiplets  with  branching  ratios  5r„.  Nonradiative  decay 
processes  within  the 4/  term  return  the  ion  to  the  ground  state. 
The  metastable  state  can  also  relax  through  direct  nonradi¬ 
ative  processes  with  a  rate  or  through  cross  relaxation 
transitions  with  neighboring  Nd!+  ions.  This  process  is 
generally  referred, to  as  concentration  quenching  and  occurs 
at  a  rate  W\.  The  quantum  efficiency,  probability  of  non- 
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lor  Nd,Y ,_.PsOh  crystals. 
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FIG  4.  Energy  levels  and  transitions  lor  Nd3*  ions. 

radiative  decay,  and  probability  of  concentration  quenching 
of  the  metastable  state  are  given  by 


QE  =  VV«/(  Wk  +  Wyit  +  kV.v ) . 

(2) 

P.v/t  =  IVWIH’g  +  W.vk  +  Wx), 

(3) 

P.x  =  VV'.v/llPg  +  W’.vft  +  W.  y ), 

(4) 

1 

where 

QE  +  P.v/r  +  t‘x  =  1-  <5) 

All  of  the  sample  properties  needed  for  Eqs.  ( 1  )-(5),  such 
as  energy  levels,  branching  ratios,"'  and  lifetimes  are  known 
for  Nd:i+  ions  in  YaAlsOi*,  YV04,  and  Nd,  Y|-,P.-,Ou  crystals. 
Since  only  the  *F. >/2  level  has  a  long  enough  lifetime  to  cause 
a  phase  shift  in  the  PA  signal,  the  evaluation  of  Kq.  ( 1 )  with 
the  model  of  Pig.  4  gives 

/„«/)  =  O'  (/’„/£„  i )  I  £„  r,  cos  ( ^  -  0) 

+  (fir, i  -  kQE)cos{ip  +  tan-l(2tr/',  r.J  -  ft] I,  (6) 

where 

k  —  fir, i  “  5-,4fi.| i  —  h (i  —  hr.jfi^i.  1 *1 

and  Tr,  is  the  lifetime  of  the  level. 

One  of  the  problems  in  making  absolute  photoacoustic  in¬ 
tensity  measurements  is  the  determination  of  the  factor  C 
representing  cell  and  system  response  characteristics.  To 
eliminate  this  factor  we  took  the  ratio  of  PA  signals  at  two 
different  excitation  wavelengths.  If  we  assume  that  the  re¬ 
sponse  characteristics  of  the  system  are  approximately  con¬ 
stant  at  least  over  a  limited  spectral  region,  the  unknown 
factor  cancels  out  and  upon  solving  for  the  quantum  efficiency 
gives 


QE  = 


AfiTscoslyl  -  <h>  +  dfiru  cos]yf  +  tan  M2trc,  r  -j  -  IF]  -  fi.^cosli/  -  I),-,)  -  E;,\  coslyt-  +  tan~'(2tri’,  r,-,) 
k\A  cos |-/  +  lan',(23ri’,  r -, )  -  0-\  -  cos|i^  +  tan_l(2Ti/,  r:,)  -  //,■.]] 


(8) 


where 

.4  =  t/,;,/:)tfi,i,/fi7l)(/V/,.|).  (9) 

The  phase  shift  due  to  equipment  response  ^  (roses  another 
problem.  This  was  handled  by  deriving  a  second  equation 
found  from  maximizing  Kq.  (6)  with  respect  to  0  and  again 
solving  for  quantum  efficiency: 

QE  m  i[ - -  +  fij  (10) 

k  'sin(^  +  tan_,t2trt',. rsl  —  0n\  / 

Similar  equations  exist  for  the  two  different  excitation 
wavelengths.  These  two  equations  along  with  Eq.  (8)  can  be 
solved  simultaneously  in  an  iterative  way  to  obtain  a  unique 
value  for  QE.  To  check  the  uniqueness  of  the  result,  the 
procedure  was  repeated  for  data  obtained  at  both  chopping 
frequencies. 

The  quantum  efficiencies  obtained  by  the  above  procedure 
are  listed  in  Table  1.  YjAljOi^Nd  was  the  only  one  of  the 
garnet  samples  analyzed  because  some  of  the  important  pa¬ 
rameters,  such  as  branching  ratios,  are  not  yet  established  for 
the  gallium  garnet  and  mixed  garnet  host  crystals.  A  value 
of  0.60  was  found  for  the  quantum  efficiency  of  Nd3*  in  this 
host.  This  is  much  smaller  than  the  value  of  0.91  that  is 
predicted  by  Judd-Ofelt  theoretical  calculations.11  However, 
it  is  quite  close  to  the  value  of  0.56  obtained  by  Singh  et  al. 12 
by  direct  measurement  techniques  and  the  values  of  0.48  and 
0.63  obtained  by  using  indirect  laser-pumping  techniques.1314 
The  variations  in  measured  values  may  be  due  to  samples 
having  different  Nd3+  concentrations  and  thus  different  levels 


I 

of  concentration  quenching.  Polycrystalline  powders  of 
Y;iAlr,0|j:Nd  have  been  found  to  have  a  quantum  efficiency 
close  to  1.0  but  this  is  greatly  reduced  when  single  crystals  are 
formed.1*’  This  has  been  attributed  to  crystal  defects  acting 
as  quenching  centers1*  and  thus  explains  the  apparent  dis¬ 
crepancy  between  measured  values  and  theoretical  predic¬ 
tions. 

The  results  obtained  for  the  2.0%  and  3.0%  Nd-doped  Y V04 
crystals  give  quantum  efficiencies  of  9.59  and  0.54.  respec¬ 
tively.  Calculations  by  Judd-Ofelt  theory16  predict  a  QE 

TABLE  I.  Quantum  efficiencies. _  _ 


QE 


Sample 

PAS  method 

Other  methods 

Ke  Terence 

Y;,AIsO|2 

0.91 

(10) 

0.56 

(11) 

(0.85%  Nd) 

0.60 

0.48 

1 121 

0.63 

(131 

YVCb 

. 

0.72 

(15) 

(2.0%  Nd)  . 

0.59 

0.53 

ia) 

(3.0%  Nd) 

0.54 

0.51 

la) 

Nd,  Y,_,  PjOu 

0.40 

(16) 

(100%  Nd) 

0.45 

0.50 

ll» 

150%  Nd) 

0.60 

0.66 

(b) 

( 10%  Nd) 

0.90 

1.00 

(b) 

(a)  From  lifetime  meusuremenu  and  the  radiative  lifetime  of  ( 15). 
(h)  From  lifetime  quenching  measurement*  of  (4). 
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value  of  0.72.  For  the  relatively  high  concentrations  of  Nd 
in  our  samples,  concentration  quenching  will  lower  the  mea¬ 
sured  QE  from  this  theoretical  value.  Use  of  the  calculated 
radiative  lifetime  and  our  measured  fluorescence  lifetimes 
predicts  values  of  QE  of  0.53  and  0.51  for  our  samples.  No 
direct  measurements  of  QE  for  Nd  in  YVO,  are  available  for 
comparison. 

For  NdFf,014  crystals  a  value  of  QE  of  0.45  was  obtained, 
which  is  close  to  the  value  obtained  by  Auzel  et  al. 17  through 
relative  PAS  measurements  on  powdered  samples.  The 
combination  of  this  result  with  the  measured  concentration 
quenching  of  the  fluorescence  lifetime4  predicts  a  QE  for  the 
4r  ,  ■>  level  in  the  absence  of  concentration  quenching  in  this 
system  of  0.00.  This  is  consistent  with  the  values  of  QE  found 
for  lightly  doped  Nd,  Y|_,  I’, On  as  listed  in  Table  I.  No  di¬ 
rect  measurements  of  quantum  efficiency  in  this  system  is 
available  for  comparison. 

III.  DISCUSSION  ANO  CONCLUSIONS 

The  results  presented  here  on  the  frequency  dependence 
of  the  PA  signal  intensity  on  a  variety  of  lightly  doped  and 
heavily  doped  samples  show  that  there  is  still  much  work  to 
be  done  in  developing  an  accurate  theoretical  description  of 
PA  signal  generation  under  typical  experimental  conditions. 
An  attempt  was  made  to  better  simulate  the  one-dimensional 
model  upon  which  theoretical  calculations  have  been  based 
by  using  a  lens  to  disperse  the  incident  laser  beam  over  most 
of  the  sample  surface.  However,  when  this  was  dqne  a  very 
different  frequency  dependence  of  the  signal  was  obtained 
which  was  less  than  linear  over  most  of  frequency  range, 
reached  a  maximum  at  about  40(1  Hz,  and  decreased  at  lower 
frequencies.  This  erratic  dependence  appears  to  lie  due  to 
scattered-light  ai  d  cell-surface  effects.  The  experiments  were 
again  repeated  with  no  lens  used,  resulting  In  an  excited  area 
of  about  1  mm  tn  radius.  The  results  obtained  under  these 
conditions  are  the  same  as  those  shown  in  Figs.  1-3.  Finally, 
to  ensure  that  scattered  light  was  not  affecting  the  results 
obtained  with  small-area  excitation,  the  laser  was  tuned  to 
48H0  A.  This  wavelength  of  light  is  not  absorbed  efficiently 
by  Nd:u  ions  and,  therefore,  the  sample  PAS  signal  should  be 
essentially  zero  while  the  scattered  light  signal  should  be  ap¬ 
proximately  the  same  as  for  the  5145-A  excitation.  It  was 
found  that  the  scattered  light  signal  was  approximately  one- 
and-one-half  order  of  magnitude  smaller  than  the  smallest 
sample  PAS  signals.  It  should  be  emphasized  that  the  spec¬ 
ular  reflection  troin  the  polished  sample  surface  came  back 
out  through  the  front  window  of  the  cell  and  was  measured  to 
be  approximately  the  same  for  both  excitation  wavelengths. 
These  experiments  indicate  that  experimental  artifacts  such 
as  scattered  light  are  not  affecting  the  data.  Also  fur  samples 
of  the  type  investigated  here  the  small-diameter  laser  beam 
appears  to  be  favorable  to  the  dispersed  excitation  source  and 
it  is  important  to  further  develop  theoretical  models  to  ac¬ 
count  for  these  types  of  experimental  conditions. 

Hecent  work  has  been  presented  which  accounts  for  some 
of  the  three-dimensional  aspects  of  the  heat  flow  problem.1"  18 
The  results  of  Cjuimby  and  Yen18  indicate  that  no  matter 
where  the  heat  enters  the  cell  from  the  sample,  the  one-  and 
three-dimensional  models  are  equivalent  as  long  as  cell-wall 
effects  are  not  important.  However,  for  the  experiments 
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described  in  this  paper  the  important  point  is  that  in  true 
three-dimensional  heat  How  Inside  the  sample,  some  of  the 
heat  that  is  predicted  to  reach  the  surface  in  a  one-dimen¬ 
sional  How  model  will  actually  never  reach  the  surface  within 
a  duty  cycle  of  the  experiment.  McDonald18  has  recently 
developed  a  three-dimensional  photoacouslic  theory  based 
on  a  point  source  lying  on  the  axis  of  a  cylindrical  sample. 
The  new  theory  developed  by  ChowJ  gives  a  rigorous  result 
lor  arbitrary  beam  profile.  This  theory  is  the  most  appro¬ 
priate  one  for  comparing  with  data  obtained  from  our  exper¬ 
imental  geometry  and  the  good  fit  that  is  obtained  indicates 
that  these  data  are  affected  by  the  three-dimensional  aspects 
of  the  experiment.  The  data  of  Quimliy  and  Yen1"  can  also 
be  fit  well  by  this  theory. 

Another  possibility  for  interpreting  the  results  could  be  the 
dominance  of  the  thermal  expansion  contribution  to  t  he  FA 
signal. -u  This  term  will  dominate  only  if  ,i/-  »  p,,/' /,  7’,,). 
where /f  /-  is  the  thermal  expansion  coefficient,  p,,  is  the  ther¬ 
mal  diffusion  length  in  the  gas,  I,  is  the  sample  thickness,  and 
7’, i  is  the  ambient  temperature.  For  an  air-filled  cell  at  a 
frequency  of  100  Hz.  p,.  ^  2.5  X  llr-  cm.  For  77,  =  300  K  and 
a  typical  sample  thickness  of  U.I5  cm,  the  right-hand  side  of 
the  inequality  was  found  to  be  about  5  X  K)~ 1  K-1 .  Since  the 
thermal  expansion  coefficients  of  the  type  of  materials  in¬ 
vestigated  here-1  are  of  the  order  of  7  X  ID-"  K-1,  the  in¬ 
equality  does  not  hold  and  thus  thermal  expansion  makes  a 
negligible  contribution  to  the  PA  signal. 

The  method  of  measuring  quantum  efficiencies  by  pho- 
toacoustic  techniques  appears  to  give  results  generallv  con¬ 
sistent  with  other  predictions  and  measurements.  The  pro¬ 
cedure  of  eliminating  unknown  system-response  properties 
by  obtaining  the  ratio  of  P.A  signals  at  two  different  excitation 
wavelengths  was  first  suggested  by  Hockley  and  Waugh--  lor 
determining  the  QE  of  organic  dye  molecules.  Two  other 
attempts  have  been  made  to  obtain  the  t)E  of  Nd1*  ions  in 
solids.  Auzel  vt  al. 17  studied  powdered  samples  of  concen¬ 
trated  stoichiometric  neodymium  materials  and  eliminated 
equipment  response  factors  by  forming  a  calibration  curve 
from  samples  of  known  quantum  efficiencies,  (juimhv  and 
Yen7’  investigated  Nd',+  in  a  glass  host  and  eliminated  the 
problem  of  equipment  response  factors  by  measuring  lifetime 
and  PA  signal  changes  on  a  series  of  samples  with  different 
concentrations.  Their  results  extrapolated  to  very  low  con¬ 
centrations  imply  that  the  QE  for  this  system  is  between  O.tio 
and  0.75  depending  on  the  Nd  sites  selectively  excited  by  the 
laser.  These  values  are  less  than  the  value  of  0.9  measured 
by  other  techniques  and  theoretically  predicted.-'  4  The 
chopping  frequency  used  in  this  previous  work  was  only  22  Hz 
and  this  eliminated  the  necessity  for  treating  the  phase-shift 
differences  of  the  PA  signal  generated  by  phonons  given  off 
in  different  relaxation  processes.  For  our  experimental  setup 
we  found  that  the  signal  to  noise  ratio  decreased  below  about 
100  Hz  and  thus  in  our  case  more  accurate  data  could  be  ob¬ 
tained  at  the  higher  chopping  frequencies.  Experimentally 
this  problem  of  phase-shift  differences  could  be  eliminated 
by  pumping  directly  into  the  *Fmi  level. 

The  most  significant  problem  in  the  method  of  obtaining 
QE’ s  described  here  is  the  accurate  determination  of  the  pa¬ 
rameter  A  given  in  Eq.  (9).  A  change  in  this  parameter  by  an 
amount  of  only  0.01  can  change  the  value  of  QE  by  about  0.10. 
For  the  type  of  samples  studied  the  ratio  of  PA  signal  inten- 
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sities  can  be  measured  accurately  and  the  ratio  of  the  energies 
of  the  excitation  lines  is  known  exactly.  The  problem  is  in 
accurately  determining  the  ratio  of  the  power  absorbed  at  the 
two  excitation  wavelengths.  The  major  difficulties  arise  in 
measuring  the  amount  of  light  scattered  at  the  cell  window 
and  sample  surface  and  in  determining  the  exact  absorption 
coefficients  for  sharp  absorption  lines  at  the  positions  of  very 
narrow  band  laser  lines.  Thus  the  values  of  QE  quoted  in 
Table  1  should  be  considered  to  have  an  error  range  of  as  much 
as  ±0.10. 

In  conclusion,  PAS  techniques  can  be  used  to  obtain  in¬ 
formation  on  radiationless  relaxation  processes  of  ions  in 
solids.  However,  it  is  important  to  develop-better  theoretical 
descriptions  of  the  signal  generation  process  when  narrow- 
beam-diaineter  laser  excitation  sources  are  used.  Quantum 
efficiencies  can  he  determined  even  in  complicated  cases  in¬ 
volving  numerous  radiationless  processes  if  enough  infor¬ 
mation  is  known  concerning  energy  levels  and  branching  ra¬ 
tios.  However,  absolute  measurements  are  difficult  and  it  is 
important  to  develop  methods  of  eliminating  cell-  and  sys¬ 
tem  -  respt mse  p  rope  r t  ies. 
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V.2 

EFFECTS  OF  PRESSURE  ON  THE  SPECTRA  AND  LIFETIMES  OF  RUBY 


AIjOjjCV’  results  and  interpretation 

The  change  in  the  fluorescence  lifetime  of  the  R  lines  of  the  ruby  sample  with  pressure 
is  shown  in  figure  5.  Despite  the  relatively  large  scatter  in  the  data  inherent  in  making 
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Figure  5.  Pressure  dependence  of  ihe  lifetime  of  the  R,  tluorescence  in  rubv. 


measurements  on  such  a  small  sample,  it  isclear  that  the  lifetime  increases  with  increasing 
pressure.  All  decay  curves  were  observed  to  be  single  exponentials.  It  is  not  clear  at  this 
time  whether  these  results  are  due  to  a  decrease  in  the  radiative  or  non-radiative  decay- 
rate  of  the  ;E  levels  of  Cr*'  with  increasing  pressure. 

The  detailed  fluorescence  spectrum  of  ruby  observed  at  one  atmosphere  has  been 
reported  previously  (Powell  and  Di  Bartolo  1972).  The  peaks  at  7017  A  and  7049  A  are 
the  N;  and  N,  lines,  respectively,  arising  from  exchange  coupled  pairs  of  Cr*+  ions.  The 
peaks  at  6977  A.  6985  A,  and  6996  A  also  arise  from  pairs,  whereas  the  longer-wave- 
length  peaks  correspond  to  vibronic  sidebands  of  the  Rt  transition  (Kushida  and  Kikuchi 
1967).  The  pressure  dependence  of  the  peak  positions  is  plotted  in  figure  6  with  approx¬ 
imate  rates  of  change  of  the  line  positions  with  pressure  listed  in  table  3.  It  can  be  seen 
that  Nj  is  more  sensitive  to  hydrostatic  pressure  than  N2,  which  is  consistent  with  the 
uniaxial  stress  data  of  Mollenauer  and  Schawlow  ( 1968).  This  difference  has  been 
attributed  to  the  details  of  the  interaction  between  third-nearest-neighbour  pairs  giving 
rise  to  Ni  and  fourth-nearest-neighbour  pairs  giving  rise  to  N:. 

For  vibronic  lines,  table  3  also  gives  the  pressure-induced  energy  shifts  relative  to 
the  R,  line  shift.  These  shifts  in  phonon  energiescan  be  used  to  estimate  mode  Griineisen 
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Figure  6.  Pressure  dependence  of  the  long-wavelength  side  peaks  of  the  R  fluorescence  of 
ruby. 


Table  3.  Pressure  shifts  and  Gruneisen  parameters  from  the  long-wavelength  side  peaks  of 
the  ruby  R  fluorescence. 


Positions  at  1  bar 

dv 

d7> 

tem'1  kbar"') 

d»H,  JV 

d  r  d  p 

tem'1  kbar  :) 

Y 

(A) 

V 

(cm  ') 

(cm;) 

ny77 

14333 

72 

-0.82  £  0.03 

6985 

14316 

89 

— 

6996 

14394 

111 

-0.88  £  0.03 

tN;)  '017 

14  351 

154 

-0.87  £  0.03 

iN,;  7twy 

14186 

:iy 

-1.08  £0.03 

”IJ76 

14133 

373 

-0.86  £  0.04 

0.11  £  0.04 

1.0  £  0.4 

7134 

14017 

388 

-0.86  £  0.03 

0.11  £  0.03 

0.~  £  0.2 

7 148 

139VO 

415 

-0.99  £  0.03 

0.34  £  0.03 

1.5  £  0.2 

7164 

13959 

446 

-1.14  £  0  04 

0.39  £  0.04 

2.2  £  0.2 

718b 

13916 

489 

-0.98  £  0  06 

0.33  £  0.06 

1.2  £0.3 

73(18 

13873 

532 

-O.y 6  £  0.08 

0.21  £  0.08 

1.0  £  0.4 

parameters  from  the  equation 

_  d  In  to,  _  Bj  dv 

Y,~  din  V  ^  ' 

where  DT  is  the  isothermal  bulk  modulus  which  is  about  2.57  Mbar  for  ruby  (Finger  and 
Hazen  1978).  The  mode  /s  obtained  from  equation  (7)  and  the  observed  data  are  listed 
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in  table  3.  These  are  difficult  to  compare  directly  with  other  types  of  data  due  to  the 
uncertainty  in  identifying  the  vibromc  peaks  with  particular  points  in  the  Rriliouin  zone. 
However,  the  values  obtained  here  fall  in  the  same  range  as  the  Al-O-,  mode  y's 
determined  from  the  pressure  derivatives  of  elastic  moduli  (Gerlich  l'v^Oi 


The  pressure  variations  of  the  lifetime  of  the  ruby  R  lines  and  the  positions  of  the  N 
lines  and  vibronic  peaks  demonstrate  the  variety  of  effects  that  can  be  studied  with 
pressure  techniques.  The  lifetime  change  most  probably  reflects  changes  in  the  nature 
of  the  crvstal  field.  The  different  variations  in  the  positions  ot  the  N,  and  N:  lines  retlect 
changes’ in  the  exchange  interactions  between  the  pairs  of  C  r'  ions.  The  change  in 
vibronic  peaks  can  be  attributed  to  changes  in  the  lattice  phonons  with  increased 
pressure  All  of  these  areas  warrant  further  study,  but  it  is  obvious  that  the  application 
o,  hydrostatic  pressure  can  provide  a  great  deal  of  useful  information  on  the  spectral 
dynamics  of  impurity  ions  in  solids. 
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with  elements  p^,  ?21  =  The  interaction  hamiltonian 

H 1 { t )  is  of  the  dipole  type,  and  can  be  written  as  H*  (t)  =  -  y£(t) .  Since 
the  dipole  transitions  are  between  states  cf  definite  parity,  p  =  = 

0;  and  the  phases  of  U, (r)  and  U^tr!  can  be  chosen  such  that  = 

u.  One  can  then  write  the  interaction  hamiltonian  as 


H*  (t) 


f  0  -pE (t)  | 


■uE(t)  0 


rhe  unperturbed  hamiltonian  is  given  by 


h  0 


f  I 

1°  '2J 


where  c1  is  the  energy  of  the  excited  state  and  is  the  energy  of  the 
ground  state. 


Thus  the  hamiltonian  H  is  given  by 


-UE (t) 


The  ensemble  average  <u>  of  the  dipole  moment  induced  by  E(t)  is 


given  by 


<u>  =  Tr(py) 


M (312+P21! 


U(021+i::21) 


Phe  density  matrix  satisfies 


i £ .  =  Li  rH  . ,1 

dt  h  lH’pJ 


which  becomes 
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1 


do  -i 
dt  h 


f  pElt) (0,-0*) 

I  21  21 


UE(t)  ^22-ou)  +  (e2-e1)o,1 


-uE(t)(p22-pu)-(e2-e1)o*l! 


.  (A-5) 


;E!tl(W' 


Defining  the  resonant  frequency  u) 


o  f> 


C2~C1 

,  it  follows  from  (A-5) 


that 


dp. 


dt 


dp 


22 


dt 


-“oC21  *  f  E<tl  l»lfC22l 


ill  (p  _p*  ) 

fi  P21  ^2l' 


.  A— 6 ) 


(A- 7) 


and  using  the  normalization  condition, 


_d_ 

dt 


(pll"°22) 


2iu 

h 


E(t)(p  -p*  )  . 

21  21 


(A- 3 ) 


The  above  equations  of  motion  for  the  density  matrix  do  not  include 
collisional  effects.  When  the  perturbing  field  E(t)  is  turned  off,  c.ne 
expects  the  off  diagonal  terms  in  the  density  matrix  to  vanish  as  the 
relative  phase  coherence  among  the  eigenfunctions  of  the  ensemble  is  lost 
via  collisions.  These  collisions  will  conserve  the  average  energy  of 
each  level,  but  cause  a  loss  of  information  about  the  phases  of  the 
vavefunctions .  The  jth  diagonal  element  of  the  density  matrix  represents 
the  fraction  of  the  systems  in  the  ensemble  that  will  give  the  answer 
Ej  when  the  energy  is  measured.  If  one  takes  into  account  only  collis¬ 
ions,  these  will  be  given  by  the  Boltzmann  factors,  and  the  equilibrium 
l 

density  matrix  C  is 
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i  <■ 


“VV 


0 


0 


'I 

I 


/k. 


_ei/kbT  "VV 

where  Q  ■  partition  function  =  l  +  l 

If  the  perturbing  field  is  turned  off,  the  density  ma trices  will  de- 
cay  back  to  equilibrium  with  decay  constant  t  ^  as 


(P3k  1  P3k] 


Because  of  the  normalization  condition,  for  a  two-level  system  the 
relaxation  times  of  the  diagonal  elements  must  be  equal  and  is  called  the 
longitudinal  relaxation  time  t^.  Since  *  2*^,  the  off  diagonal 

relaxation  times  are  also  equal  and  is  called  the  transverse  relaxation 
time  t?.  The  transverse  relaxation  time  can  be  thought  of  as  a  phase 
coherence  relaxation  time. 

After  including  collisional  effects,  (A-6)  and  (A-3)  become 


dD 


21 


dt 


"  iwoP21  +  T  (0u"P22)E(t-  " 


(A-9) 


_d_ 

dt 


(Pll'P22) 


2iUE ( t) 


(°2rp2i) 


^Pll”°22^ 


t 


1 


CA-10) 


If  one  assumes  that  the  local  perturbing  field  E(t)  is  harmonic, 

E ( t)  =  E^costut;  then  if  w  ~  ouq  the  slowly  varying  variable  •naV'  be 

defined  as 

°21(t)  "  021(t)J"1Ut  •  <A"U: 
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Substituting  into  (A-9) ,  and  throwing  out  the  nonsynchronous  term 
which  averages  to  zero,  Equation  (A-9)  becomes 


da 


21 


dt 


i(ui-<jj  )a_,  + 
o  21 


iUE 

_ c 

2  ft 


(Pir°22)  - 


21 


( A  - 1 2 ) 


terms 


Substituting 


,  2iut 
x.  and 


(A-ll)  into  (A-10)  and  throwing  out  the  nonsynchronous 
-2iut  ... 

which  average  to  zero,  one  oDtains 


_d_ 

dt 


W 


iUE 


VL1  - 


vy 


to  -o  ) 
11  22 


(A-13) 


To  obtain  the  steady  state  solutions  to  the  density  matrix,  the 
left  hand  side  of  (A— 12 )  and  (A-13)  is  set  equal  to  zero.  The  preces¬ 
sion  frequency  Q  is  defined  as 


HE 
_ o 

2ft 


and  the  steady-state  solutions  to  ( A— 12 )  and  (A-13)  are 


Se(a2i)  = 


-  fit"  (P  -p„  )  ""  (w-U  ) 

_ 2  11  i2 _ 2_ 

1  +  457*'  t,t-  +  )' 

12  2  o 


(A-14) 


Im(a?1 ) 


"t2 !?11~°22^ 


1  +  t,  t_  -  t  ) 

12  2  o 


(A-15) 


011~°22 


■  <0ll"D22)  ( 


,  ,  .  2  2 
1  +  (U>-U  )  t 
o  2 


1  +  (u-'x)  )  +  4Q'ct,t, 

0-2  2  1 


(A-16) 


If  N  is  the  density  of  active  atoms  or  ions,  then  IN  »  N (0^ -o7„)  is 
the  average  density  of  the  population  difference  between  the  two  levels. 


no 


ioi 


and  AN  »  Nlc^-c^.,)  is  the  population  difference  oetween  the  two  levels 
at  2ero  field.  From  ( A- 16 ) , 


2  2 

3  1  +  (w-u  )  t 

AN  -  ^ - 1_ - }  _ 

1  +•  (uj-u  )  t„  +  4H  t  t, 

O  2  2  1 


The  macroscopic  polarization  P  is  given  by 


P  =  N < ii '■  or,  using  (A-3)  , 


?  =  ^^21  +  J2l'  which  can  be  expressed  in  terms  of  as 


P  *  2Ny  (Re  (a^)  cosut  +  Im  < CT21  ^  sinwt) 


Substituting  in  for  Re(a^^)  and  Im((T  )  ,  one  may  write  the  polariza- 


2  i 

y  An  t  sinut  +  (u  — w)  t„ccsut 

- 2  E  t - 2 - iu _ ; 

b  o  ,  ,  .  2  2  _2 

1  +  ( w— -u  /  t  +  4..  t,t. 

O  2.  I  ~ 


;a-17) 


One  may  now  express  the  polarization  in  terms  of  the  susceptibility 


P(t)  =*  <z  (Re (X) )  E  cosojt  -  e  (Im(X))E  sinut  ; 
o  o  o  o 


comparing  with  (A- 17) ,  it  can  be  seen  from  observation  that 


Re('X) 


2  A  l 

u  t  £n 


<to  -u)  t„ 

■ _ o  2 _ i 

•>  n  ' 

1  +  (u-Mj  )‘t‘  +  41'/ 1,  t 
o  2  12 


lA-13) 
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